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The  majority  of  data  presented  in  this  work  is  for  a  helicon  plasma  discharge  driven  at 
13.56  MHz,  500  Watts  input  power,  900  Gauss  applied  magnetic  field,  10  mTorr  neutral 
Argon  gas,  and  cylindrical  plasma  of  5  cm  diameter  approximately  50  cm  long. 

High  frequency  magnetic  induction  probes  were  developed  to  measure  helieon  wave 
propagation  using  a  new  teehnique  for  frequeney  calibration  through  an  impedance 
analyzer;  up  to  100  MHz.  This  work  demonstrates  magnetic  field  measurements  in  high 
frequency  plasma  are  greatly  simplified  through  this  new  frequency  characterization 
method.  Line-lengths  and  transmission-cable-types  are  readily  identified  as  diagnostie 
limiting  factors.  The  magnetic  probe  design  enables  the  first  3-dimensional  imaging  of 
plasma  waves  through  detailed  radial  and  axial  measurements.  Strong  agreement  is 
obtained  between  the  measured  hr,  be,  and  bz  radial  profiles  with  the  numerieal  solutions 
of  helieon  waves  when  a  non-uniform  radial  density  profile  is  considered.  The  axial 
helicon  wavelength  predicted  by  the  non-uniform  radial  density  theory  also  agrees  with 
the  measured  wavelength  when  the  full  three  dimensional  wave  is  aeeurately  analyzed. 

In  some  cases,  the  differences  between  the  three  dimensional  wavelength  and  the 
numerieally  solved  values  are  less  than  30%.  This  is  in  eontrast  to  the  two  dimensional 
wavelengths  which  can  differ  from  the  numerieal  values  by  greater  than  100%.  We  show 
a  complete  visual  representation  of  helicon  waves  through  3-d  imaging  whieh  provides 
signifieantly  more  accurate  analysis  of  the  helicon  wavelength. 
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This  work  also  observed  a  density  peak  downstream  of  the  antenna  /  source  through  axial 
density  measurements  with  a  RF  compensated  Langmuir  probe  (calibrated  against  a  90 
GHz  microwave  interferometer).  Here,  the  downstream  density  peak  is  explained  in 
terms  of  a  global  energy  balance  modeled  by  an  axially  decaying  electron  temperature 
peaked  at  the  source;  Te  ->  3  -  7  eV.  This  model  does  not  require  an  assumption  of  a  RF 
plasma  or  additional  heating  by  helicon  wave  absorption;  rather  the  model  demonstrates 
excellent  agreement  with  the  measured  axial  density  profiles  when  radial  losses  are 
assumed  to  be  less  than  1 0%,  which  is  reasonably  attributed  to  the  axially  applied  static 
magnetic  field. 

The  new  diagnostic  methods  developed  in  this  work  (both  probe  characterization  and 
analysis  of  measured  data)  provide  fundamental  insight  into  laboratory  helicon  plasma. 
The  methods  and  results  derived  here  will  supplement  and  aid  the  design  of  future 
helicon  plasma  sources. 
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CHAPTER  1:  INTRODUCTION 


The  primary  objective  of  this  research  will  be  to  accurately  quantify  the  helicon 
wavelength.  This  will  be  accomplished  through  a  detailed  and  innovative  method  toward 
magnetic  probe  design  used  to  measure  the  helicon  plasma  wave  fields.  The  purpose  is  to 
provide  thorough  analysis  of  helicon  waves  through  three  dimensional  imaging  in  support 
of  future  design  and  understanding  of  helicon  plasma  sources. 

The  term  ‘helicon’  was  first  used  in  1960  by  Aigrain  [1],  to  describe  waves  that 
propagate  in  the  presence  of  a  magnetic  field  for  the  frequencies  between  the  ion  and 
electron  cyclotron  frequencies.  More  generally,  helicon  waves  are  a  subclass  of  whistler 
waves  in  the  sense  helicons  propagate  in  a  bounded  system  [2],  specific  to  this  study,  a 
cylindrical  boundary.  There  exist  other  complexities  to  helicon  waves  that  separate  them 
from  the  theory  of  whistler  waves  and  these  differences  will  be  evident  in  Chapter  2. 

This  chapter  will  introduce  helicon  waves  and  discuss  the  evolution  of  studying  helicon 
waves  since  the  1960’s.  Specifically,  this  chapter  will  highlight  the  simplicity  in 
generating  helicon  plasma  but  the  difficulty  in  characterizing  the  plasma  and 
understanding  the  results.  A  layout  for  the  remainder  of  the  dissertation  is  also  provided. 

1.1  Whistler  Waves 

Whistler  waves  were  first  investigated  by  Storey  in  the  1950’s  [3].  Storey  performed 
measurements  on  ‘whistling  atmospherics’  which  were  described  as  waves  which  follow 
“the  lines  of  force  of  the  earth’s  magnetic  field”  and  were  recorded  by  audible  “clicks”  or 
a  descending  “whistling”  tone.  The  theory  predicts  right-handed  circularly  polarized 
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waves  that  propagate  in  free  space  where  ion  effects  and  collisions  are  neglected. 
Whistler  waves  are  typically  referenced  when  helicons  are  discussed  due  to  the  similarity 
in  deriving  both  dispersion  relations. 

1.2  Helicon  History 

Helicon  plasma  is  generated  by  driving  an  antenna  at  a  frequency  between  the  electron 
and  ion  cyclotron  frequencies.  The  antenna  is  usually  wrapped  around  an  insulating 
cylinder.  As  the  power  is  increased,  the  antenna  current  induces  a  magnetic  field  which 
drives  an  electric  field  within  the  cavity  according  to  Faraday’s  law  of  induction.  The 
power  to  the  antenna  is  continually  increased  until  the  electric  field  creates  a  breakdown 
of  the  neutral  gas.  This  is  considered  inductively  generated  plasma.  To  drive  helicon 
waves,  a  DC  magnetic  field  surrounds  the  cylindrical  cavity  altering  the  properties  of  the 
plasma  dielectric  medium.  Typically,  this  is  accomplished  by  increasing  the  plasma 
density  (though  input  power)  or  increasing  the  applied  magnetic  field.  Helicon  waves 
propagate  when  these  quantities  reach  appropriate  levels  (determined  by  the  dispersion 
relation).  Typical  laboratory  parameters  are:  pressure  (mTorr),  driving  frequency  (MHz), 
input  power  (lOO’s  Watts),  DC  magnetic  field  (lOO’s  Gauss),  tube  diameter  (cm),  and 
tube  length  (lO’s  cm). 

Early  research  on  helicons  primarily  consisted  of  characterizing  the  wave  generated 
plasma  source  without  a  specific  application  for  the  device  [4-7].  It  was  not  until  the 
1970’s  and  1980’s  when  groups  in  Australia  [6,  8-12]  and  at  the  University  of  California 
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Los  Angeles  (UCLA)  [2,  13-31]  began  extensive  investigation  into  helieon  waves  for 
plasma  proeessing. 

One  benefit  of  helicons  for  plasma  processing  lies  in  the  efficient  generation  of  high 
densities  (>  10  m'  )  for  relatively  low  input  powers  (0.1-1  kW)  while  maintaining  low 
plasma  potentials  (<  100  Volts)  [12],  Additionally,  the  ability  to  generate  uniform 
density  plasma  over  large  areas  is  attractive  for  etching  onto  large  wafers  or  substrates 
[19]. 

More  recently,  helicons  have  been  considered  for  plasma  space  propulsion  [22,  32-46]. 
They  have  been  considered  either  as  an  ionization  source  with  a  secondary  stage  for 
acceleration  or  as  a  stand-alone  propulsion  device.  Concepts  such  as  VASIMR  [33,  35, 
42]  utilize  helicons  as  an  ionization  source  and  worry  about  magnetic  field  detachment 
separately.  Similarly,  annular  geometry  helicon  plasma  is  being  investigated  as  a  primary 
stage  for  space  propulsion  [45,  47]  and  will  have  to  address  detachment  and  /  or 
acceleration  scenarios  following  the  source.  Conversely,  other  devices  such  as  the  High 
Power  Helicon  (HPH)  [38]  at  the  University  of  Washington  or  the  mini  Helicon  Thruster 
experiment  (mHTX)  [46]  at  MIT  are  utilizing  helicons  as  stand-alone  propulsion  devices 
however,  these  experiments  do  not  clearly  address  the  mechanisms  for  particle 
detachment  or  thrust. 

A  proposed  method  to  circumvent  the  problem  of  detachment  is  utilizing  double  layer 
formation  [36,  37,  39,  44,  48]  in  helicon  plasma.  The  claim  for  this  acceleration 
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mechanism  is  that  a  ‘sharp’  potential  drop  occurs  downstream  at  lower  pressures  (<  1 
mTorr)  in  helicon  plasma,  which  can  accelerate  ions  to  supersonic  speeds.  This  reported 
mechanism  still  has  many  questions  as  to  what  causes  the  double  layer  formation  (helicon 
plasma/geometry  discontinuities/expansion  along  B-field  lines).  Additionally,  generating 
helicon  plasma  at  such  low  pressures  may  prove  more  difficult  in  terms  of  power  or  B- 
field  requirements  since  the  medium  will  have  a  relatively  lower  plasma  density  than  at 
higher  pressure  (mTorr  vs.  sub-mTorr).  In  any  event,  it  is  the  author’s  belief  that 
adapting  helicons  to  space  propulsion  stands  a  greater  chance  of  success  when  utilized 
solely  as  an  ionization  source.  The  mechanism  of  laboratory  helicon  double  layers  has 
yet  to  be  explained  whereas  the  use  of  helicons  as  a  high-density  ionization  source  is  well 
documented. 

1.3  Helicon  Research 

Helicon  waves  were  initially  studied  in  metals  and  semiconductors  by  Libchaber  [49]  and 
Rose  [50]  in  the  early  1960’s.  They  observed  wave  propagation  at  frequencies  much  less 
than  the  electron  cyclotron  frequency  with  wave  vector  k  parallel  to  the  applied  DC 
magnetic  field.  However,  it  was  Lehane  and  Thoneman  [4]  who  performed  the  first 
comprehensive  study  of  helicon  waves  in  gaseous  Xenon  where  neutral  gas  pressures 
ranged  10-70  mTorr.  The  experiment  consisted  of  a  cylindrical  Pyrex  tube;  10  cm 
diameter  -  100  cm  length.  Power  was  fed  to  an  exciting  coil  by  a  3  kW  generator  at 
frequencies  variable  from  15-17  MHz.  An  axial  Bo-field  of  up  to  500  Gauss  was 
applied  over  the  middle  50  cm  of  the  discharge  cylinder.  A  goal  of  the  work  was  to 
compare  experimental  results  with  the  theory  of  plasma  wave  propagation  in  an  insulated 


4 


cylinder  for  a  uniform  radial  density  profile  originally  derived  by  Legendy  in  1964  [51, 
52]  and  Klozenberg,  McNamara,  and  Thoneman  in  1965  [53]  (KMT  theory).  Good 
agreement  was  shown  with  wave  field  measurements  through  b-dot  probes  seanned 
aeross  the  diameter  of  the  eylinder.  Peak  densities  of  3-15x10  m'  were  obtained 
with  a  double  Langmuir  probe  as  the  magnetie  field  was  inereased  from  70-510  Gauss. 
Results  at  the  higher  magnetic  fields  did  not  agree  as  well  with  the  lower  field  tests; 
however,  the  radial  density  profiles  also  indieated  a  greater  non-uniformity  whieh  eould 
aeeount  for  the  diserepaneies  [4]. 

Kuckes  (1968)  [5]  and  Nyack  (1975)  [7]  also  observed  good  agreement  with  uniform 
radial  density  theory  observing  right  handed  eireularly  polarized  waves.  However,  it  was 
Nyaek  and  Christiansen  [7]  in  1974  who  first  considered  a  non-uniform  radial  density 
profile.  They  eonsidered  a  “top-hat”  (step  funetion)  and  parabolie  profile  whieh  were 
solved  numerically  and  compared  with  experimental  results  for  the  wave  fields.  They 
found  good  agreement  with  results  provided  the  density  profile  used  was  similar  to  the 
one  measured  experimentally. 

Some  of  the  more  extensive  measurements  on  helicons  consisted  of  the  work  by  Boswell 
[6,  8-11]  at  Flinders  University  in  Australia  and  Chen  [2,  13-31]  at  UCLA.  Both  groups 
have  thoroughly  eharacterized  a  range  of  sources  with  variables  ineluding;  antenna 
eonfiguration,  souree  diameter,  input  power,  matehing  network  configuration,  applied  Bo 
field,  driving  frequeney,  and  neutral  gas  pressure.  The  primary  diagnosties  used  by  both 
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groups  have  been  b-dot  probes  and  RF  compensated  Langmuir  probes  with  the  intent  of 
measuring  the  wave  fields,  plasma  density,  and  electron  temperature. 

Beginning  in  1970  [6]  and  later  in  1984  [8],  Boswell  performed  plasma  wave  field 
measurements  on  a  5  cm  and  10  cm  diameter  Pyrex  tube;  although  measurements  on  the 
5  cm  tube  were  “highly  distorted  due  to  the  plasma  being  significantly  perturbed  by  the 
probes”.  Experiment  conditions  were  driving  frequency  of  8.8  MHz,  pressure  1 .5 
mXorr,  applied  DC  magnetic  field  up  to  1500  Gauss.  Radial  hr,  be,  and  bz  measurements 
were  taken  and  agreement  with  the  theoretical  profiles  was  shown.  Axial  hr  and  bz 
measurements  down  the  centerline  of  the  tube  were  also  taken  and  the  axial  wavelength 
was  reported  by  phase  changes  of  180  degrees  being  observed  at  successive  minima  in 
wave  amplitude.  During  density  measurement,  Boswell  also  observed  density  jumps  at 
distinct  but  not  all  operating  conditions  as  the  antenna  length  and  power  were  altered. 

The  observed  density  jumps  will  be  discussed  later  in  this  work  in  Chapter  5. 

In  1987,  Boswell  published  further  density  ‘n’  vs.  Bq  magnetic  field  measurements  for 
the  Waves  On  Magnetized  Beams  And  Turbulence  (WOMBAT)  experiment  again 
showing  density  ‘jumps’  at  some  operating  conditions  [10].  This  time  the  experiments 
were  run  on  a  20  cm  diameter  tube  with  antenna  driven  at  7  MHz,  pressure  30  mXorr,  and 
power  up  to  1.5  kW.  Later  in  1989,  Boswell  performed  similar  density  measurements  on 
a  6  cm  diameter,  13.56  MHz,  1200  Gauss  source  tube  [11].  The  point  being  conveyed  for 
work  on  helicon  plasma  is  the  large  variety  of  source  variables  which  can  often  make 
comparisons  between  different  research  groups  difficult.  Most  groups  have  now  settled 
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on  operation  at  13.56  MHz;  which  will  be  the  frequency  used  throughout  the  current 
study. 

Beginning  in  the  early  1990’s,  UCLA  and  Chen  began  publishing  work  on  their  helicon 
efforts.  In  1991,  Chen  provided  a  detailed  derivation  for  wave  field  solutions  due  to 
uniform  radial  density  profiles  in  terms  of  Bessel  function  solutions  [2].  Although  this 
work  had  been  previously  considered,  a  comprehensive  analysis  for  experimental  design 
aspects  was  also  considered.  It  was  also  in  this  work  that  Chen  introduced  a  thorough 
analysis  including  a  collisional  damping  term  in  Ohm’s  law  when  solving  for  the  plasma 
wave  field.  Collisional  effects  gave  rise  to  a  second  branch  in  the  theoretical  derivation 
which  was  deemed  an  electrostatic  wave;  in  helicon  literature,  now  commonly  referred  to 
as  the  radial  Trivelpiece-Gould  (TG)  wave. 

In  1993,  Chen  thoroughly  derived  numerical  solutions  for  non-uniform  radial  density 
profiles  based  off  of  parabolic  distributions  [15]  for  obtaining  the  hr,  be,  and  bz  wave 
fields.  Although  this  work  was  first  considered  by  Shoji  [54]  in  1993,  the  results  were 
not  as  comprehensive,  though  Shoji  did  provide  similar  numerical  results  and  describes 
the  electrostatic  wave  branch  in  addition  to  the  helicon  branch.  Shoji  also  discusses  the 
merging  of  the  two  branches  at  low  frequency  or  low  magnetic  field  as  the  parallel  and 
perpendicular  wave  numbers  become  equal.  Further  numerical  solutions  and  wave 
absorption  effects  were  considered  by  Kramer  in  1999  and  2007  [55-57]  with  similar 
results.  The  TG  branch  will  be  described  in  further  detail  in  Chapter  2. 
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In  1994,  Sudit  and  Chen  published  RF  compensation  techniques  for  making  density  and 
temperature  measurements  in  RF  plasma  [58].  Their  technique  is  similar  to  earlier  work 
by  Godyak  [59,  60]  where  sheath  impedances  are  considered  with  respect  to  probe 
impedances  for  filtering  or  removing  unwanted  RF  pickup.  Previous  density  and 
temperature  measurements  on  helicons  had  not  thoroughly  described  whether 
compensation  had  been  considered  and  the  impact  on  quantitative  results. 

ULCA’s  experiment  consisted  of  a  2.5  cm  and  5.0  cm  diameter  tube  with  an  antenna 
driven  at  27.12  MHz.  Pressures  considered  were  approximately  10  -  20  mTorr,  powers 
up  to  2  kW,  and  applied  magnetic  fields  from  0  -  900  Gauss.  In  the  mid  1990’s  UCLA 
performed  extensive  investigation  into  radial  and  axial  density,  temperature,  and  wave 
field  measurements.  Axial  measurements  to  determine  the  helicon  wavelength  were 
based  on  single  axis  bz  scans  down  the  discharge  centerline  using  magnetic  induction  (b- 
dot)  probes.  Wavelengths  on  the  order  of  lO’s  cm  (or  approximately  the  length  of  the 
antenna)  were  reported  based  off  phase  measurements  where  local  wavelengths  agreed 
with  the  dispersion  relation  [29,  31].  The  work  of  UCLA  in  1995  was  also  one  of  the 
first  to  report  on  a  “downstream  density  peak”  up  to  several  wavelengths  away  from  the 
antenna  (source).  This  phenomenon  has  been  explained  as  a  consequence  of  decay  in 
electron  temperature  away  from  the  source  [31].  At  this  location  “the  ionization  rate  is 
low  and  the  diffusion  rate  is  large,  so  that  production  and  loss  mechanisms  would  not  be 
expected  to  produce  a  peak”.  However,  performing  an  energy  balance,  they  were  able  to 
show  a  density  rise  due  to  the  temperature  decay.  Additionally,  since  nearly  all  power  is 
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absorbed  in  the  near  field  antenna  region  of  the  helieon  diseharge  the  density  peak  was 
not  due  to  any  signifieant  wave  absorption  downstream  [31]. 

Transition  between  eapaeitive  -  induetive  -  helieon  modes  (E-H-W  modes  respeetively) 
were  investigated  in  1996  by  Ellingboe  and  Boswell  [61]  whieh  were  deseribed  by 
differences  between  the  radial  wave  profiles  and  radial  density  profiles.  Radial  density 
profiles  of  helicon  plasma  are  more  centrally  peaked  whereas  capacitive  and  inductive 
plasma  exhibit  more  uniform  profiles. 

Work  on  helicon  plasma  prior  to  2000  had  suggested  the  high  efficiency  and  absorption 
of  the  discharge  was  due  to  Eandau  damping,  however,  no  significant  population  of  high 
energy  electrons  had  been  observed  and  this  effect  has  since  been  ruled  out  [21]. 
hollowing  2000,  some  of  the  more  basic  research  has  been  on  identifying  the  existence  of 
TG  waves  (predicted  by  including  electron  mass  effects)  in  the  region  of  the  antenna. 
Inclusion  of  these  effects  reveals  the  electrostatic  radial  wave;  considered  radial  because 
the  perpendicular  wave  number  is  approximately  equal  to  the  total  wave  number  k. 
The  first  experimental  evidence  of  TG  waves  was  reported  by  Blackwell  in  2002  [26] 
using  a  J-dot  current  probe  (miniature  Rogowski  coil).  Success  in  observing  helicon  TG 
waves  was  observed  for  an  m  =  0  antenna  at  low  (<  100  Gauss)  magnetic  fields.  The 
lower  magnetic  fields  ‘expand’  the  TG  wavelength  by  lowering  the  perpendicular 
wavenumber  whereas  at  higher  magnetic  fields  the  TG  waves  have  been  too  localized 
(sub-cm)  to  accurately  measure. 
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The  vast  majority  of  theoretical  work  has  neglected  ion  effects;  the  ions  being  massive 
and  slow  allow  this  to  be  a  reasonable  assumption.  However,  some  groups  have  recently 
considered  this  effect  to  be  non-negligible  due  to  discharge  operation  near  the  lower 
hybrid  frequency  where  it  is  possible  for  the  ions  to  exchange  energy  with  the  wave  [62- 
66].  The  impact  of  ion  effects  on  the  helicon  wave  profiles  is  not  yet  known,  especially 
considering  how  strongly  the  theoretical  consideration  of  non-uniform  radial  density 
profiles  agrees  with  experimental  results.  Additionally,  some  groups  have  considered 
radially  localized  helicon  modes  (RLH  modes)  as  opposed  to  TG  waves  to  explain  the 
high  power  absorption  within  helicon  plasma.  However,  experimental  evidence  of  this 
has  been  limited  [67,  68]  and  the  implication  of  this  theory  is  not  yet  understood. 

1.4  Relevant  Previous  Work 

Previous  research  relevant  to  analyzing  the  helicon  wavelength  and  plasma  discharge 
focus  on  two  primary  aspects:  measurements  of  the  bz-field  on  the  z-axis  centerline  and 
understanding  the  downstream  density  peak  away  from  the  antenna  /  source.  Much  work 
to  date  has  taken  an  axial  location  (typically  underneath  an  end  ring  of  the  antenna)  and 
performed  radial  scans  here.  It  has  been  fortuitous  that  this  location  for  radial  scans  was 
selected  because  helicon  waves  are  truly  three-dimensional  helices  (this  will  be  shown 
later).  Consequently,  the  radial  wave  profiles  will  ‘turn  over’  (rotate  in  phase)  as  the 
axial  distance  away  from  the  source  is  traversed.  Additionally,  axial  measurements  to 
date  have  been  performed  along  the  centerline  of  the  plasma  discharge  for  quantifying  the 
helicon  wavelength.  This  method  will  be  shown  to  be  inaccurate  because  valuable 
information  regarding  the  axial  wavelength  is  lost  in  the  radial  direction.  Inclusion  of  the 
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radial  helix  component  which  contributes  to  the  axial  wavelength  may  assist  in  the 
explanation  of  why  theoretical  predictions  apparently  overestimate  the  helicon 
wavelength.  A  discussion  of  this  overestimation  was  provided  by  Franck  in  2005  [69] 
where  helicon  wavelengths  were  identified  based  on  bz  scans  for  a  10  cm  diameter  tube 
as  shown  in  Figure  1.1  along  with  the  axial  density  profile. 


axial  distance  to  antenna  (mm) 

Figure  1.1.  Axial  and  density  measurements  from  Franck  2005  Ref.  [69].  Plots  (a) 
and  (b)  are  for  different  Bo  field  strengths  and  the  line-types  correspond  to  the  density 

Boswell  showed  similar  measurements  in  1984  [8]  on  a  10  cm  diameter  tube  for  hr  and  bz 
with  axial  scans  taken  down  the  discharge  centerline  shown  in  Figure  1 .2  while  Chen 
performed  a  similar  measurement  in  1995  [16,  17]  on  a  5  cm  diameter  tube  with  axial  bz 
profiles  shown  in  Figure  1.3.  All  results  were  much  less  than  the  predicted  values  as 
calculated  from  the  dispersion  relation  and  boundary  conditions. 
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Figure  1.2.  Axial  br  and  bj  measurements  from  Boswell  1984  Ref.  [8]. 


Figure  1.3.  Axial  b^  measurements  from  Cben  1995-1996  Refs.  [15,17]. 

An  additional  aspect  of  this  work  will  focus  on  the  downstream  density  peak  observed  by 
Sudit  and  Chen  in  1995  [16,  17,  31].  This  is  shown  in  Figure  1.4  for  the  same  discharge 
described  above.  Another  group  to  have  observed  this  effect  is  Tysk  [70,  71]  from  the 
University  of  Wisconsin  as  shown  Figure  1.5  but  little  is  mentioned  regarding  the 
observation.  Perusing  much  helicon  literature  reveals  two  popular  presentations  of  data: 


12 


power  vs.  density  plots  and  radial  density  or  radial  wave  profiles.  The  axial  data  reported 
is  surprisingly  limited  to  the  few  groups  summarized  here.  It  will  be  the  primary  foeus  of 
this  work  to  investigate  the  axial  variation  in  plasma  wave  profiles  and  axial  density 
profiles.  Ultimately,  three  dimensional  mapping  of  the  plasma  waves  will  reveal  new 
insight  into  helicons.  The  axial  measurements  will  establish  a  new  definition  for  the 
helicon  wavelength. 


0  20  40  60  80 

Distance  from  Center  of  Antenna  (cm) 


Figure  1.4.  Axial  density  profile  illustrating  downstream  density  peak  from  Chen  1995-1996 
Refs.  [16-17].  The  antenna  end  rings  are  shown  hy  the  shaded  lines  straddling  the  z=0  location. 
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Figure  1.5.  Axial  density  profile  with  downstream  density  peak  from  Tysk  2003-2004 
Ref.  [71].  Density  scale  is  logarithmic. 


1.5  Objective  of  This  Work 

This  work  was  performed  with  the  following  quote  from  Stephen  Hawking’s  “A  Brief 
History  of  Time”  kept  elose  in  mind  throughout: 

“A  theory  is  a  good  theory  if  it  satisfies  two  requirements.  It  must  accurately  describe  a 
large  class  of  observations  on  the  basis  of  a  model  that  contains  only  a  few  arbitrary 
elements,  and  it  must  make  definite  predictions  about  the  results  of  future  observations.  ” 

Therefore,  we  begin  by  applying  this  idea  to  helieon  plasma  waves.  The  model  that  will 
be  utilized  throughout  this  work  is  for  a  non-uniform  radial  density  distribution  where 
eleetron  and  ion  mass  effeets  ean  be  ignored  to  prediet  the  helieon  radial  wave  profiles; 
ions  assumed  infinitely  massive  with  respeet  to  eleetrons  and  the  eleetron  eyelotron 
frequeney  is  mueh  higher  than  the  driving  frequeney.  The  arbitrary  elements  included  in 
this  model  are  driving  frequency/  static  magnetic  field  Bq,  tube  radius  a,  operating  gas 
(Argon  here),  and  radial  density  profiles  n(r)  (Gaussian  profiles  are  used  in  this  work). 
Finally,  the  definite  predictions  this  work  will  provide  include:  helicon  wavelength 
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predictions  based  off  of  three  dimensional  imaging  of  the  hr,  be,  and  bz  profdes  and  the 
ability  to  predict  axial  density  profiles  provided  an  input  temperature  profile.  This  work 
will  also  demonstrate  some  important  features  of  diagnostic  implementation  in  short  time 
scale  environments;  not  just  limited  to  radio  frequency  plasma  environments  but 
inclusive  of  magnetic  fields  due  to  fast-pulsed  currents  i.e.,  ringing  discharges. 

The  format  of  this  dissertation  is  as  follows: 

Chapter  2  analyzes  the  theory  of  helicon  waves  beginning  from  a  simplified  uniform 
density  profile  to  a  theory  that  predicts  TG  waves  to  explain  the  high  absorption  of 
helicon  plasma,  and  finally  where  non-uniform  radial  density  distributions  predict  radial 
helicon  wave  profiles. 

Chapter  3  details  the  facility  and  apparatus  used  in  this  work. 

Chapter  4  thoroughly  characterizes  magnetic  induction  probes  for  use  in  high  frequency 
plasma  and  describes  the  careful  implementation  of  using  these  probes  in  plasma. 
Additionally,  RF  compensated  Langmuir  probes  are  examined  with  a  full  frequency 
spectrum  characterization  of  the  probe  impedance. 

Chapter  5  illustrates  how  strongly  the  numerically  solved  wave  profiles  agree  with  those 
measured  by  b-dot  probes.  The  initial  radial  results  are  expanded  into  a  three 
dimensional  contour  image  of  helicon  waves  providing  visual  representation  of  a  helicon 
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wave.  The  results  will  also  identify  the  downstream  density  peak  and  provide  a 
quantitative  model  based  on  an  exponentially  decaying  electron  temperature  profile 
peaking  at  the  antenna  /  source  ‘exit’. 

Chapter  6  will  identify  the  successes  and  contributions  reported  in  this  study  as  well  as 
the  areas  where  this  work  (both  observational  and  predictive)  may  exhibit  shortcomings. 
This  will  provide  a  basis  for  where  additional  investigation  is  warranted,  specifically 
obtaining  better  accuracy  in  determining  the  electron  temperature  and  areas  where  the 
model  which  predicts  axial  density  can  be  improved  and  further  validated. 
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CHAPTER  2:  THEORY  OF  HELICON  WAVES 


The  theory  of  helieon  waves  has  steadily  evolved  over  the  past  several  deeades.  The 
foundation  is  rooted  in  the  theory  of  plasma  waves  in  an  infinite  uniform  magnetie  field 
whieh  results  in  the  eold  plasma  dispersion  relation  (epdr).  This  simple  theory  ean 
quiekly  beeome  mathematieally  eomplex  and  so  the  results  manifest  themselves  in  the 
form  of  a  Clemmow-Mullaly- Allis  (CMA)  diagram.  Sinee  helieon  waves  are  bounded 
whistler  waves  whieh  otherwise  propagate  in  free  spaee,  the  CMA  diagram  provides  a 
good  starting  point  for  analysis.  However,  for  the  aetual  bounded  plasma  eonsidered  in 
this  work,  the  theory  will  eorrespondingly  evolve  and  beeome  more  eomplex. 

This  ehapter  will  traek  the  evolution  of  helieon  waves  from  simplified  uniformity 
assumptions  where  solutions  for  the  wave  field’s  by,  ba  bz  are  solved  analytieally.  When 
eollisional  effeets  and  finite  eleetron  mass  are  eonsidered,  the  internal  plasma  fields 
reveal  a  seeondary  wave  braneh  (the  Trivelpieee-Gould  mode)  whieh  ean  also  be  studied 
analytieally.  Finally,  when  a  non-uniform  radial  density  distribution  is  eonsidered,  the 
differential  equation  to  be  solved  beeomes  non-linear  and  eomplex  enough  that  only 
numerieal  solutions  are  valid.  However,  these  will  later  be  shown  to  provide  a  superbly 
aeeurate  deseription  of  the  helieon  wave  fields  under  investigation  in  this  work.  While 
some  groups  have  eonsidered  more  elaborate  theories  sueh  as  ineorporating  ion  mass 
effeets  [1-3]  or  wave  propagation  based  on  radially  loealized  modes  [4,  5],  the 
experimental  results  aeeompanying  these  deseriptions  is  ambiguous  or  ineomplete  as  far 
as  deseribing  the  wave  patterns.  This  is  not  to  say  that  they  may  not  also  aeeurately 
deseribe  helieon  wave  propagation  when  studied  more  thoroughly.  However,  the 
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mathematical  complexities  introduced  do  not  provide  for  any  more  ease  to  the  physical 
understanding  or  significantly  alter  the  results  obtained  than  those  when  considering  a 
non-uniform  radial  density  distribution. 

2.1  Plasma  Waves 

Due  to  the  multitude  of  parameters  that  can  change  the  operational  modes  within  plasma, 
the  subject  of  plasma  waves  is  often  termed  a  ‘plasma  zoo’.  Considering  a  two  species 
plasma,  ions  and  electrons,  immersed  in  a  uniform  magnetic  field  and  homogeneous  in 
space.  Maxwell’s  equations  with  the  linearized  equation  of  motion  yield  the  cold  plasma 
dispersion  relation.  The  mathematical  derivation  of  the  equations  for  the  cutoffs  and 
resonances  of  various  plasma  waves  is  given  in  Appendix  A.  The  graphical 
representation  was  first  given  by  Clemmow  and  Mullaly  in  1955  [6]  and  later  updated  by 
Allis  in  1959  [7];  commonly  referred  to  as  a  CMA  diagram  in  Figure  2.1.  To  accurately 
show  the  entire  frequency  spectrum,  an  ion-electron  mass  ratio  of  5  was  selected. 

The  axes  on  the  CMA  diagram  in  terms  of  the  plasma  frequency  and  electron  cyclotron 

frequency  can  be  transformed  according  to  their  definitions  as  density  and  magnetic  field. 

2 

Here  the  electron  plasma  frequency  is  given  by  co  = - ^  and  the  electron  cyclotron 

cB 

frequency  is  given  by  =  — -  .  The  n  vs.  Bq  diagram  is  given  in  Figure  2.2  for  Argon 

gas  and  a  driving  frequency  of f=  13.56  MHz.  Within  each  region,  only  certain  waves 
are  allowed  to  propagate  within  the  dielectric.  A  detailed  description  of  the  various 
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regions  which  allow  right-  or  left-  handed  waves  and  either  ordinary  or  extraordinary 
waves  to  propagate  is  given  by  Stix  [8]  and  Swanson  [9]  among  others. 


P  =  0  ^  Increasing  Frequency 


We  will  elaborate  on  one  result  of  relevance  to  the  current  work  and  that  is  the  dispersion 
relation  obtained  for  right-handed  wave  propagation  at  an  angle  0  with  respect  to  a 
magnetic  field  in  a  +z  direction.  The  result  drawn  on  from  Appendix  A  is  given  by 


n^=R  = 


=  \- 


CO [co-co^^  cos  0) 


(2.1) 


where  ion  effects  have  been  neglected  and  the  index  of  refraction  ‘n’  is  defined  by 

kc 


n  =  - 


CO 


(2.2) 


and  k  =  kj-^j.  is  the  total  wave  number  according  to 
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k. 


Simplifying  where  k^=kcos6  and  where  yields 


,  0)pe 

Ktot 


k.  coc 


(2.3) 


This  is  the  basic  dispersion  relation  for  right-handed  whistler  wave  propagating  in  free 
space.  Later,  the  dispersion  relationship  between  whistler  waves  described  by  this  result 
and  those  for  helicon  wave  propagation  will  be  discussed. 


«  [m 

Figure  2.2.  CMA  diagram  in  terms  of  n  vs.  B. 

Thus  far,  we  have  described  plasma  waves  that  propagate  in  an  unbounded  medium  and 
have  only  put  a  restraint  on  the  operating  gas  and  driving  frequency.  In  the  following 
sections,  the  subject  of  bounded  laboratory  helicon  waves  will  be  addressed.  We  will 
describe  waves  that  propagate  between  the  lower  hybrid  frequency  and  electron  cyclotron 
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2 

frequency <a)<a)^^.  This  frequency  space  within  the  -  ffifce  and  n  -  B  space  of 

the  CMA  diagram  envelopes  the  conditions  where  helicon  waves  propagate.  Although 
the  CMA  diagram  is  a  cursory  analysis  for  studying  bounded  plasma,  it  provides  one  of 
the  simplest  representations  toward  understanding  the  restrictions  when  designing  a 
laboratory  system  to  drive  helicon  plasma. 

2.2  Bounded,  Collisionless,  Uniform  Density  Dispersion  Relation 

Consider  a  cylindrical  plasma  discharge  of  radius  ‘a’,  uniform  radial  density 
distribution Hq  (r)  =  Hq  ,  subject  to  an  infinite  homogenous  magnetic  field  Bq  in  the  +z 
direction.  Maxwell’s  equations  are  then 


VxE  = - 

dt 

(2.4) 

Vxb  =  pJ 

(2.5) 

< 

II 

o 

(2.6) 

Here,  we  have  neglected  the  displacement  current  on  the  assumption  that  the  wave  phase 
velocity  «  the  speed  of  light,  i.e., 

\=j«c  (2.7) 

Assuming  for  now  that  the  entire  plasma  current  is  carried  by  the  E  x  B  motion  of  the 
electrons  where  the  driving  frequency  is  much  smaller  than  the  electron  cyclotron 
frequency  co  <K  and  is  much  larger  than  the  lower  hybrid  frequency 

CO  »  co^^  =  ^Jco^e^^i  such  that  ion  motions  can  also  be  neglected,  the  Ohm’s  law 
expression  to  be  evaluated  is  given  by 
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enoE  =  jxB^, 


(2.8) 


Additionally,  we  have  assumed  here  that  the  resistivity  is  zero,  i.e,  r]  =  0,  where  r]  is 
defined  by 


'  2 

ne 

Essentially  eollisional  effeets  have  been  negleeted  where  the  eleetron  mass 
effeets  will  be  ineluded  in  the  subsequent  seetion.  Returning  to  Maxwell’s 
(2.4)  -  (2.6)  and  Ohm’s  law  (2.8),  E*  order  wave  perturbations  of  the  form 

exp  i{m9  +  k^z  -  cot) 


(2.9) 

zero.  These 
equations 


where  ::$>  b  simplify  the  expressions  for  the  wave  magnetie  fields  sueh  that 


b  = 


f  A"' 

CO 

V  Bq  j 


Vxb 


(2.10) 


or 

ab  =  Vxb  (2.11) 

where  the  total  wave  number  is  defined  as 


CO  jUo^% 
k,  Br, 


(2.12) 


Inserting  the  definitions  of  the  eleetron  plasma  and  eyelotron  frequeney  yields 


COC 


CO  //penp  _  CO  (Ope 

k,  B„  k,  coc^ 


(2.13) 


whieh  is  the  identieal  result  to  (2.3)  for  right-handed  whistler  wave  propagation  in  free 
spaee.  Taking  the  eurl  of  (2.10)  yields  a  seeond  order  differential  equation  for  the  plasma 
wave  fields 
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(2.14) 


where 

T^=kl=  -kl=a^-kl  (2.1 7) 

The  differential  equation  (2.16)  is  Bessel’s  equation  subjeet  to  the  finite  boundary 
condition  at  the  origin  so  that 

K{'-)  =  CJ,(Tr)  (2.18) 

The  r  and  ^components  of  the  wave  magnetic  field  are  obtained  from  (2.10)  as 

j  /  \  iC  TVICX  j  /  \  j  j  '  /  \  /r\  -t  r\\ 

^r{r)  =  ^  - K{r)-kA{r)  (2.19) 

I  y  r  ; 

=  (2.20) 

I  y  r 

The  second  boundary  condition  is  given  by 

b^(a)  =  0  (2.21) 

regardless  of  whether  the  plasma  dielectric  is  surrounded  by  a  conducting  or  insulating 
boundary  [10].  This  is  demonstrated  by  the  fact  that  substitution  of  (2.10)  into  (2.5) 
yields 

MoJ  =  (2.22) 
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so  that  on  an  insulating  boundary  (a)  =  0  ^  (a)  =  0  .  On  a  perfect  conducting 

boundary  (a)  =  0  ,  and  (2.4)  gives 

KE,[r)  =  -a)b^[r)  (2.23) 

so  that  £'^(a)  =  0^b^(a)  =  0.  Therefore,  the  second  boundary  condition  is  obtained 
from  (2.19) 

^K{a)  =  Kb'^{a)  (2.24) 

which  assures  a  unique  value  to  kz  provided  other  inputs  for  tube  radius  ‘a’,  plasma 
density  no,  static  magnetic  field  Bq,  antenna  configuration  mode  ‘m’,  and  driving 
frequency  For  example,  consider  a  2.5  cm  radius  tube,  with  a  uniform  900  Gauss 
magnetic  field,  and  a  m  =  +1  helical  antenna  driven  at  13.56  MHz.  The  boundary 
condition  given  by  (2.24)  provides  akz  =  11.58  Az  =  54.25  cm  via 

4=2^^  (2.25) 

defining  the  axial  helicon  wavelength.  Finally,  the  internal  plasma  wave  fields  can  be 
solved  from  (2.18)  thru  (2.20)  and  are  shown  in  Figure  2.3. 
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Figure  2.3.  Internal  b„  bg,  and  b^  fields  for  a  =  2.5  cm,  m  =  +1  antenna  driven  at  13.56  MHz  in  a  900 
Gauss  uniform  magnetic  field  with  a  uniform  radial  density  of  no  =  1  x  10^*  m'^. 


2.3  Uniform  Density  with  Collisions:  The  Trivelpiece-Gould  Mode  Revealed 
Previously,  we  assumed  the  component  Ez  was  zero  and  helicon  wave  propagation  in  the 
+z-direction  was  undamped;  i.e.,  the  plasma  current  j  was  driven  hy  E  x  B  motion  of 
electrons.  However,  when  considering  wave  damping  due  to  collisions  or  Landau 
damping,  an  Ez  field  arises  to  move  electrons  in  the  +z  direction.  This  is  taken  into 
consideration  in  the  jz  component  of  the  plasma  current  for  motion  parallel  to  the  Bq  field. 
When  this  is  considered,  the  Ohm’s  law  equation  (2.8)  has  to  be  altered  to  include  this 
component.  Specifically, 

en^E  =  jxB^ - ico  +  iv\j  (2.26) 

e 


or 


„  jxB^  im  ,  s- 
E  =  - ^ ^ — [(o  +  iv)  j 


en^ 


2 

e  Wn 


(2.27) 
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Maxwell’s  equations  (2.4)  thru  (2.6)  are  then  used  with  Ohm’s  law  (2.27)  and  a  similar 
formulation  to  the  uniform  density  ease  previously  considered  is  followed.  Taking  the 
curl  of  (2.27)  and  substituting  into  (2.4) 


VxE 


en^ 


-^^((»  +  zV)Vx  j  =  icob 
e 


(2.28) 


en^ 


+  zV)  V  X  j  =  icob 

e  zzo 


(2.29) 


taking  the  curl  of  (2.5) 


V  X  V  X  6  =  //qV  X  j 


(2.30) 


substituting  (2.5)  and  (2.30)  into  (2.29)  for  the  plasma  current  j  yields 


mia  +  iv)  -  k  B„  r  r 

-A - ^-VxVxb - ^^Vxh  +  h  =  0 


e  HqJUqCO 


eriQjUQCO 


(2.31) 


utilizing  the  previous  definition  of  a  from  (2.12)  with  the  electron  cyclotron  frequency 

COce 

(<jo  +  iv\  -  1  -  - 

- ’-VxVxb - Vxh  +  h  =  0  (2.32) 

co^gCtk^  a 


additionally  defining 


7  = 


(nz  +  zV) 
cok. 


(2.33) 


so  that  (2.32)  becomes 


xV  xb-V  xb  +  ab  =  0 
and  the  characteristic  equation  to  solve  is 


+  «  =  0 


(2.34) 


(2.35) 


where  roots  are  given  by 
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and  the  square  root  term  has  been  Taylor  series  expanded  to  seeond  order.  Taking  the 
‘bottom’  root  of  (2.36)  gives 


y02  =«(!  +  «?')  (2.37) 

and  when  the  eleetron  mass  is  zero,  i.e,  mg  =  0;  y  0  and  (2.37)  gives  P2  =  a,  whieh  is 
the  previously  derived  helicon  root;  (2.12).  When  finite  electron  mass  effects  are 
considered  P2  becomes 


A  =  « 


1  +  a 


~h  I 


77 


(2.38) 


where  the  resistivity  is  defined  as 


7  = 


mv 


n^e 


(2.39) 


and  the  collision  frequency  here  is  assumed  to  the  electron-ion  collision  frequency  Vg, 


=^2.906x10 ‘Mn A  nX 


(2.40) 


Returning  to  the  characteristic  equation  (2.35),  the  second  root  is  revealed  when  we 
expand  a  and  y  using  (2.13)  for  a  and  (2.33)  for  y  and  taking  the  real  part  of  y 


Q)  T 

13^-13  +  — r^  =  0 


cok. 


co.„kc 


(2.41) 


k.  = 


CO/  f 


p 


V 


^  y 


(2.42) 


Here,  /3  =  kj-^j  analogous  to  where  a  =  kjQj  for  the  uniform  density  undamped  case. 


The  ‘new’  total  wavenumber  P  now  includes  the  damping  term.  The  effects  are  best 


illustrated  as  demonstrated  by  the  dependence  of  kz  on  P  at  different  magnetic  field 
strengths.  For  the  first  example  a  density  of  Hq  =  1  x  is  used.  The  graphical  results 


are  shown  in  Figure  2.4. 
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Figure  2.4.  Separation  of  TG  branch  and  Helicon  branch  for  no  =  1  x  10  m  . 


Analysis  demonstrates  that  as  the  magnetic  field  is  increased  a  second  branch  to  the 
dispersion  relation  is  observed.  To  the  left  of  the  minimum  kz  value  is  the  helicon  branch 
and  to  the  right  is  the  Trivelpiece-Gould  branch.  The  difference  between  the  two 
branches  can  be  seen  in  Figure  2.5. 
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Figure  2.5.  Difference  when  finite  electron  mass  effects  are  considered.  The  secondary  TG  branch 
is  revealed.  This  is  for  Bo  =  900  Gauss  and  no  =  1  x  10^*  m'^. 


The  diagonal  line  where  k^=  (3  occurs  from  the  axial  wavenumber  being  required  to  be 
less  than  the  total  wavenumber,  i.e,  +kl=  k^^j.  =  /3^  as  k^  ^  0  and  determines  a 
maximum  kz  for  helicon  waves.  Additionally,  the  minimum  value  of  kz  is  obtained  by 
differentiating  (2.42)  for  kz,  equating  the  result  to  zero  and  obtaining  the  value  of  P 
according  to 


dk^ 


CO 


1  o 


dp  0)^^  P  co^^  c 


=  0 


CO 


pe 


=min 


so  that  a  minimum  in  kz  occurs  at 


,  .  .  I - 

=2— = 


(2.43) 

(2.44) 


(2.45) 
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Therefore,  knowing  the  minimum  and  maximum  values  of  the  axial  wavenumber  for 
helicon  waves  at  a  designated  density  and  magnetic  field  places  an  additional  restraint  on 
the  axial  wavelength  range  according  to 


(2.46) 


Different  densities  for  Bo  =  900  Gauss  yield  the  following  dispersion  curves,  Figure  2.6. 


Figure  2.6.  TG-H  dispersion  curves  different  densities  at  Bo  =  900  Gauss. 


Here  the  minimum  kz  value  is  decreasing  as  the  density  is  decreased  according  to  (2.45) 
and  the  axial  wavelength  range  is  increasing  according  to  (2.46). 


The  TG-branch  can  be  better  understood  by  considering  that  as  the  magnetic  field  is 
decreased,  the  two  branches  (TG-  and  H-)  begin  to  come  together,  until  the  magnetic 
field  is  so  low  that  only  the  TG-branch  can  propagate.  Referring  back  to  the  CMA 
diagram  for  no  vs.  Bo,  Figure  2.2  will  demonstrate  that  the  TG-branch  is  an  electron 
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cyclotron  wave.  This  can  be  viewed  for  n,,  =  1  x  1 0'“*  m  ^  as  beginning  on  the  no  vs.  Bo 
diagram  at  ~  1  kGauss  and  ‘traveling’  vertieally  downward  (deereasing  magnetie  field) 
toward  the  R  oo  resonanee,  which  is  an  electron  cyclotron  wave  at  5  Gauss  for 
Argon  at  this  density.  However,  for  stronger  magnetic  fields  (on  the  order  of 1  kGauss) 
Figure  2.4  shows  that  the  two  branehes  are  well  separated.  When  this  is  the  ease,  for 
each  value  of  kz  there  exist  two  values  for  p.  Regardless,  for  the  TG-braneh 
where  ,  the  total  wavenumber  P  is  dominated  by  the  perpendieular  eomponent 

and  henee  the  TG  wave  is  a  radial  eleetron  eyelotron  wave  of  short  wavelength  (large 
k^  )  given  by 


A,  =2^/ 


/ 

(^±  ~  PjQ  ) 


(2.47) 


As  an  example,  assuming  the  above  eonditions  for  density  =lxl0'^m  ^and  magnetie 
field  Bo  =  900  Gauss),  P  varies  from  43.78  m"'  to  8090.16  m"'  with  the  helieon  braneh 
inereasing  until  k^  at  P  =  595.07  m'\  Therefore,  for  a  plasma  with  an  axial 

wavenumber  of  kz  =  25  m'^  the  two  values  of  P  are  77.53  m'^  and  4567.25  m'\ 

Aeeording  to  (2.25)  the  axial  helieon  wavelength  will  be  kz  =  25. 13  em.  For  Ptg  = 
4567.25  m'^  the  total  wavenumber  again  is  dominated  by  the  perpendieular  eomponent, 
i.e.  P  =  kj.Qj.  =k^  +k^Kk^  and  the  TG  wave  has  a  perpendieular  wavelength  of  = 
0.14  em.  Conversely,  when  the  branehes  are  not  well  separated,  i.e.,  at  a  low  magnetie 
field.  Bo  =  100  Gauss,  P  varies  between  134.26  m'^  and  2637.31  m'^  with  the  helieon 
braneh  again  inereasing  until  k^  at  P  =  595.07  m'^  If  kz  =  60  m'^,  then  (2.25)  gives  kz  = 

10.47  em  while  the  eorresponding  Ptg  =  790.84  m'^  and  X^  =  0.79  em. 
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The  wave  fields:  br,  be,  and  bz,  with  finite  eleetron  mass  now  ineluded,  will  not  greatly 
alter  the  wave  patterns  [11].  This  is  a  eonsequenee  of  the  additional  terms  only  altering 
the  total  wave  number  kjoT^y  a  numerical  factor  as  shown  in  (2.37)  and  (2.38).  This  can 
be  recognized  that  in  the  limit  of  strong  magnetic  fields  ('-'  1  kGauss)  the  dispersion 
relation  returns  only  the  helicon  wave  portion,  namely  (3  =  a  .  It  is  not  until  we  consider 
any  perpendicular  or  non-uniform  effects  that  the  plasma  wave  patterns  will  substantially 
change. 


2.4  Non-uniform  Radial  Density  Distribution:  Numerical  Solutions 

Up  to  this  point  a  uniform  radial  density  distribution  has  been  assumed.  In  what  follows 
we  will  again  neglect  the  wave  damping  aspect  in  order  to  simplify  the  analysis,  however 
a  non-uniform  radial  density  will  be  used.  Maxwell’s  equations  (2.4)  thru  (2.6)  remain 
unaltered  except  that  the  displacement  current  will  be  included  as 


Vx6 


BE 


juJ-icD^^s^E 


(2.48) 


Additionally,  Ohm’s  law  must  now  change  to  include  the  radial  distribution  according  to 


en 


{r)E  =  jxB^ 


(2.49) 


or 

E  =  (2.50) 

en[r ) 

The  analysis  follows  the  exact  formulation  as  the  uniform  density  distribution  except  that 
some  of  the  curl  and  vector  dot  products  must  take  into  account  the  radial  density 
component.  Therefore,  combining  Maxwell’s  equations  with  Ohm’s  law  yields 
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where 


ab  = 


+  1 

V 


z  +  ^{klb><z-a'{r^b^z^ 


(2.51) 


and 


CO  jUQen[r) 


(2.52) 


(2.53) 


Looking  at  (2.51)  the  displacement  current  gives  rise  to  the  terms  including  ko  where 


,  CO/ 

_ 0_ _  /  C  _ phase. 

k  ^  k  ~  /' 


(2.54) 


and  is  typically  on  the  order  of  0.01  for  laboratory  helicon  plasma  and  can  be  neglected 


tx  / 

by  letting  %  ^  0  since  v  <K  c ;  similar  to  the  assumption  in  (2.7).  Similarly,  the 


density  gradient  gives  rise  to  the  last  term  in  (2.5 1)  where  a  '(r)  oc  n  '(r)  such  that  the 

derivative  would  be  zero  for  constant  radial  density.  Returning  to  (2.51),  neglecting  the 
displacement  current  and  assuming  a  constant  density  will  return  equation  (2.10)  as 


ab  =  {v xb^  +^Vxhj  z  =  V xb  (2.55) 

To  obtain  the  individual  components  of  b,  (2.51)  is  expanded  as 


a/  =—b^-ikjbg 
r 

(2.56) 

(^bg=-l)^+ikjb^ 

(2.57) 

(^K=y{^^{rbg)  -—b^-ya'b^ 
yr  r  )  k^ 

(2.58) 

and  (2.56)  -  (2.57)  can  be  reduced  in  terms  of  b^ 
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(2.59) 


/3b,=—ab,+ikjt;, 


I3b,=-ah,--kjb 


Z  Z!  z 

r 


(2.60) 


where 


P  =  oc^-ky 


(2.61) 


r  =  \-  ^ 

k. 


(2.62) 


To  obtain  the  differential  equation  for  the  plasma  fields,  (2.59)  and  (2.60)  are  substituted 
into  (2.58)  and  reduced  to 


+  f{r)biy)  +  g{r)by)  =  Q 


(2.63) 


where 


N  1  2aa' 

J  \r)  = - 

^  ’  r  p 


(2.64) 


/X  P  m  a'  (  ^kly^ 

S{r)  =  ^ - Y-- - 1  +  ^- 

Y  r  k^yA  p 


(2.65) 


The  same  boundary  conditions  as  before  are  utilized,  namely  the  wave  field  is  finite  at  the 
origin  and  the  radial  component  6^  (a)  =  0  vanishes  on  the  boundary;  (2.21)  which 
follows  from  (2.51)  for  the  vanishing  radial  plasma  current  jr 


ab  =  {Vxb^^= 


(2.66) 


Therefore,  (2.59)  can  be  evaluated  for  6^  (r  =  a)  =  0 


fYl  ' 

—  a{r  =  a)b^  (r  =  a)  +  k^yb^  (r  =  a)  =  0 


(2.67) 
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The  radial  density  profiles  studied  for  this  analysis  were  given  by  Gaussian  profiles 
according  to 

=  (2.68) 
where  ‘w’  is  the  waist  size  of  the  profile  or  V2  the  width  of  the  radial  distance  for  the 
density  to  drop  to  1/e  of  no.  Taking  =  lx  the  radial  density  profiles  are  shown 

in  Figure  2.7. 


Radius  [cm] 

Figure  2.7.  Gaussian  profiles  of  waist  size  V  for  Uo  =  1  x  10^’  m'^  and  tube  radius  a  =  2.5  cm 

The  following  conditions  were  used  as  a  baseline  in  the  numerical  simulation  of  the  bz 
profiles  according  to  (2.63):  Bo  =  900  Gauss,  m  =  +1  antenna, /=  13.56  MHz,  and  a  =  2.5 
cm.  The  differential  equation  was  solved  using  Mathematica:  NDSolve  function  (see 
Appendix  B  for  the  code).  It  uses  an  interpolation  scheme  with  a  specified  step  size.  In 
this  case,  the  step  size  was  0.005  cm  which  provided  approximately  500  points 
between  r(0,a^  .  The  bz  results  for  each  of  the  radial  density  profiles  are  given  in  Figure 
2.8. 
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Radius  \cm\ 

Figure  2.8.  bz  Helds  for  Gaussian  radial  density  profiles. 

Bfl  =  900  Gauss,  m  =  +1  antenna,/  =  13.56  MHz,  and  a  =  2.5  cm 

The  derivative  of  bz  was  numerically  solved  in  the  same  manner  such  that  the  remaining 

br  and  be  components  could  be  solved  according  to  (2.56)  and  (2.57).  These  results  are 

shown  in  Figure  2.9  and  Figure  2.10. 


Radius  [cm] 

Figure  2.9.  -ibr  fields  for  Gaussian  radial  density  profiles. 

Bo  =  900  Gauss,  m  =  +1  antenna,/  =  13.56  MHz,  and  a  =  2.5  cm 
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Figure  2.10.  be  fields  for  Gaussian  radial  density  profiles. 

Bfl  =  900  Gauss,  m  =  +1  antenna,/  =  13.56  MHz,  and  a  =  2.5  cm 


The  accuracy  of  the  results  for  the  numerical  simulation  in  Figure  2.8  through  Figure  2.10 
is  verified  when  the  uniform  density  fields  are  compared  against  those  obtained  in  Figure 
2.3  and  complete  agreement  is  achieved.  Additionally,  both  boundary  conditions  are 
satisfied  when  bz  is  solved  for,  namely  the  radial  component  of  the  b,-  fields  go  to  zero  for 
each  profde  studied;  Figure  2.9. 

Returning  to  the  boundary  condition  (2.67)  and  expanding  a,  y  and  rearranging 


(2.69) 


shows  that  when  the  boundary  condition  is  applied  for  a  given  density  profile,  the  only 
free  parameter  is  the  axial  wavenumber  kz.  In  particular,  kz  is  adjusted  until  the  condition 
is  satisfied  for  the  plasma  wave  fields.  The  impact  of  the  radial  density  profiles  can  be 
seen  in  Figure  2.1 1. 
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Figure  2.11.  Influence  of  density  profiles  on  wavenumber  between  two  wave  branches, 
no  =  1  X  10*’  m  *,  Bo  =  900  Gauss,  m  =  +1  antenna,/  =  13.56  MHz,  and  a  =  2.5  cm 


Essentially,  the  values  act  as  eigenmodes  between  the  branches.  The  P  value  for  the 
TG  branch  then  becomes  k\=  0^  -  0  =  0  .  This  again  shows  the  TG-branch  to  be  a 
dominantly  radial  wave.  For  these  values  of  P  the  radial  wavelength  is  extremely  short; 
varying  between  (0. 15cm  -0.26cm)  whereas  the  axial  helicon 

wavelength  ^(26.79 cm -46.23 cm) . 


2.5  Chapter  Summary 

This  chapter  has  summarized  the  derivation  of  the  helicon  wave  radial  b-profiles.  The 
two  primary  considerations  were:  finite  vs.  zero  electron  mass,  and  uniform  vs.  non- 
uniform  radial  density  distributions.  These  derivations  have  been  independently 
considered  in  the  variously  referenced  journal  articles;  however  this  chapter  compiles  the 
results  into  the  comprehensive  layout  presented  here. 
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Inclusion  of  electron  mass  revealed  an  additional  wave  of  very  high  perpendicular  wave 
number;  the  TG-wave.  The  associated  wavelength  with  this  radial  TG-wave  is  typically 
less  than  1  cm.  When  considering  a  non-uniform  radial  density,  the  helieon  wave  b- 
profiles  notably  alter  their  radial  structure.  As  the  density  profile  more  centrally  peaked, 
so  too  do  the  radial  b-profiles. 

This  chapter  has  much  more  thoroughly  quantified  the  results  for  the  axial  wavelength 
than  has  been  previously  considered.  We  have  shown  the  helicon  wavelength  varies  over 
the  lO’s  cm  range  dependent  on  the  density  and  applied  magnetie  field.  Additionally,  we 
demonstrated  that  for  each  Gaussian  density  profile  considered,  the  axial  wavelength 
must  be  varied  until  the  boundary  condition  is  satisfied.  Finally,  the  signifieance  of  this 
chapter  is  manifest  in  the  derived  measureable  quantities:  br,  ba  bz,  and  n(r)  whieh 
will  be  verified  in  the  laboratory. 
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CHAPTER  3:  EXPERIMENT  AND  TEST  SETUP 


To  maintain  consistency,  the  experiments  performed  in  this  work  were  designed  and 
scaled  based  off  of  earlier  helicon  efforts  [1-4],  This  also  allowed  validation  of 
preliminary  work  and  later  allowed  for  comparison  with  current  results.  The 
experimental  work  performed  here  was  conducted  at  the  Air  Force  Research  Laboratory’s 
(AFRL)  Electric  Propulsion  (EP)  facility  in  Edwards  Air  Eorce  Base,  California.  This 
chapter  will  describe  the  apparatus  used  for  testing. 

3.1  Vacuum  System 

All  testing  was  conducted  with  a  quartz  tube  attached  to  a  large  diffusion  chamber  (the 
details  of  each  tube  will  be  discussed  with  the  antennae  tested).  The  diffusion  chamber  is 
a  cylindrical,  stainless  steel,  0.5  meter  diameter  by  1.0  meter  long  vacuum  chamber.  It 
has  24  access  ports  /  flanges  of  varying  diameter  for  diagnostic  feed-throughs,  pumping 
connections,  and  internal  viewing  and  access.  The  schematic  is  shown  in  Eigure  3.1. 


Foreline  Convectron  Turbo  Pump 


Roughing  Pump 

Figure  3.1.  Schematic  of  vacuum  diffu$iou  chamber  with  pump  locatious 
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The  chamber  is  pumped  by  a  Varian  TV-300  HT  turbomolecular  vacuum  pump  with  a 
pumping  speed  of  250  1/s  backed  by  a  dry  scroll  pump.  The  chamber  and  foreline 
pressure  are  monitored  with  two  Varian  convectron  gauges  from  atmosphere  to  1  x  10' 
Torr.  Below  this  pressure  a  Varian  cold  cathode  gauge  is  used  to  monitor  the  chamber  to 
a  base  pressure  of 2.4  x  10'^  Torr.  For  tests  in  this  work,  the  pressure  range  varied  from 
0.5  -  10  X  10'  Torr.  Neither  the  convectron  nor  cold  cathode  gauge  provided  reliable  or 
accurate  pressure  measurement  in  this  range.  Therefore,  a  MKS  629D  Baratron 
capacitance  manometer  with  a  full  scale  range  of  10'"^  -  10'^  Torr  was  used  for  the 
intermediary  pressure  monitoring.  This  gauge  was  positioned  on  the  ‘downstream’  end  of 
the  antennae  /  source  at  a  location  in  between  the  quartz  tube  and  diffusion  chamber.  All 
vacuum  equipment  is  mounted  on  a  120  x  300  x  30  cm  optics  table  supported  by 
vibration  suppressing  pneumatic  legs.  Additionally,  the  turbomolecular  pump  is 
vibrationally  isolated  with  an  elastomer  flange  inserted  between  the  foreline  and  diffusion 
chamber.  The  optics  table  not  only  provided  a  ‘good’  environment  for  suppressing 
mechanical  noise  but  also  provided  a  solid  electrical  ground  for  the  vacuum  chamber  and 
other  equipment. 

3.2  Experiment 

Power  to  the  antennae  was  supplied  by  an  ENI  AlOOO  broadband  power  amplifier.  The 
unit  is  an  all  solid  state  linear  Class  A  amplifier  which  takes  a  0  -  1  Volt  input  signal  and 
outputs  up  to  1000  Watts  over  0.3  -  35  MHz  and  50  fl  input  /  output  impedance.  The 
amplifier  input  is  taken  from  a  Tektronix  AFG  3102  100  MHz  arbitrary  waveform 
function  generator.  The  forward  and  reflected  power  to  the  amplifier  was  monitored  with 
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a  Bird  4021  directional  power  sensor  (1.8  -  32  MHz;  0.3  -  1000  Watts)  connected  to  a 
Bird  4421  RF  power  meter.  The  matching  network  was  a  Manitou  Systems  ACM- 
1500L-12-MC-VC.  It  is  an  L-network  capable  of  delivering  up  to  1500  Watts  at  13.56 
MHz  with  automated  feedback  to  tune  two  variable  capacitors  to  ensure  the  total  load  to 
the  amplifier  is  matched  to  50  fl.  The  capacitor  in  series  with  the  plasma  is  a  COMET 
50  -  500  pF  8  kV  and  the  shunt  capacitor  to  ground  is  a  COMET  20  -  2000  pF  5  kV 
capacitor.  There  is  also  an  inductor  in  series  with  the  plasma  load  that  is  variable  to  14 
|uH  but  during  operation  remains  at  a  set  value.  It  is  only  variable  in  the  sense  that  the 
internal  connectors  to  the  inductor  can  be  moved  between  7  consecutive  taps  which 
increase  or  decrease  the  total  number  of  turns  which  RF  current  is  driven  through.  The 
inductance  varies  ~  2  pH  /  tap,  however,  for  all  measurements  and  testing  conducted  in 
this  work,  the  lowest  value  of  inductance  was  used  at  2  pH.  A  schematic  of  the  RF 
power  system  is  shown  in  Figure  3.2. 


Figure  3.2.  RF  power  delivery  system  incorporating  the  auto-match  network  circuit. 

All  tests  were  conducted  using  UHP  Grade  (99.999  %)  Argon  gas.  The  mass  flow  rate 
was  controlled  through  an  Alicat  Scientific  MC-IOOSCCM-D  (Ar)  flow  controller.  Flow 
rates  were  varied  for  operating  test  pressures  from  0.5  -  10  mTorr  and  were  performed  for 
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100  -  500  Watts  input  power  at  300,  600,  and  900  Gauss  static  Bo-field.  A  schematic  of 
the  test  setup  with  electrostatic  shield  and  diagnostic  entry  ports  is  shown  in  Figure  3.3. 

90  GHz  Microwave  0-1200  Gauss 


0-100  MHz 
Waveform  Generator 


Figure  3.3.  Schematic  of  AFRL  Experimental  Helicon  research  station. 

The  static  magnetic  field  is  generated  by  4  glycol  cooled  magnets  that  can  supply  a 
continuous  uniform  field  up  to  1200  Gauss  (custom  designed  by  Arnold  Magnetics). 

They  operate  using  a  100  Volt  100  Amp  EMS  DC  power  supply.  The  electromagnets 
span  a  distance  of  45  cm  and  provided  uniform  magnetic  field  over  approximately  30  cm. 
The  DC  magnetic  field  was  measured  with  a  F.W.  Bell  Gauss  meter  (Model  7030).  The 
determination  of  the  static  field  as  measured  on  the  quartz  tube  centerline  for  varying 
input  currents  is  shown  in  Figure  3.4. 
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Figure  3.4.  Axial  applied  magnetic  field. 


The  majority  of  experimental  measrrrement  was  eonducted  for  a  driving  frequency /= 
13.56  MHz  at  500  Watts  input  power  and  900  Gauss  axially  applied  magnetic  field.  The 
gas  used  was  Argon  at  a  flow  rate  of  47  seem  which  corresponded  to  an  operating 
pressure  of  approximately  10  mTorr.  These  were  the  conditions  where  previous 
measurements  for  the  helicon  wave  profiles  matched  those  predicted  analytically  [2,  3,  5, 
6]  as  shown  in  Chapter  2.  Additionally,  these  test  conditions  provided  signals  of 
adequate  strength  for  internal  probe  measurement  so  that  diagnostic  intrusion  was 
limited;  this  will  be  discussed  in  further  detail  in  Chapter  4-5. 


3.S  Antennae 

Previous  helicon  work  has  typically  utilized  three  different  antenna  configurations  [7]  . 
These  being  the  Nagoya,  Boswell,  or  twisted  Nagoya  (also  known  as  a  full  turn  or  half 
turn  helical  antenna).  The  three  antennae  are  shown  in  Figure  3.5  a,b,c. 
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a.  Nagoya  III  b.  1/2  turn  Twisted  Nagoya  c.  Boswell 


Figure  3.5  a,b,c.  Examples  of  different  antennae  used  in  helicon  research. 

Throughout  this  study,  the  V2  turn  twisted  Nagoya  or  V2  turn  helieal  antenna  was  tested. 
This  antenna  strueture  has  been  extensively  tested  and  provided  the  most  data  for 
eomparison  with  the  eurrent  work.  Additionally,  sinee  it  was  wrapped  as  a  right-handed 
antenna  with  respeet  to  the  statie  magnetie  field,  the  antenna  strueture  is  known  to 
propagate  the  m  =  +1  helieon  wave  fields  deseribed  throughout  Chapter  2.  If  the 
direetion  of  the  magnetie  field  is  reversed  then  the  antenna  would  be  left-handed  and  the 
m  =  -1  wave  fields  should  propagate. 

The  antennae  eonsisted  of  wrapping  eopper  strap  around  the  outside  diameter  of  quartz 
tubes  that  were  glass  welded  to  stainless  steel  vaeuum  feed-throughs.  This  meant  that  the 
antennae  eurrent  loop  was  in  atmosphere  whereas  the  indueed  eurrents  from  the  antenna 
were  within  the  eylindrieal  quartz  eavity  under  vaeuum.  The  hardware  eonsisted  of  three 
quartz  eylinders,  eaeh  45  em  length.  Two  are  straight  eylinders  with  an  outside 
diameter  of  1.5  inehes  (3.8  em)  and  2.5  inehes  (6.4  em),  while  the  last  is  a  eonieal  design 
starting  at  1.5  inehes  (3.8  em)  outside  diameter  and  ending  at  4.0  inehes  (10.2  em) 
diameter  as  shown  sehematioally  in  Figure  3.6.  This  quartz  tube  and  antenna  strueture 
ean  either  be  eonverging  or  diverging  with  respeet  to  the  magnetie  field  in  a  manner 
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similar  to  the  right-  and  left-  handed  helicity  of  the  helical  antennae.  Antenna  lengths 
tested  ranged  from  2.0  inches  (5.1  cm)  to  8.0  inches  (20.3  cm)  between  the  end  rings. 


Figure  3.6.  Conical  geometry  tube;  can  be  converging  with  Bo  or  diverging  with  Bq. 

The  primary  tube  tested  (the  6.35  cm  OD  cylindrical  tube  ~  5  cm  ID)  consisted  of  V2  turn 
helical  antenna  lengths  of  approximately  20.3  ±0.1  cm,  10.2  ±0.1  cm,  and  5.1  ±0.1  cm. 
The  antenna  lengths  and  their  location  within  the  static  magnetic  field  are  shown  in 
Figure  3.7.  Here,  the  antenna  length  measurement  is  reported  as  the  distance  between  the 
two  end-rings.  Typically,  this  has  been  how  the  anteima  length  has  been  defined; 
however,  we  will  refine  this  to  more  accurately  represent  the  true  antenna  length  in 
Chapter  5. 


Direction  of  Applied  Field  AAA 


Figure  3.7.  Applied  Magnetic  Field  with  Location  and  Lengths  of  Tested  Antennas. 
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CHAPTER  4:  DIAGNOSTICS 


As  stated  in  Chapter  2,  the  two  parameters  that  predominate  the  structure  of  helicon 
waves  are  the  plasma  bz  fields  and  the  radial  density  profiles  n(r).  It  will  be  these  two 
quantities  that  will  be  the  focus  of  measurement.  Additionally,  since  interest  lies  in  the 
three-dimensional  nature  of  the  radial  cross-section  and  axial  propagation,  internal  probes 
are  necessary.  Toward  this  end,  many  considerations  owing  to  the  intrusive  nature  of 
plasma  probes  will  be  taken  into  account  as  well  as  the  precarious  measurements  that  can 
persist  in  a  high  frequency  plasma  environment  generated  by  MHz  plasma.  This  chapter 
will  specifically  address  using  b-dot  probe  measurements  to  obtain  the  plasma  br,  be,  and 
bz  fields  as  well  as  discuss  the  use  of  a  RF  compensated  Langmuir  probe  and  its 
calibration  against  a  90  GHz  microwave  interferometer.  Finally,  an  error  analysis  will  be 
discussed. 

4.1  High  Frequency  b-dot  Probes^ 

Magnetic  induction  probes  have  been  used  as  a  diagnostic  to  measure  plasma  and  other 

transiently  induced  fields  since  the  1960’s  [2-4].  When  sized  appropriately  with  respect 

to  the  measuring  area,  a  spatial  distribution  of  internal  fields  can  be  measured  without 

much  disturbance  to  the  plasma  or  media  of  interest.  However,  when  designing  a  probe 

for  use  in  a  high  frequency  environment  the  probe’s  calibration,  sensitivity,  and 

frequency  response  must  all  be  considered.  This  is  especially  important  since  a  probe’s 

behavior  is  not  only  dependent  upon  its  physical  geometry  (coil  size,  etc)  but  also  by  non- 

*  The  work  on  High  Frequency  b-dot  probes  has  been  previously  published  by  Reilly  [1]  M.  P.  Reilly,  W. 
Lewis,  and  G.  H.  Miley,  “Magnetic  field  probes  for  use  in  radio  frequency  plasma,”  Review  of  Scientific 
Instruments,  vol.  80,  no.  5,  pp.  053508-5,  2009.  W.  Lewis  is  thanked  for  his  review  of  the  manuscript  and 
assist  is  data  collection. 
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ideal  factors  such  as  cable  lengths,  cable  types,  capacitive  pickup  due  to  plasma  or  other 
sources,  and  impedance  matching  of  the  probe  to  an  oscilloscope  or  other  measuring 
device.  More  thorough  qualitative  considerations  on  probe  characteristics  are  given  in 
References  [2-4],  In  this  paper,  each  of  these  will  be  addressed  and  quantified 
specifically  for  b-dot  probes  that  were  designed  to  take  measurements  on  13.56  MHz  RF 
cylindrical  plasma. 


4.1.1  Magnetic  Probe  Theory 

Inductive  probes  in  the  most  general  sense  consist  of  a  loop  of  wire  which  when  subject 
to  a  time  varying  magnetic  flux  obey  Faraday’s  Law 


F 


dt 


(4.1) 


where  F  is  the  induced  voltage  in  the  wire  and  Oh  is  the  magnetic  flux  through  the  loop 
or  coil.  We  consider  the  coil’s  cross  sectional  area  small  enough  that  the  flux  and 
magnetic  field  do  not  vary  over  this  area  and  obtain 


V=-A  ^ 
‘  dt 


(4.2) 


where  b  is  the  component  of  the  field  that  lies  along  the  axis  of  the  coil  and  ^e^is  the 
effective  area  over  which  the  magnetic  flux  is  measured.  Aeff  can  be  described  by 

4#  =  nA  (4.3) 


where  n  is  the  number  of  turns  in  the  sensing  coil  and  A  is  the  cross  sectional  area  of  each 
turn.  Substituting  for  yields, 

Tr  A  dB  , . 

F  =  -nA —  (4.4) 

dt 
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However,  this  will  most  likely  not  be  the  voltage  present  at  the  measuring  instrument 
(oscilloscope)  due  to  many  non-ideal  factors  that  cause  signal  attenuation  through 
impedance  mismatch.  A  thorough  analysis  and  theory  on  the  non-ideal  characteristics  of 
b-dot  probes  is  given  in  References  [2,  5]. 


To  obtain  a  large  voltage,  it  is  desirable  to  have  a  coil  with  a  large  cross-sectional  area 
and  many  turns.  However,  if  the  area  of  the  coil  is  large,  the  probe  can  disturb  the 
plasma  and  measurements  become  spatially  limited.  A  large  number  of  turns  in  the 
sensing  coil  with  a  large  radius  will  cause  the  inductance  to  rise.  This  is  seen  by  looking 
at  the  inductance  of  a  single  layer  solenoidal  coil  [3] 

^  f.  8w  w"'  5w®  35w* 

L  = -  1 - + - + - 

/  ^  3;r  2  4  16  64 

Clearly,  as  the  radius  of  the  coil  increases,  so  does  the  inductance.  Consequently,  when 
the  inductance  L  rises,  the  cutoff  frequency  fco  of  the  coil  decreases  according  to  [5] 


+  ... 


where  w  = 


(4.5) 


7 

f  —  0 

J  CO 


InL 


(4.6) 


Therefore,  a  compromise  must  be  found  that  provides  an  adequate  induced  voltage  (probe 
sensitivity;  large  n  &A)  while  maintaining  good  spatial  resolution  (small  A)  and  an 
adequate  frequency  response  (small  L  small  n  &  A). 


Considering  a  time  varying  field 6 (i)  =  hsin(<ui) ,  equation  (4.4)  for  the  resulting 

magnitude  of  the  magnetic  field,  as  measured  by  the  probe,  becomes 

V. 

(4.7) 

nAo) 
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Therefore,  in  order  to  eharaeterize  a  magnetie  induction  probe  for  use  at  a  specific 
frequency/  we  must  determine  the  product  nA.  This  is  also  known  as  the  calibration 
factor. 

4.1.2  Magnetic  Probe  Calibration 

Previous  work  has  shown  that  there  are  typically  three  ways  in  which  magnetic  probe 
calibration  is  accomplished  [2-4]: 

1)  direct  geometrical  inspection  of  the  coil  dimensions  (count  n,  measure  A) 

2)  measuring  the  output  voltage  when  the  coil  is  put  in  a  known  pulsed  field 

3)  comparison  of  the  probe  output  with  that  of  a  coil  of  known  dimensions,  i.e.,  a 
Helmholtz  coil,  and  calculate  what  should  be  generated  compared  with  what 
‘b’  is  measured. 

More  recently,  characterization  of  magnetic  induction  probes  have  also  utilized  network 
analyzers  in  conjunction  with  one  of  the  above  methods  to  gain  a  more  detailed 
measurement  of  the  probes  behavior  [5,6]  .  However,  there  exist  limitations  to  all  three 
methods  regardless  of  the  instrumentation  used. 

For  example,  with  direct  geometrical  inspection,  accurate  measurement  when  the  probe  is 
nominally  small  can  become  difficult.  Method  (2)  requires  having  a  known  field  near  the 
frequency  of  interest,  which  is  often  difficult.  Method  (3)  has  been  found  to  be  the 
preferred  method  of  calibration  amongst  many  authors  [2-4,  6];  however,  careful 
attention  must  be  used  when  this  method  is  considered. 
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A  Helmholtz  coil  is  utilized  most  frequently  due  to  the  uniform  field  at  the  center.  This 
field  is  given  by 


B.  = 


r  8  ^ 

V  «  J 

(4.8) 


where  ‘a  ’  is  the  coil  radius  and  is  equal  to  the  coil  separation  distance,  N  is  the  number  of 
turns  in  each  coil,  and  /  is  the  current  in  the  coils.  Immediately,  we  see  that  geometrical 
inspection  is  used  in  determining  the  number  of  turns  and  the  radius  of  the  calibrating 
Helmholtz  coil.  Additionally,  we  must  accurately  know  the  current  in  the  calibration 
circuitry  as  well  as  the  frequency  response  of  the  Helmholtz  coil  based  upon  its 
capacitance  and  inductance.  As  an  example,  utilizing  a  low  inductance  resistor  to 
measure  the  current  through  a  Helmholtz  coil,  we  find  deviations  for  the  ideal  impedance 
response  resulting  in  a  3  dB  point  of  66.5  MHz  for  a  120  Q  resistor  with  <0.1  qH  at  100 
kHz  as  shown  in  Figure  4.1.  (Noted  at  100  kHz  because  the  resistive  /  inductive  values 
will  vary  at  higher  frequencies  due  to  factors  such  as  skin  depth). 


Figure  4.1  Low  inductance  120  Q  resistor;  impedance  frequency  response. 
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Attention  must  also  be  paid  to  the  connection  wires  for  the  Helmholtz  coil  (as  well  as  in 
the  case  of  the  magnetic  induction  probe).  For  example,  the  current  source  for  the 
Helmholtz  coil  will  most  likely  have  some  connection  lead  wires  both  prior  to  the 
Helmholtz  windings  as  well  as  exiting  from  the  Helmholtz  windings.  Also,  if  the  current 
through  the  calibrating  coil  is  measured  as  the  current  through  a  low  inductance  resistor, 
the  connection  wires  that  sense  the  resistor  current  will  have  some  inherent  inductance  to 
them  and  likewise  a  cutoff  frequency  which  will  dictate  the  frequency  to  which  a 
Helmholtz  coil  can  be  utilized;  the  limitation  either  being  the  resistor  itself,  or  the  cutoff 
frequency  due  to  the  coils  inductance.  For  example,  coaxial  (BNC)  cable  is  typically 
rated  for  50  fl  use  to  a  few  GHz;  however,  if  the  load  is  not  50  fl,  or  breakouts  are  used, 
the  useful  frequency  range  drastically  changes.  If  alligator  clip  coax  breakouts  are  used 
to  either  measure  the  resistor  current  or  the  connectors  for  the  Helmholtz  source  current, 
then  the  useful  frequency  range  of  these  breakouts  is  limited  to  10  MHz,  as  shown  in 
Figure  4.2.  This  measurement  was  not  made  with  an  impedance  analyzer  to  illustrate  the 
frequency  limitation  of  connection  wires.  This  resistor  has  a  similar  inductance  to  that  of 
the  120  Q  resistor  previously  used;  however,  when  the  frequency  response  is  measured 
with  a  function  generator  and  oscilloscope,  we  find  that  the  frequency  is  limited  by  the 
inductance  in  the  connection  wires.  Therefore,  great  care  must  be  taken  when  utilizing 
method  (3)  as  a  method  of  calibration  because  errors  in  the  calibration  circuitry  and  setup 
can  easily  propagate  into  the  magnetic  probe’s  circuitry. 


62 


Figure  4.2.  Frequency  response  of  50  Q  resistor  and  coaxial  alligator  breakouts  as  measured 
with  a  function  generator  and  oscilloscope. 

Due  to  the  limitations  of  each  of  the  calibration  methods  described  above,  it  is  the 
author’s  preference  that  when  using  prefabricated  coils,  the  manufacturer's  geometrical 
specifications  be  used  for  the  number  of  turns  and  area  based  on  the  bobbin  size.  This  is 
typically  the  case  when  using  a  surface  mount  inductor  as  the  sensor  coil  as  is  used  in  this 
work. 


4.1.3  Magnetic  Probe  Capacitive  Pickup 

Due  to  the  nature  of  the  plasma  environment,  there  are  two  sources  from  which  a  probe 
will  measure  pickup;  inductive  signals  (which  are  due  to  the  magnetic  flux  lines  which 
pass  through  the  sensing  coil)  and  capacitive  pickup  (which  is  due  to  a  potential 
difference  between  the  plasma  discharge  and  coil).  These  differences  are  illustrated  in 
Figure  4.3. 
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Figure  4.3.  Illustration  of  inductive  and  capacitive  signals  present  when  magnetic  field 

probes  are  used  in  plasma. 

For  the  purpose  of  measuring  the  plasma  magnetie  field  ‘b’,  the  only  signal  of  interest  is 
that  induced  in  the  sensing  coil.  To  remove  unwanted  capacitive  pickup,  we  consider 
putting  two  separate  sensing  eoils  next  to  each  other  but  geometrieally  orientated  180 
degrees  with  respect  to  the  other.  The  two  signals  are  then  combined  through  a  Center- 
Tapped  Transformer  (CTT)  such  that  the  eapacitive  voltage  signal  is  subtracted  out  when 
the  leads  from  the  two  probes  are  conneeted  to  opposite  ends  of  the  CTT.  The  two  probe 
orientation  is  shown  in  Figure  4.4. 
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Figure  4.4.  Two  B-dot  probes  orientated  180  degrees  with  respeet  to  one  another; 
signal  is  subtracted  through  CTT. 

The  CTT  used  with  these  probes  is  a  Mini-Cireuits  T16-6T+,  whieh  has  a  flat  frequeney 
response  to  75  MHz.  Additionally,  it  is  an  unbalanced  to  50  fl  balanced  load  transformer 
with  an  impedanee  ratio  of  16  [7].  To  further  reduce  unwanted  noise  and  piekup,  a 
differential  voltage  probe  is  used.  The  differential  voltage  probe  used  in  this  setup  is  the 
Tektronix  P6246,  whieh  has  a  400  MHz  bandwidth,  <  1  pF  input  capaeitance,  200  kQ 
input  resistanee,  and  a  60  dB  eommon  mode  rejeetion  ratio  (emrr);  the  b-dot 
probe/CTT/differential  voltage  probe  sehematic  is  shown  in  Figure  4.5. 


65 


Probe  1 


Probe  2 


Q 


V  =  v.  +v 

^1  ^  ind  '  cap 


0 


v^=- 


-V.  +v 

ind  cap 


0 


> 


Secondary 


O 


V  =v  -  V  =2V 

V  tot  '2  ind 


Primary 


O 


To  Differential  Voltage  Probe  and  Oscilloscope 


Figure  4.5.  Center-Tapped  Transformer  with  b-dot  probes  and  differential  voltage  probe. 


4.1.4  Probe  Characterization 

4. 1.4.1  Selection 

With  the  knowledge  that  the  leads  from  the  magnetic  induction  probe  lowers  the  overall 
frequency  response,  it  was  desirable  to  design  a  coil  with  significantly  lower  inductance 
such  that  when  the  connection  wires  to  the  coil  are  made,  allowance  for  use  at  13.56  MHz 
(desired  operating  frequency)  is  accomplished.  With  this  in  mind  and  the  desire  for  good 
probe  sensitivity,  a  surface  mount  inductor  was  selected.  This  also  allowed 
miniaturization  of  the  sensing  coil  while  providing  a  larger  number  of  turns  for  a  small 
cross  sectional  area. 


The  inductor  selected  was  a  Vishay  Dale  high  frequency  surface  mount  inductor,  P/N 
IMC1008ERR39J.  This  surface  mount  inductor  has  20  turns  wound  on  a  rectangular 
bobbin  measuring  0.83  mm  x  1.29  mm.  Physically,  the  inductor  is  2.5  mm  x  2.0  mm  x 
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1.6  mm  with  a  ceramic  core  and  solder-able  end  connector  plates.  The  specified 
inductance  is  390  nH  with  a  self-resonant  frequency  of  530  MHz  [8]. 


4. 1.4. 2  Frequency  Response 

The  most  important  consideration  when  using  a  magnetic  induction  probe  is  identifying 
the  resonant  frequency  of  the  entire  measuring  circuit.  This  is  based  upon  many  factors 
beyond  just  the  coil  inductance  and  the  parasitic  capacitance  between  successive  turns  of 
wire  as  shown  in  Figure  4.6.  They  include  effects  such  as  transmission  line  lengths,  the 
types  of  cable  used,  and  the  circuit’s  termination  impedance.  An  impedance  analyzer 
quantified  these  effects  on  the  probe's  frequency  response. 


I^series  (SHldll) 

Figure  4.6.  Coil  inductance  and  parasitic  capacitance  that  will 
determine  probe  /  coil  resonance  frequency. 

The  impedance  analyzer  used  in  this  work  is  the  Agilent  4294A  Precision  Impedance 
Analyzer.  The  analyzer  can  be  swept  over  any  specified  frequency  range  from  40  Hz  - 
110  MHz  with  a  user-defined  number  of  points  from  2-801  between  the  frequency 
intervals.  Sweeps  can  be  specified  to  be  linear  or  logarithmic  and  typically  take  <  1  s  [9]. 
The  16047E  test  fixture  was  utilized  with  the  4294A  unit.  This  fixture  is  attached  to  the 
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front  of  the  unit  and  allowed  for  easy  test  and  measurement  of  cables  and  axial  leaded 
components,  such  as  those  utilized  in  this  work.  All  tests  were  conducted  with  a  500  mV 
sinusoidal  signal.  The  analyzer  simultaneously  displays  the  impedance  magnitude  Z  and 
the  phase  delay  6.  While  it  has  the  ability  to  display  a  variety  of  other  parameters  such  as 
resistance,  inductance,  and  capacitance  (both  series  and  parallel  values)  as  a  function  of 
frequency,  the  one  of  primary  interest  in  probe  design  is  the  magnitude  of  the  impedance 
at  the  frequency  under  investigation; /=  13.56  MHz. 

The  frequency  response  of  the  probe  was  investigated  and  found  to  heavily  rely  on  two 
things:  1)  the  length  of  cable  from  the  sensing  coils  to  the  CTT;  and  2)  the  type  of  cabling 
used  (in  this  case  either  a  twisted  shielded  pair  or  coaxial).  The  length  of  cable  from  the 
50  Q  balanced  side  of  the  CTT  to  the  oscilloscope  differential  voltage  probe  was  found  to 
have  a  small  effect  on  the  overall  frequency  response  in  comparison  to  the  cable  length 
on  the  unbalanced  side  of  the  CTT;  this  is  the  result  of  the  CTT  being  a  50  Q  balanced 
load  in  conjunction  with  the  transmission  line  to  the  differential  voltage  probe. 

The  relationship  between  frequency  response  and  transmission  line  length  for  a  twisted 
shielded  pair  of  wires  and  coaxial  wire  was  tested  using  a  single  sensing  coil.  The  results 
are  shown  in  Figure  4.7. 
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Figure  4.7.  Impedance  frequency  response  for  twisted  shielded  pair  and  coaxial  cable. 

Both  cables  exhibit  nearly  identical  resonance  frequeneies  for  the  same  length 
transmission  line.  However,  at  the  resonant  frequeney,  the  coaxial  cable  has  an 
impedance  value  that  is  five  times  smaller  than  the  twisted  shielded  pair.  For  the  intent 
of  this  probe,  the  frequency  of  interest  is  13.56  MHz,  so  that  if  the  transmission  line 
length  is  kept  to  less  than  1  meter,  the  type  of  cable  chosen  should  not  have  much  effect 
on  the  passing  signal.  It  is  prudent  to  note  that  probes  used  in  environments  such  as 
plasma  discharge  will  be  subjeet  to  harmonic  signals  as  well  as  the  fundamental. 
Therefore,  depending  upon  the  application,  higher  or  lower  impedance  at  frequeneies 
larger  than  the  fundamental  may  be  desirable. 


In  this  work,  the  twisted  shielded  pair  of  cable  was  chosen  specifieally  for  the 
aforementioned  reason.  Since  13.56  MHz  was  the  desired  frequency  of  study,  higher 
order  harmonics  are  viewed  as  spurious  signals.  Additionally,  the  twisted  shielded  pair 
of  cable  was  used  on  the  primary  winding  side  of  the  CTT  so  that  when  eonnected  to  the 
differential  voltage  probe,  a  differential  measurement  eould  take  plaee  between  two  lines 
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not  directly  connected  with  an  instrument  ground.  The  line  length  effect  on  the  twisted 
shielded  pair  from  the  probe  to  the  CTT  was  then  investigated  (with  the  CTT  now 
inserted  in-line  whereas  the  CTT  was  not  present  previously).  These  results  are  shown  in 
Figure  4.8  and  clearly  illustrate  that  shorter  line  lengths  yield  a  higher  resonant  frequency 
and  allow  the  probe  to  have  a  wider  operating  regime. 


4.1.5  Discussion 

The  final  probe  characteristics  were  selected  based  upon  some  experimental  and 
laboratory  limitations.  First  of  all,  the  2-probe  combination  (180  degree  spatial 
orientation)  was  measured  to  be  approximately  5.0  mm  in  length  and  2.5  mm  in  height. 
The  probe  was  sized  specifically  to  be  close  to  this  dimension  so  that  spatial  ‘h’  profiles 
could  be  resolved  when  probing  the  -  5  cm  diameter  cylindrical  plasma.  This  ensures 
that  the  physical  disturbance  to  the  plasma  is  an  order  of  magnitude  smaller  than  the 
discharge.  The  probe  is  shown  in  Figure  4.9. 
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Figure  4.9.  Two  probe  180  degree  spatial  orientation; 
shown  with  twisted  shielded  leads  to  the  CTT. 


Second,  the  size  of  the  CTT  was  significantly  larger  than  the  sensing  inductors  and  if  they 
were  placed  close  together  (so  as  to  increase  the  circuit  resonant  frequency),  the  probe  / 
CTT  combination  would  become  large  and  disturb  the  magnetic  environment  under 
investigation.  As  such,  the  probe-to-CTT  length  selected  was  0.67  meters  with  a  resonant 
frequency  of  27.54  MHz.  For  the  investigation  of  13.56  MHz,  this  resulted  in  an 
impedance  of  approximately  10  Q’s.  The  calibration  factor  nA,  based  off  manufacturer 
specifications  was  21.27  mm  . 

The  method  of  using  an  impedance  analyzer  eliminated  many  of  the  previous  difficulties 
in  magnetic  induction  probe  design.  Previously,  discrete  data  points  would  be  taken  over 
a  large  frequency  range  rather  than  swept  over  the  range  of  interest;  this  allows  distinct 
features  in  the  probe  /  circuit  to  be  readily  identified.  It  is  not  only  a  more  time  efficient 
method  of  probe  characterization  but  it  also  allows  the  user  to  measure  the  entire  circuit 
as  a  function  of  frequency.  One  can  clearly  see  the  effect  transmission  line  length  has  as 
well  as  cable  type  on  the  circuit  resonance. 
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For  example,  connecting  a  1  meter  cable  to  a  390  nH  inductor  with  a  self  resonant 
frequency  of  530  MHz,  quickly  drops  the  circuit  resonance  to  17  MHz,  more  than  an 
order  of  magnitude.  The  only  way  to  circumvent  this  problem  is  to  move  a  balanced  - 
unbalanced  CTT  closer  to  the  inductor  (probe).  However,  for  the  13.56  MHz  design 
point  in  this  work,  we  were  fortuitous  that  transmission  line  lengths  did  not  cause 
significant  concern.  While  this  work  does  approach  the  limitation  of  small  physical  size 
(small  A),  with  good  signal  strength  (nA),  design  for  much  higher  frequencies  may  prove 
difficult.  Although,  if  higher  frequency  operation  is  desired,  one  technique  recently 
investigated  has  been  that  of  frequency  mixing  where  a  high  signal  frequency  is  mixed 
with  a  local  oscillator  (LO)  resulting  in  a  more  manageable  intermediate  frequency  (IF) 
[10].  This  method  proved  unnecessary  for  the  work  done  here,  though  a  comparison 
between  the  two  methods  is  probably  warranted. 

4.1.6  Probe  Analysis  in  Plasma 

As  previously  mentioned,  often  times  in  work  with  radio  frequency  plasma  discharge, 
there  exists  the  potential  for  probes  to  pickup  ‘noise’  or  spurious  signal  contributions 
from  harmonics  of  the  driving  frequency.  In  the  case  of  plasma  driven  at  13.56  MHz, 
harmonics  of  this  fundamental  can  be  found  at  27.12  MHz,  40.68  MHz,  54.24  MHz,  etc. 
These  are  contributions  to  the  primary  signal  that  are  not  filtered  out  by  proper  probe 
characterization  and  not  necessarily  due  to  plasma  oscillations.  To  illustrate  this  point, 
we  consider  a  raw  magnetic  probe  signal  obtained  from  the  oscillating  plasma  fields 
(driven  at  500  Watts  with  a  900  Gauss  applied  static  magnetic  field)  as  shown  in  Figure 
4.10.  Additionally,  we  consider  the  Fast-Fourier  Transform  (FFT)  of  the  raw  signal  in 
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order  to  view  eontributions  due  to  harmonics;  the  FFT  in  Figure  4.1 1  was  calculated  as 


N-\ 


waveppj.[n]  =  Y^wave^^.g[k]^e 

k=0 


2m-knl  N 


(4.9) 


where  N  =  25,000  is  the  total  number  or  data  points  sampled  at  the  interval  At  =  1 .6  x  10 
''  s.  In  order  to  accurately  report  what  the  fundamental  (13.56  MFIz)  contribution  is,  a 
Gaussian  fit  was  then  applied  to  the  FFT  in  the  vicinity  of  the  fundamental  and  the  peak 
value  was  interpolated  with  the  amplitude  calculated  as 
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Immediately  evident  from  Figure  4.10  is  that  the  signal  obtained  is  not  a  pure  sinusoid.  It 
has  slight  distortions  due  to  frequency  harmonic  contributions.  Similarly,  Figure  4.1 1 
shows  clear  contributions  at  the  and  5*  harmonics;  40.68  MHz  and  67.80  MHz 
respectively.  In  this  case,  what  appears  to  be  a  greater  than  4  Volt  amplitude  signal  from 
the  raw  data  is  actually  a  3.43  Volt  amplitude  signal  when  the  FFT  is  analyzed  according 
to  (4.9)  and  (4.10);  or  about  22.5%  signal  measurement  error  due  to  frequency 
harmonics. 

For  a  more  demonstrative  example  of  frequency  harmonic  contribution,  we  consider  a 
second  raw  data  trace  obtained  in  a  500  Watt,  600  Gauss  plasma  shown  in  Figure  4.12. 
The  signal  appears  much  more  distorted  due  to  harmonic  contributions  and  when  the  FFT 
(Figure  4.13)  of  the  raw  data  is  analyzed,  contributions  due  to  the  2“‘*,  4*,  5*,  6*,  and 

7*  harmonics  are  found.  Again,  fitting  a  Gaussian  to  the  fundamental,  we  can  calculate 
the  signal  amplitude  to  be  0.42  Volts  where  upon  inspection  the  amplitude  may  be 
(~  0.90+ ~  0.50) 

reported  as  ^ - -  «  0.70  Volts  ;  or  about  a  67%  signal  measurement  error  if 

harmonic  contribution  is  ignored. 
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Figure  4.12.  Raw  Voltage  Data  Trace. 


Obviously,  in  order  to  properly  interpret  results  obtained  from  magnetic  probes  in 
plasma,  it  is  necessary  to  remove  any  spurious  contributions  due  to  frequency  harmonics. 
The  most  accurate  method  to  accomplishing  this  is  to  record  the  raw  signal,  perform  a 
FFT,  and  take  the  peak  of  a  Gaussian  fit  to  the  area  of  interest.  Although,  this  can  often 
be  more  time  consuming  than  utilizing  an  RC  integrating  circuit  before  data  acquisition 
in  order  to  directly  record  the  presumed  magnetic  field  amplitude,  it  is  the  more  accurate 
approach  and  allows  the  user  to  see  harmonic  contributions  to  the  signal  that  will 
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introduce  error.  The  routine  used  to  analyze  all  magnetic  probe  data  in  this  work  was 
performed  and  coded  in  IGOR  and  ean  be  found  in  Appendix  C. 

4.2  RF  Compensated  Langmuir  Probe  and  Microwave  Interferometer 
Metallic  probes  in  plasma  have  been  used  for  deeades  as  one  of  the  most  readily 
implemented  diagnosties  of  obtaining  plasma  parameters  such  as  densities,  temperatures, 
and  potentials.  They  provide  a  method  toward  obtaining  spatially  resolved  measurements 
of  plasma  discharge;  however,  they  are  also  an  intrusive  diagnostic  to  the  plasma 
environment  if  care  is  not  taken  to  minimize  this  effect.  The  classieal  theory  of  Langmuir 
probes  extends  back  to  the  1920’s  although  significant  advancements  in  theory  have  been 
adapted  to  plasma  regions  where  classical  theory  is  no  longer  valid,  i.e.,  interpretation 
and  use  in  strong  magnetic  fields,  impact  of  sheath  thiekness  relative  to  probe  size 
(Lafromboise),  and  most  importantly  for  this  work,  the  implementation  for  use  in  a  RF 
plasma  environment.  This  section  will  discuss  the  important  results  of  construeting  a  RF 
compensated  Langmuir  probe,  whieh  inelude  characterization  of  the  probe’s  cireuit 
impedance  and  ealibration  against  a  90  GHz  mierowave  interferometer. 

4.2.1  RF  Compensation 

RF  compensated  Langmuir  probes  were  constructed  in-house  with  four  noteh  filters;  two 
at  the  fundamental  and  two  at  the  second  harmonic.  Additionally,  a  1  nF  capacitor  was 
placed  between  the  measuring  probe  and  eompensation  electrode  as  shown  in  Figure 
4.14. 
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Figure  4.14.  Schematic  of  Langmuir  Probe  Notch  Filters  and  Compensation  Circuitry. 

The  compensation  electrode  allows  the  probe  tip  to  follow  the  RF  fluctuations  in  the 
potential  while  the  filters  provide  a  large  impedance  for  the  AC  voltage  component  with 
respect  to  the  plasma  sheath  impedance  [11,  12].  The  rationale  for  this  is  that  the  ratio  of 
the  sheath-probe  impedance  should  be  made  as  small  as  possible  so  that  the  overall 
circuit  (probe  and  sheath)  acts  as  a  voltage  divider  and  the  RF  voltage  component  is 
primarily  across  the  large  probe  impedance;  since  the  DC  current  is  the  quantity  of 
interest  when  obtaining  the  IV  characteristic.  The  relationship  that  needs  to  be  satisfied 
when  performing  RF  Langmuir  probe  measurements  is  given  by 


^ probe  I  '^sheath  \ 


(4.11) 


This  is  the  simplified  definition  when  designing  RF  probe  circuitry.  A  more  thorough 
analysis  is  given  in  [1 1,  12]  but  the  result  of  minimizing  the  RF  potential  distortion  to  the 
sheath  can  be  accomplished  by  following  (4.1 1).  When  considering  inductive  plasma  this 
requires  that  the  compensation  electrode  be  on  the  order  of  a  few  nF’s  while  the  choke 
impedance  be  »  1  kQ  [12-14],  This  is  because  the  sheath  resistance  can  be  defined  as 
[12] 

kT 

Rsbeo,b=^^  (4.12) 

which  in  practice  is  on  the  order  of  a  few  kQ’s.  Therefore,  the  probe’s  circuit  or  notch 
filters  must  have  an  impedance  much  larger  than  a  few  kQ’s. 
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The  filters  used  in  the  construction  of  the  probe  are  high  frequency  inductors  from  Vishay 
Dale  where  the  series  inductance  and  parallel  parasitic  capacitance  between  the  number 
of  successive  turns  results  in  the  resonance  curve  for  each  notch  filter.  Due  to  the  fact 
that  each  inductor  is  slightly  different  in  the  values  of  series  inductance  and  parallel 
capacitance,  the  resonant  frequency  for  each  individual  inductor  will  vary.  Additionally, 
the  manufacturer’s  listed  self  resonant  frequency  is  not  always  exact  due  to  non-testing  of 
every  inductor.  However,  for  our  purpose,  this  was  very  important  and  consequently, 
hundreds  of  inductors  impedance  curves  were  analyzed.  To  do  this  quickly,  efficiently, 
and  accurately,  frequency  sweeps  were  again  performed  with  the  Agilent  4294A 
Precision  Impedance  Analyzer;  the  same  one  used  in  the  calibration  of  the  b-dot  probes. 
This  made  it  possible  to  rapidly  test  many  different  inductors  to  determine  which  would 
yield  the  highest  impedance  at  the  fundamental  and  second  harmonic  frequencies.  The 
frequency  response  of  the  four  notch  filters  in  series  is  shown  in  Figure  4.15. 


Figure  4.15.  Impedance  Curve  for  Four  Series  Notch  Filters. 
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There  should  be  caution  that  in  actual  probe  implementation  the  impedance  curves  in 
Figure  4.15  are  higher  than  the  overall  circuitry  impedance  response.  Once  the 
connection  wires  and  the  compensation  electrode  are  connected,  the  impedance  values 
here  drop  by  nearly  an  order  of  magnitude  due  to  shunt  capacitances  and  line  length 
effects.  This  can  often  be  a  limiting  factor  on  the  accuracy  of  Langmuir  probe 
measurements  in  RF.  For  3  eV  plasma  temperature  with  1.0  mA  electron  current,  the 
sheath  resistance  according  to  (4.12)  is  approximately  3  kQ’s.  Therefore,  at  13.56  MHz, 
where  the  total  circuit  impedance  (shown  in  Figure  4.16)  measured  is  85.38  kQ,  the  best 
accuracy  due  to  estimated  sheath  impedances  would  be  about  5%,  or  +  0.15  eV  for  a  3 
eV  temperature.  However,  this  accuracy  will  vary  in  accordance  with  (4.12).  It  is 
noteworthy  to  point  out  that  to  the  best  of  the  author’s  knowledge,  impedance  traces  for 
circuit  filters  such  as  those  shown  in  Figure  4.15  have  been  produced  in  the  literature; 
however,  it  is  the  first  time  that  the  entire  measuring  circuit’s  frequency  response  has 
been  identified. 


Frequency  [MF[z] 

Figure  4.16.  Overall  Langmuir  probe  circuit  frequency 
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The  resonances  in  Figure  4.16  are  for  the  same  inductors  used  in  producing  the  series 
resonance  curves  of  Figure  4.15,  the  only  difference  being  that  the  curves  are  distorted 
and  overall  circuit  impedance  lowered  due  to  line  and  shunt  capacitances.  A  basic  rule  to 
keep  in  mind  is:  capacitors  ‘pass’  high  frequency  signals  and  inductors  ‘pass’  low 
frequency  signals.  This  is  similar  to  the  RF  chokes  where  the  inductor’s  parasitic 
capacitance  in  parallel  with  the  coil  inductance  acts  as  a  bandpass  or  notch  filter.  The 
relations  discussed  here  are  simple  enough  in  theory,  though  complex  when  actual 
inductor  and  capacitor  values  are  considered  with  line  lengths  and  series  /  parallel 
combinations.  The  use  of  the  impedance  analyzer  greatly  reduced  this  difficulty  in  the 
probe  circuit  characterization. 

4.2.2  Obtaining  the  RF IV  Characteristic 

In  order  to  take  actual  plasma  measurements,  the  probe  tip  was  swept  from  ±100  Volts 
over  200  ms  and  the  DC  component  of  the  current  was  measured  across  a  1  kQ  resistor. 
The  measuring  circuitry  consisted  of  a  Kepco  Bi-polar  voltage  power  supply,  while  the 
current  was  measured  with  a  Tektronix  P5205  100  MHz  high  voltage  differential  probe  to 
minimize  any  common  mode  signal.  The  probe  tip  and  compensation  electrode  were 
made  from  0.009”  diameter  (0.02286  cm)  tungsten  wire.  A  typical  probe  tip  was  3.0  mm 
in  length  while  the  compensation  electrode  was  wound  20  times  on  a  0.025”  OD  (0.0635 
cm)  alumina  tube  spanning  6.0  mm.  This  yielded  an  effective  probe  surface  area  of  2.91 
mm  and  a  compensation  electrode  surface  area  of  18.87  mm  ;  such  that  the 
compensation  electrode’s  surface  area  was  approximately  6.5  times  larger  than  the 
probe’s  surface  area.  The  electrical  schematic  is  shown  in  Figure  4.17. 
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Figure  4.17.  Electrical  schematic  of  Langmuir  probe  data  acquisition. 

A  representative  example  of  two  raw  fV  traees  along  with  aeeompanying  eubie  spline  fits 
are  shown  in  Figure  4.18  and  Figure  4.19.  The  eubie  spline  fit  is  1000  data  points 
approximating  the  25,000  obtained  for  a  raw  oseilloseope  traee.  This  was  done  for  two 
reasons:  1)  to  remove  some  low-level  RF  noise  in  the  trace;  and  2)  to  reduce  the  data  set 
to  a  more  manageable  number  of  points.  As  the  data  shows,  no  relevant  information  is 
lost  in  performing  the  cubic  spline  fit.  The  data  analysis  follows  [15]  for  the  Lafromboise 
method  and  a  complete  representative  analysis  can  be  found  in  Appendix  D.  A  difficulty 
in  analyzing  the  RF  IV  characteristic  is  ambiguity  in  identifying  the  plasma  potential. 
Typically,  the  difference  between  the  floating  and  plasma  potential  will  be  several  times 
higher  than  the  electron  temperature  due  to  differences  between  treating  the  thick  vs.  thin 
sheath  scenarios.  However,  most  IV  characteristics  analyzed  in  this  work  did  not  exhibit 
this  effect  and  consequently  the  plasma  potential  is  most  likely  higher  than  what  is 
reported.  This  will  result  in  an  error  for  the  electron  temperature.  The  value  of  the 
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density  reported  will  only  be  affected  by  the  sheath  to  probe  impedance  ratio  with  the 
associated  error  from  the  electron  temperature  calculation. 
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Figure  4.18.  Langmuir  Probe  IV  trace 
47  seem;  10  mTorr;  100  Watt;  900  Gauss;  6.35  cm  OD 


Figure  4.19.  Langmuir  Probe  IV  trace 
47  seem;  10  mTorr;  500  Watt;  900  Gauss;  6.35  cm  OD 

In  Chapter  2,  the  structure  of  helicon  waves  (br,  be,  and  bz  fields)  changed  with  the 
anisotropic  nature  of  the  plasma  cavity  based  on  the  radial  density  distribution. 

Therefore,  the  primary  focus  of  Langmuir  probe  measurements  will  be  based  off  of  the 
radial  and  axial  density  distributions.  The  electron  temperature  obtained  from  the  IV 
characteristic  can  often  be  a  very  subjective  process;  however,  results  for  the  electron 
temperature  will  later  be  reported  as  the  role  of  the  axial  temperature  distribution  plays  an 
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important  part  in  the  axial  plasma  density  distribution.  Sinee  the  density  is  of  immediate 
importance,  it  is  necessary  to  verify  that  the  density  measurements  obtained  from  the  IV 
characteristic  are  validated  against  another  method.  For  this,  a  90  GHz  microwave 
interferometer  was  used. 

4.2.3  Microwave  Interferometer 

Microwave  interferometry  measurements  have  previously  been  used  to  calibrate 
Langmuir  probe  results  [16,  17].  The  diagnostic  has  both  its  advantages  and 
disadvantages  with  respect  to  electrostatic  probes.  While  it  is  a  non-intrusive  diagnostic 
in  the  sense  that  no  perturbation  to  the  plasma  environment  occurs,  this  also  means  that 
no  spatially  resolved  plasma  quantity  is  obtained.  The  result  is  a  radially  integrated 
electron  density.  For  collisionless  plasma  where  the  collision  frequencies  are  much 
less  than  the  electron  plasma  frequency  ,  the  electron  density  is  given  by 

(4.13) 

where  is  the  phase  change  in  the  path  length  7’  of  the  interferometer.  A-o  is  the 

wavelength  of  the  oscillator  frequency,  in  this  case  for /=  90  GHz,  A,o  =  3.33  mm. 
Therefore,  the  two  quantities  that  must  be  determined  are  the  path  length  7’  and  the  phase 
change  of  the  wave  which  passes  through  the  plasma  .  Each  individual 

interferometer  can  have  a  unique  method  toward  properly  calibrating  its  system  such  that 
thermal  and  mechanic  vibrations  are  quantified  or  calibrated  to  minimize  error  or  noise  in 
the  phase  change.  For  this  particular  interferometer,  a  unique  calibration  analysis  was 
developed  so  that  each  branch  (reference  arm  and  plasma  arm)  could  individually  be 
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characterized  in  terms  of  the  signal  amplitude  and  phase  ehange  assoeiated  with  eaeh 
arm.  A  sehematie  of  the  interferometer  setup  is  given  in  Figure  4.20  and  the  technique 
used  to  ealibrate  the  microwave  interferometer  is  given  in  Appendix  E. 
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Figure  4.20.  Schematic  of  the  millimeter  I/Q  Phase  Bridge  Interferometer. 


The  relevant  expression  for  the  phase  ehange  due  to  the  plasma  is  then  given  by 


j  |y0cos(^g)-sin(0,) 

<  - ^ -  > 

cos(^/)  +  y0sin(00) 


(4.14) 


where 


A  V  -V 

13=^'^  OUT,  ^OFF, 

^OUTq  ~^OFFg 

and  Aq  ,  Ai  are  the  sinusoidal  signal  amplitudes  of  each  arm  and  the  voltages  are  the 
output  and  offset  voltages  of  eaeh  arm  (the  details  of  the  definitions  are  also  found  in 
Appendix  E).  Similarly,  the  quantities  0q  and  0i  are  the  phases  of  each  arm  during 
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testing,  which  due  to  manufacturing  tolerances,  are  not  necessarily  taken  at  the  zero 
location.  The  results  show  that  for  an  ideal  interferometer,  where  there  is  no  offset 
voltage,  the  phase  changes  due  to  vibrations,  thermal  effects,  or  small  changes  in  signal 
path  length  are  zero,  and  the  signal  amplitude  within  each  arm  is  the  same,  i.e., 

V  =V  =0  =0  =0 

Aq=A, 

then  the  expression  for  the  plasma  phase  change  reduces  to 

(4.16) 

[  ^OUTg  j 

which  is  the  ratio  of  the  voltage  deflection  in  each  arm  due  to  the  presence  of  plasma. 

This  is  of  course  the  description  for  an  ideal  situation  when  in  practice  and  for  the 
purposes  of  this  experiment,  the  technique  in  Appendix  E  according  to  (4.14)  and  (4.15) 
were  followed. 

4.2.4  Testing  and  Calibration 

Calibration  was  conducted  on  a  500  Watt  inductive  plasma  discharge  with  no  magnetic 
field  while  ten  consecutive  IV  characteristics  were  measured  against  the  microwave 
interferometry  phase  measurement.  The  average  density  value  from  the  Langmuir  probe 
analysis  was  9.97  x  10  m'  while  the  average  phase  deflection  was  7.9  degrees  or  0.14 
radians  resulting  in  a  chord  length  averaged  density  of  2.94  x  10  m'  ;  these 
measurements  yielded  a  34  %  difference  between  average  density  values  and  a 
calibration  factor  of  2.95  for  density  measurements.  Again,  the  densities  in  both  cases 
were  obtained  through  the  analysis  found  in  Appendix  D  and  Appendix  E.  Typical 
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microwave  interferometry  measurements  consisted  of  measuring  the  DC  voltage 
deflection  in  each  arm  of  the  interferometer  as  shown  in  Figure  4.21. 
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Figure  4.21.  Calibration  deflections  in  the  I  and  Q  branches  of  the  microwave  interferometer. 


4.3  Error  Analysis 

The  data  presented  in  the  following  chapters  will  incorporated  an  error  analysis  which 
will  focus  on  two  primary  reported  measurements,  plasma  density  n  and  electron 
temperature  Te.  The  uncertainty  in  these  measurements  will  be  discussed  in  terms  of  the 
dominant  sources  of  error. 

The  uncertainty  in  electron  temperature  will  be  due  to  the  sum  of  the  standard  deviation 
of  10  consecutive  IV  characteristic  analysis’  plus  the  error  due  to  the  ratio  of  the 
estimated  sheath  impedance  to  the  measured  probe  impedance  i.e., 

AT  =(7^  sheath/  (A  1  y') 

^^e  '-’lOtraces^  /7 

/  probe 
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An  example  of  10  conseeutive  Te  measurements  with  corresponding  error  bars  calculated 
from  (4.17)  is  shown  in  Figure  4.22. 


Figure  4.22.  10  Tp  values  for  the  correspondiug  IV  characteristics  with  error  bars  shown. 

The  electron  temperature  is  then  reported  according  to 

T,=T^{avg)  +  M^{avg)  (4.18) 

For  the  example  case  considered  here:  =3.5  ±0.6  eV .  While  the  uncertainty  to  the 

standard  deviation  does  not  change  as  it  is  calculated  from  the  10  Te  values,  the 
uncertainty  due  to  sheath  /  probe  impedance  ratios  does  because  the  sheath  impedance  is 
not  a  fixed  value  and  is  calculated  according  to  [12] 

kT 

(4-19) 

eL 


varying  with  the  collected  saturation  current  at  each  data  point. 
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The  density  error  analysis  was  calculated  in  much  the  same  way  although  we  also  include 
the  error  present  in  the  microwave  interferometer  which  manifests  in  the  calibration 
factor  of  2.95.  The  error  in  density  measurements  can  be  calculated  according  to 

=  ^Laces  +  (//wave)  (4.20) 

/  probe 

The  standard  deviation  of  the  10  traces  is  a  straightforward  calculation  and  the  impedance 
ratio  error  is  the  same  percent  error  calculated  for  the  electron  temperature.  The 
uncertainty  in  the  calibration  factor  can  be  obtained  from  the  voltage  data  in  Figure  4.21. 
The  I-branch  and  Q-branch  have  cr^  values  of  0.15  mV  and  0.09  mV,  respectively.  This 
translates  to  an  error  in  the  phase  measurement  of  4.3  degrees  and  an  error  in  the 
calibration  factor  of  1.6.  Therefore,  the  calibration  factor  of  2.95  becomes  2.95  ±  1.6. 

The  error  in  density  can  then  be  calculated  utilizing  (4.20).  The  density  values  with 
corresponding  error  bars  are  shown  in  Figure  4.23. 


Figure  4.23.  10  n  values  for  the  correspondiug  IV  characteristics  with  error  bars  shown. 
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17  3 

For  this  example,  the  density  is  reported  as  2.94  ±  1.4  x  10  m'  .  There  exist  other 
sources  of  error  such  as  uncertainty  in  the  exact  probe  tip  length  of  3.0  mm;  however, 
these  will  likely  be  much  less  than  those  due  to  the  sheath  /  probe  compensation  and 
standard  deviations  which  are  present.  In  summary,  the  average  error  in  the  density 
measurements  was  46  %  where  approximately  14  of  this  uncertainty  is  due  to  calibration 
against  the  90  GHz  microwave  interferometer.  The  average  error  in  the  electron 
temperature  measurements  was  a  more  modest  18  %;  however  this  is  likely  the 
consequence  of  not  calibrating  against  a  well  known  electron  temperature  diagnostic  or 
source. 


4.4.  Chapter  Summary 

This  chapter  described  the  impedance-frequency  characterization  of  magnetic  induction 
probes.  A  novel  technique  was  developed  to  calibrate  the  entire  probe’s  frequency 
response;  up  to  100  MHz  using  an  impedance  analyzer.  The  effects  of  line-lengths  and 
cable-types  were  investigated  were  the  key  results  showed  that  longer  line-lengths 
between  the  probe  and  data  acquisition  device  significantly  decrease  the  overall  probe 
circuit  frequency  response;  limiting  the  probe’s  useful  frequency  range.  A  proper 
technique  to  ensure  the  correct  induced  voltage  in  the  sensing  probe  was  also  developed 
such  that  plasma  capacitive  voltage  was  not  being  measured.  Finally,  we  illustrated  the 
necessity  of  performing  FFT’s  on  the  raw  data  signals  in  order  to  remove  harmonics 
where  errors  were  found  to  be  as  great  as  67%  had  the  correct  data  analysis  through 
FFT’s  not  been  performed. 
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The  same  impedanee-frequency  analysis  was  applied  to  our  RF  eompensated  Langmuir 
probe  sueh  that  we  ensured  a  mueh  higher  overall  probe  impedanee  with  respeet  to  the 
plasma  sheath  impedanee  at  the  driving  frequeney  and  2"^^  harmonic.  Additionally,  the 
compensation  electrode  allowed  the  probe  tip  to  float  with  the  13.56  MHz  oscillating 
plasma  potential.  This  form  of  probe  compensation  provides  a  more  accurate 
measurement  of  the  probe  current  (for  density  calculations)  as  well  as  a  more  accurate 
measurement  of  the  plasma  potential  (for  electron  temperature  calculations)  had  no 
compensation  been  implemented.  Finally,  the  RF  Langmuir  probe  density  measurements 
were  calibrated  against  a  90  GHz  microwave  interferometer  to  obtain  the  chord  averaged 
calibration  factor  of  2.95  ±  1.6. 
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CHAPTER  5:  RESULTS  AND  DISCUSSION 


An  important  aspect  of  helicon  wave  plasma  is  the  parameter  space  of  operation.  Chapter 
2  described  helicon  conditions  dependent  on  frequency/  density  n,  and  applied  magnetic 
field  Bo.  These  quantities  are  intricately  related  through  the  helicon  dispersion  relation, 
which  also  incorporates  some  geometrical  aspects  of  the  experiment  such  as  tube  radius 
a,  antenna  configuration  m,  and  to  a  lesser  extent,  the  antenna  length  /.  In  addition,  there 
is  a  pressure  and  power  aspect  to  the  operational  parameter  space  because  density  is  also 
a  function  of  the  neutral  gas  pressure  and  how  well  the  antenna  couples  or  transfers 
energy  to  the  plasma 

In  light  of  these  many  quantities  it  is  prudent  to  experimentally  determine  the  range  of 
parameters  that  satisfy  helicon  wave  propagation;  or  rather,  determine  the  conditions 
where  the  onset  of  helicon  waves  occur  and  then  to  further  investigate  those  conditions. 
Additionally,  in  accordance  with  Chapter  2,  helicon  waves  are  described  by  radial 
variations  in  the  wave  hr,  be,  and  bz  -  profdes  which  are  dependent  on  radial  density 
profiles  n(r).  Therefore,  this  chapter  will  begin  with  determining  and  discussing  the 
identification  of  helicon  waves  at  various  operating  conditions  and  comparing  the 
measured  and  theoretically  predicted  profiles.  The  chapter  will  progress  to  follow  the 
evolution  of  radial  hr,  be,  and  bz  -  scans  to  two-  and  three-  dimensions  illustrated  by 
contour  maps  of  helicon  waves  over  an  entire  wavelength.  The  results  of  mapping 
helicon  waves  in  three  dimensions  reveals  important  information  regarding  the 
wavelength  of  helicon  waves  and  quantifying  them  as  a  true  3-d  helix  whereas  single  axis 
measurements  only  reveal  what  would  otherwise  appear  to  be  a  2-d  plane  wave. 
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Finally,  this  chapter  will  show  axial  density  profiles  which  correspond  to  the  wave  b-field 
measurements.  In  some  cases,  the  axial  density  measurements  reveal  an  unexpected 
density  peak  downstream  of  the  antenna.  The  explanation  of  this  effect  will  be  shown  to 
be  independent  of  helieon  wave  plasma  and  the  model  provided  is  solely  based  off  of 
performing  an  energy  balanee  within  a  plasma  volume  while  assuming  a  decaying 
(cooling)  electron  temperature  away  from  the  source  (antenna). 

5.1  Identification  of  Helicon  Waves 

Identification  of  helicon  wave  propagation  was  originally  based  upon  the  magnetic  field, 
hr,  be,  and  bz  -  profile  measurements  conducted  with  b-dot  probes  [1].  The  results  were 
then  compared  to  the  hr,  be,  and  bz  -  field  solutions  in  a  cylindrical  plasma  discharge 
subject  to  a  uniform  plasma  density  distribution  [2].  For  a  prescribed  driving  frequency, 
tube  radius,  and  gas,  the  power  and  the  static  magnetic  field  Bo  could  be  varied  (typically 
increased)  until  they  very  nearly  satisfied  the  dispersion  relation.  At  that  point,  the 
internal  measured  magnetic  field  profiles  would  nearly  resemble  the  caleulated  profiles 
and  the  conclusion  was  that  the  dielectric  could  support  the  propagation  of  helicon  waves. 
However,  there  have  been  several  other  methods  commonly  used  to  determine  helicon 
wave  propagation.  These  methods  and  their  accuracy  will  be  briefly  addressed. 

5.1.1  Forward  -  Reflected  Power  Spectrum 

The  first  method  does  not  utilize  internal  plasma  diagnostics.  Identification  is  based  on 
monitoring  the  forward-to-reflected  power  ratios  at  a  given  magnetic  field.  As  the  power 
is  increased,  a  sharp  decline  in  the  reflected  power  level  is  observed  and  helicon  waves 
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are  assumed  to  propagate.  An  example  of  this  situation  is  given  for  a  5  em  cylindrical  m 
=  +1  half  turn  helical  antenna  (20.32  cm  between  end  rings).  The  gas  is  Argon  at  47 
seem  (~  10  mTorr  neutral  gas  pressure).  The  forward  and  reflected  powers  were  recorded 
at  different  input  powers  and  magnetic  field  strengths  and  the  results  are  shown  in  Figure 
5.1. 


Figure  5.1.  Forward  vs.  Reflected  Power  at  47  seem;  10  mTorr  All  pressure. 


During  each  condition,  the  reflected  power  begins  as  a  minimum  and  then  slowly 
increases  as  the  input  power  is  raised  until  some  point  where  there  is  a  ‘sharp’  decline  in 
the  reflected  power.  At  this  point,  the  power  is  said  to  couple  more  efficiently  to  the 
plasma  and  could  be  used  as  a  method  to  identify  the  onset  of  helicon  wave  propagation. 
Also,  Figure  5.1  shows  that  as  the  magnetic  field  is  increased,  the  decline  in  reflected 
power  occurs  at  a  lower  input  power  level.  The  transitional  power  value  is  equivalent  to 
analyzing  the  dispersion  relation  for  helicon  wave  propagation  [3-6]  at  a  fixed  frequency 
and  magnetic  field;  the  relation  is  satisfied  once  the  density  reaches  a  critical  value. 

Since  density  is  proportional  to  input  power  [5],  once  the  forward  power  level  generates 
this  critical  density  value  for  a  given  magnetic  field,  helicon  waves  can  propagate. 
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Utilizing  the  power  speetrum  as  a  method  to  identify  helicon  wave  propagation  can  be 
vague  and  inaccurate.  Since  power  loss  can  come  in  many  forms,  i.e.,  transmission  line 
loss,  radiated  loss,  antenna  coupling  efficiency,  etc,  using  the  power  spectrum  as  a  tool  to 
identify  helicon  waves  can  be  ambiguous.  Although,  as  a  qualitative  observation,  some 
revealing  information  can  be  gained,  i.e.,  there  exists  a  magnetic  field  strength  threshold 
where  little  is  gained  by  further  increasing  the  power.  This  is  shown  by  the  similarities  in 
the  900  Gauss  and  1050  Gauss  curves  of  Figure  5.1.  The  forward  power  threshold  is 
nearly  identical  in  each  case  suggesting  strengthening  the  magnetic  field  beyond  900 
Gauss  for  these  conditions  does  not  generate  a  transition  to  helicon  mode  at  lower 
forward  powers. 

Finally,  the  rise  in  reflected  power  can  also  be  attributed  to  the  initial  settings  of  the 
capacitors  within  the  auto-match  network  and  the  tolerances  set  within  the  feedback 
circuit.  For  example,  referring  to  the  circuit  of  the  matching  network  with  plasma  as 
given  in  Chapter  3  and  repeated  here  in  Figure  5.2. 


Figure  5.2.  RF  power  delivery  system  incorporating  the  auto-match  network  circuit. 
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The  RF  voltage  souree  will  output  power  into  an  equivalent  50  Q  load  at  point  A.  The 
total  applied  power  to  the  matehing  network  and  antenna  will  be  given  by  the  differenee 
between  the  forward  and  reflected  power  measured  at  point  A  according  to 

P  =  P  -P  t5 

TOT  ^  FWD  ^  RFL  V-^  • 


Also,  the  impedance  in  each  “branch”  is  given  by 


z  .  =  z,  ^  +  z,  +z  , 

senes  Load  L  plasma 


Z,  =  Z 


TUNE 


(5.2) 

(5.3) 


For  an  equivalent  circuit  impedance  of 


Z 


eq 


series'^  \\ 
'^series  '^ll 


(5.4) 


Therefore,  regardless  of  the  plasma  impedance,  there  exist  capacitance  values  in  which 
the  equivalent  load  can  be  set  to  50  fl  resistive.  The  restraint  on  the  matching  network 
lies  with  the  dynamic  range  of  the  variable  capacitors.  To  match  the  50  Q  load,  the 
plasma  complex  impedance  may  attempt  to  force  the  capacitors  to  values  outside  their 
working  range.  In  this  case,  the  auto-match  network  searches  for  a  minimum  value  or 
minimum  reflected  level  within  the  feedback  tolerance.  If  no  suitable  values  are  found, 
the  reflected  power  will  continue  to  rise  until  the  plasma  load  again  returns  to  the 
operational  range  of  the  variable  capacitors.  This  explains  the  rise  in  reflected  power  as  a 
consequence  of  the  matching  network  and  not  necessarily  a  mode  transition,  although  the 
power  spectrum  result  will  be  shown  to  closely  coincide  with  the  onset  of  helicon  wave 
fields. 
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5.1.2 Density  ‘Jumps’ 

Another  method  utilized  to  identify  helicon  waves  has  been  through  observation  of 
apparent  density  ‘jumps’  as  the  input  power  is  increased.  This  is  a  sharp  increase  in 
density  from  a  capacitive  or  inductive  mode  to  helicon  wave  mode.  For  the  case  of  900 
Gauss  static  magnetic  field  on  the  same  cylindrical  tube  as  above,  radial  density  profiles 
were  analyzed  for  various  input  powers  as  shown  in  Figure  5.3. 


Radius  [cm] 

Figure  5.3.  Radial  plasma  density  profiles  at  900  Gauss  for  various  input  power 
levels  on  a  5  cm  diameter  quartz  tube  at  13.56  MHz. 


The  density  profiles  and  power  spectrum  analysis  show  helicon  wave  plasma  exhibits  a 
radial  parabolic  profile  whereas  the  inductive  (lower  power)  discharge  maintains  a 
relatively  fiat  density  profile.  As  can  be  seen  for  the  200  Watt  case,  the  ‘small’  Langmuir 
probe  used  to  measure  the  density  destroys  or  interferes  with  the  helicon  wave  from  about 
-1  cm  to  +1  cm  as  the  probe  is  swept  radially  through  the  discharge.  However,  at  higher 
powers,  the  density  is  large  enough  to  overcome  the  presence  of  the  probe  in  order  to 
maintain  helicon  wave  propagation  around  the  probe.  This  also  provides  a 
complementary  explanation  to  the  observation  of  a  density  ‘jump’  which  was  explained 
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in  terms  of  overcoming  parasitic  circuit  losses  [7]  (sharp  drop  in  reflected  power  due  to 
plasma  load).  The  rationale  here  is  that  helicon  waves  would  exist  at  this  condition  if  the 
probe  were  not  present  (or  if  the  capacitors  had  a  wider  dynamic  range).  It  is  the 
presence  of  the  measuring  probe  in  this  case  which  inhibits  the  propagation  of  the  wave 
until  the  plasma  discharge  conditions  are  sufficient  to  overcome  the  probe’s  presence. 


Additionally  observed  in  Figure  5.3,  the  peak  density  appears  to  be  shifted  slightly  left  of 
center  (~  0.5  cm);  this  can  be  attributed  to  several  factors.  First,  the  alignment  of  the 
probe  in  the  center  of  the  quartz  tube  has  some  associated  error  as  this  was  not  given 
significant  importance.  Similarly,  the  quartz  tube  may  not  have  been  perfectly  aligned 
within  the  center  of  the  magnets.  In  order  to  ensure  that  the  ‘peak’  did  not  shift  or 
oscillate,  radial  profiles  at  different  downstream  locations  from  the  antenna  were  taken 
and  shown  to  consistently  be  shifted  to  the  left  of  the  presumed  center  as  shown  in  Figure 
5.4. 
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Figure  5.4.  Downstream  radial  density  profiles  at  500  Watt,  900  Gauss. 
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The  density  jumps  exhibited  in  Figure  5.3  have  been  previously  modeled  off  of  the 
matehing  network  configuration  (L-,  T-,  Ft-),  its  component  values  (capacitances  and 
inductances),  and  transmission  lines  [7].  The  experimental  data  presented  here  validates 
the  model’s  results,  yet  also  demonstrates  the  effect  of  probe  intrusion  at  lower  power 
levels,  which  can  inhibit  helicon  waves. 

Consequently,  in  both  the  case  of  analyzing  the  power  spectrum  and  utilizing  plasma 
density  measurements,  both  methods  can  be  inaccurate  at  identifying  the  propagation  of 
helicon  waves.  Furthermore,  when  starting  the  analysis  of  helicon  waves,  the  only 
derived  measurement  is  that  for  the  internal  magnetic  fields.  Therefore,  the  next  section 
will  address  identifying  helicons  through  the  measurement  of  the  internal  radial  plasma 
wave  fields. 

5.1.3  Radial  -  Profiles 

Helicon  wave  fields  can  have  many  different  structures  as  shown  in  Chapter  2.  These 
different  profiles  are  dependent  on  the  radial  density  profiles.  Therefore,  once  the  radial 
density  distribution  is  measured,  the  wave  profdes  (b-dot  probes)  should  match  the 
theoretical  ones.  The  corollary  of  this  statement  should  also  be  true;  measuring  the 
internal  wave  profiles  i.e.,  the  component  of  the  wave,  the  radial  density  distribution 
can  be  predicted. 

This  section  will  examine  the  radial  magnetic  field  profiles  for  the  bz  wave  component. 
Measurements  are  made  from  pressures  of  0.5  mTorr  -  10  mTorr  on  a  5  cm  quartz  tube 
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with  a  m  =  +1  half  turn  helical  antenna  (20.3  ±0.1  cm  length  between  end  rings)  driven 
at  13.56  MHz.  Magnetic  fields  of  300,  600,  and  900  Gauss  are  tested  at  input  power 
levels  from  100  -  500  Watts.  The  onset  of  helicon  wave  propagation  in  each  instance  is 
examined  and  compared  with  theoretical  profiles  derived  from  Chapter  2  and  correlated 
with  a  radial  density  distribution. 

The  first  test  conditions  were  47  seem  (~  10  mTorr)  Argon  gas,  input  power  300,  400, 
and  500  Watts,  and  peak  magnetic  fields  of  300,  600,  and  900  Gauss.  Those  results  are 
shown  in  Figure  5.5.  For  the  nine  conditions  tested  at  this  pressure,  live  were  found  to  be 
in  excellent  agreement  with  theoretical  profiles.  Those  five  are  indicated  in  Figure  5.5  by 
a  GREEN  border  and  are  listed  here  as 

500  Watts  -  600  &  900  Gauss 
400  Watts  -  600  &  900  Gauss 
300  Watts  -  900  Gauss 

The  ones  identified  with  a  YEEEOW  border  are  radial  scans  where  the  measured  profiles 
near  the  tube  radius  would  seem  to  indicate  the  start  of  a  helicon  bz  profile.  However, 
similar  to  the  200  Watt  density  scan  performed  in  Eigure  5.3,  the  probe’s  presence 
interferes  with  the  wave  propagation  at  these  conditions,  i.e.,  the  wave  field  intensity  is 
too  low  to  overcome  the  presence  of  the  probe.  Einally,  the  RED  border  case, 

300  Watts  -  300  Gauss 

indicates  no  wave  field.  This  is  a  condition  where  the  density  and  magnetic  field  do  not 
sufficiently  create  a  dielectric  medium  to  support  helicon  waves. 


101 


An  alternative  argument  for  the  RED  border  cases  still  propagating  helicon  waves 
considers  the  following; 

Suppose  at  the  condition  300  Watts  -300  Gauss  helicon  waves  are  present.  And  perhaps 
the  data  do  not  show  identifiable  wave  profiles  due  to  the  measuring  b-dot  probe 
exhibiting  insufficient  signal  detection 

The  response  is  helicon  waves  either  do  not  exist  below  some  critical  levels  or  the  wave 
information  carried  by  the  plasma  is  so  small  that  it  is  undetectable  by  the  means 
presented  here.  To  support  this  rationale,  the  results  show  b-dot  probes  which  “pick-up” 
the  start  of  a  helicon  wave  profile  near  the  radius  of  the  cylinder;  however,  the  probe 
interference  destroys  these  signals,  altering  the  dielectric  medium.  Helicon  waves  exist  at 
the  YELLOW  border  conditions  though  the  probe  presence  obscures  this  result. 


Figure  5.5.  profiles  at  10  mTorr  and  various  input  power  and  magnetic  field  strengths. 
Green  borders  indicate  well  matched  with  theory  where  ‘w’  indicates  the  Gaussian  density 
profile  width.  Yellow  borders  indicate  probe  interference  with  helicon  wave  mode.  Red 
borders  indicate  insufficient  conditions  for  wave  propagation. 
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Similar  power  (density)  /  magnetie  field  matriees  are  presented  in  Appendix  F  for  the 
pressures  of  0.5,  2,  4,  6,  and  8  mXorr.  In  all  cases,  as  the  pressure  is  further  decreased, 
the  wave  profiles  only  match  the  theoretical  predictions  at  higher  powers  and  magnetic 
field  strengths.  For  example,  in  the  2  mXorr  analysis,  there  is  only  one  condition  where 
the  wave  profile  nearly  matches  a  theoretical  profile;  500  Watts  -  900  Gauss. 
Consequently,  for  each  of  the  neutral  gas  pressures  tested  where  helicon  wave-like 
structures  were  observed,  the  measured  wave  profdes  were  correlated  (overlaid)  with  the 
theoretical  profiles  and  the  onset  of  helicon  waves  was  determined  by  their  agreement. 


Returning  to  the  GREEN  border  cases,  the  theoretical  curves  based  on  radial  density 
profiles  were  fit  to  the  measured  bz  data.  Xhe  radial  density  profiles  which  identify  the 
conditions  above  are  shown  in  Figure  5.6  where  the  densities  have  been  set  to  a  center 

}Q  0 

value  of  no  =  1  X  10  m  which  closely  approximates  the  measured  centerline  density. 


Figure  5.6.  Radial  density  profiles  corresponding  to  the  data. 

Recalling  Figure  5.3  for  the  measured  radial  density  profiles  at  900  Gauss  and  300,  400, 
and  500  Watts,  the  measured  density  profiles  are  compared  to  the  theoretical  profiles 
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derived  from  the  bz  data.  This  is  shown  in  Figure  5.7.  However,  the  theoretieal  density 
profiles  do  not  neeessarily  identieally  match  those  measured.  The  Gaussian  profiles 
defined  in  Chapter  2  and  shown  in  Figure  5.6  for  w  =  1.0  are  much  narrower  than  those 
measured.  In  fact,  the  best  Gaussian  fit  is  for  w  =  1.7;  also  shown  in  Figure  5.7  for 
comparison  with  the  measured  density  profiles. 
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Figure  5.7.  Theoretical  density  profiles  compared  with  experimentally 


The  differences  in  the  wave  profiles  are  demonstrated  using  the  two  theoretical  density 
profiles  (w  =  1.0  and  w  =  1.7)  and  numerically  solving  for  bz  as  outlined  in  Chapter  2. 
The  differences  between  both  are  illustrated  in  Figure  5.8. 
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Figure  5.8.  Differences  in  calculated  wave  profiles  for  radial  density  profiles  w  =  1.0, 1.7  as 
compared  with  the  experimentally  measured  profiles. 

While  the  theoretieal  eurves  for  w  =  1.0  better  represent  the  experimental  bz  values  as 
opposed  iow  =  1. 7,  the  differenees  in  the  calculated  wave  fields  are  minimal.  An 
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explanation  to  resolve  the  small  difference  may  be  due  to  a  non-vacuum  gap  between 
plasma  and  cylindrical  tube  boundary  when  considering  the  theoretical  analysis.  If  the 
model  were  further  developed  to  include  an  additional  vacuum  layer  between  the 
boundary  and  plasma,  this  may  account  for  the  difference. 

Regardless,  the  results  clearly  show  helicon  wave  theory  which  incorporates  a  non- 
uniform  radial  density  distribution  provides  the  most  accurate  representation  of  the 
measured  plasma  wave  fields.  Additionally,  a  suitable  operational  space  for  helicon 
waves  has  been  determined,  which  will  allow  more  in  depth  study.  The  500  Watt  -  900 
Gauss  operating  condition  at  10  mTorr  neutral  gas  pressure  will  be  used  as  the  control 
condition. 

Each  density  profile  also  provides  a  unique  value  of  kz  for  the  propagating  wave.  This 
value  was  provided  when  the  wave  profiles  and  boundary  conditions  were  solved. 
Calculating  these  value  for  both  the  w  =  1.0  and  w  =  1.7  profiles  yields  a  wavelength 
range  of  =  29.95-33.85  cm  for  the  500  Watt  -  900  Gauss  plasma.  The  expected 
wavelength  (as  will  be  shown)  is  most  likely  an  underestimate.  This  is  due  to  the 
Gaussian  density  profile  not  remaining  constant  in  the  axial  direction.  As  the  wave 
propagates  away  from  the  source,  the  density  profile  will  no  doubt  ‘relax’.  In  fact,  this 
was  already  shown  in  Figure  5.4,  when  the  downstream  radial  density  profiles  were 
examined.  For  each  profile,  one  would  expect  a  different  value  for  the  axial 
wavenumber,  and  hence  the  axial  wavelength.  In  essence,  the  wave  should  be  expected 
to  be  damped  in  the  axial  direction.  This  point  was  already  implied  at  the  end  of  Chapter 
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2  where  the  various  Gaussian  density  profiles  returned  a  different  axial  wavelength; 
however,  the  results  imply  this  oeeurs  in  the  axial  direction  due  to  the  non-uniformity  in 
the  +z  direction.  This  also  explains  some  reported  results  in  which  the  wavelength  in  the 
axial  direction  varies  [8]  which  should  be  expected  from  the  decreasing  axial  density 
where  the  wavelength  will  increase.  Additionally,  the  measured  wavelengths  reported  do 
not  exhibit  a  smooth  increase  in  the  axial  wavelength.  In  fact,  the  wavelength  initially 
decreases  before  smoothly  increasing.  This  has  been  an  unexplained  feature  in  the  axial 
wavelength  to  date;  however,  we  will  later  show  the  initial  decrease  in  axial  wavelength 
is  due  to  a  downstream  rise  in  density.  Regardless,  this  provided  enough  motivation  to 
further  investigate  the  helicon  wavelength  and  its  relation  to  the  plasma  dielectric 
properties  through  axial  measurements. 


5.2  Axial  br,  be,  and  b^  -  Profiles 

The  location  of  the  three  separate  antennae  tested  as  shown  in  Chapter  3  are  repeated  here 
as  shown  in  Figure  5.9. 


Direction  of  Applied  Field  AAA 


Figure  5.9.  Applied  Magnetic  Field  with  Location  and  Lengths  of  Tested  Antennas. 
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It  is  necessary  to  redefine  the  antenna  lengths  in  terms  of  their  full  three  dimensional 
helix.  Where  previous  lengths  have  been  defined  in  terms  of  the  distance  between  the 
antenna  end  rings,  the  full  antenna  length  must  include  the  radial  wrapped  component. 
We  define  this  for  a  full  turn  and  V2  turn  helical  antenna  as  it  relates  to  the  anticipated 
helicon  wavelength  by 


f  =  ^d^  ±(2;ra)^  full  turn 

(5.5) 

h  ~  ±(^a)^  turn  — >  2f 

(5.6) 

where  d  is  the  distance  between  the  end  rings  and  a  is  the  radius  of  the  tube  the  antenna  is 
wrapped  around.  For  the  previous  defined  antenna  lengths  of  20.3  cm,  10.2  cm,  and  5.1 
cm  wrapped  around  a  V2  turn  helical  6.4  cm  OD  quartz  tube,  the  new  antenna  lengths 
according  to  (5.6)  are  22.6  ±0.1  cm,  14.2  ±0.1  cm,  and  1 1.2  ±  0.1  cm  respectively. 
Shorter  distances  between  antenna  end  rings,  where  d,  has  the  effect  of  significantly 
miscalculating  the  anticipated  helicon  wavelength  2^^  if  the  three  dimensional  lengths  are 

not  considered.  For  example,  in  the  case  of  the  5.1  cm  distance  between  end  rings,  the 
true  antenna  length  given  by  (5.6)  is  1 1.2  cm,  or  an  error  of  55  %  between  the  antenna 
lengths. 

5.2.1  Single  Axis  Profiles 

A  routinely  used  method  toward  measuring  the  helicon  wavelength  has  been  to  take  an 
axial  scan  and  report  either  the  phase  difference  between  the  sensing  probe  and  the 
antenna  current  or  to  take  the  length  between  successive  maxima  and  minima  in  the  wave 
amplitude.  The  successive  maxima  amplitude  method  was  repeated  in  this  study; 
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however,  it  is  inaccurate  in  reporting  the  helicon  wavelength  even  though  the  results 
appear  correct  upon  initial  inspection.  For  the  22.6  cm  antenna,  the  results  of  a  single 
axis  scan  for  the  bz  wave  field  at  three  different  radial  ‘x’  values  are  shown  in  Figure 
5.10. 


Figure  5.10.  Single  axis  scan  for  the  wave  fields  taken  at  three  different  radial  locations. 

For  each  scan,  a  different  value  for  the  helicon  wavelength  is  reported,  often  differing  by 
as  much  as  ~  5  cm.  The  same  measurement  for  the  14.2  cm  and  1 1.2  cm  antenna  lengths 
are  shown  in  Figure  5.11  and  Figure  5.12.  In  all  three  cases,  the  single  axis  measured 
two  dimensional  helicon  wavelengths  are  ambiguous.  The  scans  are  out  of  ‘phase’ 
depending  on  where  you  are  performing  the  measurement  and  the  resulting  length  is  seen 
to  vary  by  up  to  50  %.  While  the  scans  do  appear  to  represent  a  wave  with  an  associated 
wavelength,  three  dimensional  scans  of  the  helicon  wave  provide  a  much  more 
meaningful  representation  of  the  true  wavelength  and  the  apparent  phase  difference 
between  successive  locations  is  readily  explained  when  the  three  dimensional  wave  is 
considered. 


no 


Axial  Distance  From  Antenna 

Figure  5.11.  Single  axis  scan  for  the  wave  fields  for  the  14.2  cm  antenna 


Axial  Distance  From  Antenna 

Figure  5.12.  Single  axis  scan  for  the  wave  fields  for  the  11.2  cm  antenna  length. 


5.2.2  Two-  and  Three-  Dimensional  Profiles 

The  scans  at  different  radial  locations  being  out  of  ‘phase’  are  the  result  of  the  two- 
dimensional  representation  of  a  helix  as  shown  in  Figure  5.13.  The  result  for  the  helix 
(which  is  three-dimensional)  is  a  single  wave  with  only  one  associated  wavelength. 


Ill 


Figure  5.13  provides  an  explanation  for  why  previous  results  of  the  helicon  wave  appear 
to  be  ‘phase’  shifted  when  measured  at  different  radial  locations. 


0  degree  shift 


Figure  5.13.  Phase  shift  for  2-d  representation  of  3-d  helix. 


When  wavelengths  have  been  previously  reported  based  off  of  the  two-dimensional 
single  axis  scans  for  the  profdes  [1,  4,  9-13]  they  have  been  reported  based  off  of 
measurements  as  shown  in  Figure  5.10  -  Figure  5.12.  Flowever,  utilizing  the  3-d 
representation  of  an  attenuating  wave  as  modeled  in  Figure  5.14  (in  2-d)  and  Figure  5.15 
(in  3-d)  as  one  that  linearly  decreases  in  the  radial  direction  from  1  cm  to  0.5  cm  over  two 
wavelengths,  the  actual  length  of  this  modeled  wave  can  be  then  calculated  from 
Equation  (5.5)  or  (5.6)  which  will  include  the  azimuthal  component  of  the  wave.  The 
waves  in  Figure  5.14  -  Figure  5.15  are  the  exact  same  waves;  with  Figure  5.14  being  the 
2-d  projection.  This  clearly  illustrates  how  the  2-d  measurement  will  underestimate  the 
full  wave  helix.  The  2-d  wavelength  is  observed  to  be  6  cm  while  the  3-d  wavelength  is 
6.6  cm;  9  %  error. 
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Figure  5.14.  2-d  representation  of  single  phase  helix. 


Figure  5.15.  3-d  helix  represented  hy  the  series  of  2d  phase  shifted  waves. 

Consequently,  to  image  the  entire  wave,  3 -dimensional  profiles  were  made  at  22 
axially  separate  eross-sections.  Each  cross  section  consisted  of  scans  in  the  XY  plane  for 
approximately  75  data  points  per  cross-section  or  about  1650  data  points  per  3- 
dimensional  contour  map.  The  distance  between  data  points  in  a  cross-section  varied 
between  1.2  -  1.4  mm  while  each  axial  cross-section  was  separated  by  approximately  1.5 
cm.  The  data  point  /  node-numbering  scheme  and  contour  plots  were  constructed  in 
Techplot  10.  A  typical  2-dimensional  cross-section  scheme  is  shown  in  Figure  5.16. 
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The  compiled  3 -dimensional  contour  plots  are  shown  in  Figure  5.17  for  the  cases  of 

22.6  cm,  14.2  cm,  and  1 1 .2  cm  antenna  lengths;  as  previously  defined.  Additionally,  2- 
dimensional  (x,z)  contour  plots  in  the  y  =  0  plane  are  shown  in  Figure  5.18  for  b^.  In  the 
case  of  all  3  antenna  lengths,  the  edge  of  the  antenna  is  set  at  z  =  0.  This  is  to  facilitate 
the  determination  of  the  helicon  wavelength.  The  actual  location  of  each  antenna  edge 
with  respect  to  the  applied  magnetic  field  was  given  in  Figure  5.9.  As  illustrated  in 
Figure  5.17  and  Figure  5.18,  the  shortest  antenna  length  does  not  necessarily  correspond 
to  the  shortest  axial  wavelength  helicon.  This  is  likely  due  to  the  14.2  cm  antenna 
propagating  at  a  resonance  with  the  cylindrical  tubes  geometric  radial  dimensions.  The 
antenna  couples  energy  to  the  plasma  more  efficiently  and  exhibits  a  more  defined  wave 
pattern.  This  interpretation  will  later  be  complemented  where  the  14.2  cm  antenna  is 
shown  to  exhibit  the  highest  plasma  density. 
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Figure  5.17.  Three  dimensional  helicon  fields.  Antenna  lengths:  22.6,  14.24,  and  11.19  cm. 
Antenna  edge  located  at  z  =  0  cm.  Wave  propagation  and  static  magnetic  field  in  the  +z  direction. 


Figure  5.18.  Two  dimension  wave  fields  taken  across  center  of  tube.  Antenna  lengths  of  22.6, 
14.2,  and  11.2  cm.  Antenna  location  is  at  z  =  0  cm  and  wave  propagation  and  static  magnetic 
field  are  in  the  +z  direction. 


From  inspection  of  Figure  5.17  and  Figure  5.18,  the  helicon  wavelengths  can  be 
determined.  These  values  in  comparison  to  the  antenna  length  and  expected  wavelength 
are  given  in  Table  5.1.  Since  each  antenna  is  a  V2  turn  helical  (m  =  +1)  antenna,  the 
expected  helicon  wavelengths  are  defined  by  lu  =  24. 
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Table  5.1.  Summary  of  Expected  and  Measured  Wavelengths. 


3-d  Antenna 

Expected  Wavelength 

Measured  Axial  Wavelength 

Length  [cm] 

/i^i  =  2/'  [cm] 

(z-direction  only)  [cm] 

22.6 

45.2 

~25 

14.2 

28.4 

~15 

11.2 

22.4 

~19 

In  contrast,  the  expected  wavelength  may  also  be  calculated  from  the  density  profiles  at 
different  downstream  locations.  As  an  example  of  this,  the  radial  density  profiles  from 
Figure  5.4  are  matched  with  a  theoretical  Gaussian  profile  (Figure  5.19)  and  the 
wavenumbers  and  wavelengths  for  each  case  are  numerically  solved. 


I - 1 - O.OOE+00 — - 1 - 1 - 1 

-2.5  -1.5  -0.5  0.5  1.5  2.5 

Radius  [cm] 

Figure  5.19.  Radial  density  profiles  matched  with  Gaussian  to  be  numerically  solved  for 
the  wavelength.  The  full  antenna  length  is  22.6  cm. 

The  resulting  wavelengths  for  each  profile  are  shown  in  Figure  5.20.  As  would  be 
expected,  lower  densities  and  broader  Gaussian  profiles  yield  longer  anticipated 
wavelengths.  However,  the  expected  wavelengths  are  still  longer  than  the  measured  ones 
reported  above.  Recall,  the  wavelengths  reported  in  Table  5.1  and  from  inspection  of 
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Figure  5.17  and  Figure  5.18  are  based  solely  on  the  +z  directed  values  and  do  not  yet 
account  for  the  helical  nature  (r  -  component). 


Figure  5.20.  Wavelength  solution  based  on  axial  varying  Gaussian  density  profiles. 


5.2.3  Three  Dimensional  Wavelengths 

The  14.2  cm  antenna  is  used  as  a  template  for  the  following  analysis.  From  the  contour 
plots  of  Figure  5.17,  the  results  indicate  the  waves  have  a  finite  ‘thickness’  or  radius  to  its 
structure.  This  is  illustrated  more  accurately  in  Figure  5.21  for  the  14.2  cm  antenna  as 
the  wave  decreases  in  radial  structure  from  ~  1  cm  at  the  axial  z  =  0  cm  location  to  ~  0.75 
cm  one  ‘wavelength’  ‘downstream’  at  z  =  15  cm. 
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Figure  5.21.  2-d  cross  sections  showing  the  radial  decrease  in  the  wave  over  an  axial  length. 
Radial  length  decreases  from  ~  1  cm  to  ~  0.75  cm  over  ~  15  cm  axially. 

The  full  wavelength  can  then  be  calculated  by  using  a  linear  decrease  in  the  radial 

structure  (which  is  valid  from  the  slow  damping  observed  from  the  2-d  profiles  of  Figure 

5.10  -  Figure  5.12)  and  taking  the  average  radius  over  1  wavelength  to  obtain 


(5.7) 


=  15.98  cm 


'Helicon 


This  differs  from  the  2-d  wavelengths  measured  in  Figure  5.11  (depending  on  the  spatial 
location  of  measurement)  by  anywhere  from  16  %  to  37  %.  Additionally,  it  differs  from 
the  expected  value  of  28.48  cm  by  78  %  based  off  of  the  antenna  length  expectation. 
Similar  results  are  obtained  for  the  other  two  antennae  tested  where  the  full  three 
dimensional  wavelengths  are  ~  25.1  ±  0.5  cm,  16.0  ±  0.5  cm,  and  19.08  ±  0.5  cm  in 
reference  to  the  22.6  cm,  14.2  cm,  and  1 1.2  cm  antennae  respectively. 
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Although  the  full  measured  wavelengths  are  shorter  than  those  predicted  by  either  the  3-d 
antenna  length  or  calculated  from  radial  density  profiles,  the  results  demonstrate  the 
helicon  wavelength  is  most  accurately  represented  and  measured  by  3-d  imaging.  The 
remaining  challenge  is  to  resolve  the  difference  between  the  theoretically  predicted 
values  and  those  measured  in  the  lab  by  examining  the  causes  of  this  discrepancy. 

The  method  explained  in  Chapter  2,  where  a  unique  value  of  kz  was  calculated  which 
satisfied  the  boundary  condition  and  density  profile,  employs  two  main  assumptions.  The 
first  is  that  the  plasma  density  profile  extends  all  the  way  to  the  walls  and  the  second  is 
the  radial  profile  is  constant  in  the  axial  direction.  We  have  already  shown  the  radial 
density  profiles  in  the  axial  direction  are  not  constant,  and  as  a  result  the  axial  wavelength 
will  not  be  either.  However,  considering  the  22.6  cm  antenna  length  where  the  measured 
density  profde  atz  =  0  predicts  a  33.78  cm  wavelength  (Figure  5.20)  and  the  measured 
wavelength  from  3-d  imaging  is  25.1  cm,  the  difference  is  35  %.  This  is  in  contrast  to  the 
2-d  measured  wavelengths  from  Figure  5.10  which  differ  from  the  predicted  result 
anywhere  from  54  %  to  118%,  depending  on  the  chosen  2-d  scan. 

5,2,4  b„  b e Measurements 

In  addition  to  the  bz  field  measurements,  br  and  be  measurements  were  taken  in  the  y  =  0 
plane.  By  taking  the  scans  in  the  y  =  0  plane,  we  can  redefine  the  off-axis  plasma  wave 
fields  br  and  be  as  bx  and  by,  respectively.  Here,  br  or  bx  correspond  to  field  lines  which 
point  radially  inward  whereas  be  or  by  correspond  to  field  lines  perpendicular  to  ‘r’  and 
point  in  the  azimuthal  direction.  These  measurements  encompass  radial  scans  from  -a  to 
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+a  in  the  x  -  y  plane  similar  to  the  bz  seans  of  Figure  5.16.  The  results  for  the  14.2  em 
antenna  are  shown  in  Figure  5.22.  The  bz  sean  is  shown  again  to  referenee  the  signal 
magnitude  against  the  other  two  direetions.  For  this  antenna,  the  magnitude  of  the  bz 
signal  amplitude  is  signifieantly  higher  than  either  the  hr  or  be  signals  (up  to  an  order  of 
magnitude  larger  in  some  loeations).  Some  of  these  eontour  plots  will  be  revisited  later 
when  the  data  are  eorrelated  with  axial  density  profiles.  For  now,  the  wave  profiles  on 
eaeh  eross  seetion  mateh  those  predieted;  be  is  symmetrie  about  r  =  0  while  hr  is  peaked 
on  axis  and  radially  deeays.  The  unantieipated  result  is  the  axial  variation  of  the  hr 
profile,  hr  is  not  at  the  maximum  value  under  the  antenna  (where  one  might  expeet)  but 
rather  peaks  within  one  wavelength  downstream  from  the  antenna. 
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Figure  5.22.  br,  be,  b^  wave  fields  for  14.24  cm  antenna.  All  scans  were  taken  in  the  y=0 
plane  and  therefore  br  and  be  can  also  be  designated  as  b^  and  by,  respectively. 


This  unexpected  feature  was  investigated  for  the  1 1 .2  cm  antenna  length  and  similarly 
compared.  The  two  hr  comparisons  are  shown  in  Figure  5.23  with  a  noticed  similarity; 
the  peak  signal  amplitude  appears  downstream  of  the  antenna  edge.  It  is  this 
unanticipated  feature  that  inspired  a  more  thorough  analysis  of  the  axial  density  profiles. 
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Figure  5.23.  Comparison  of  br  contour  plots  for  14.2  cm  and  11.2  cm  antenna  lengths. 


5.3  Axial  Density  Measurements 

Axial  density  scans  were  conducted  for  all  three  of  the  previously  discussed  antennae; 
22.6  cm,  14.2  cm,  and  1 1 .2  cm.  The  test  conditions  were  500  Watt  input  power  and  900 
Gauss  axial  magnetic  field  with  the  antennae  location  given  in  Figure  5.9.  Density 
profiles  were  constructed  along  the  centerline  of  a  cylindrical  plasma  column  as  shown  in 
Figure  5.24  thru  Figure  5.26. 
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Figure  5.24.  Axial  density  profile  down  tube  centerline  for  22.6  cm  antenn: 
Dashed  lines  indicate  end  rings  of  antenna. 
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Figure  5.25.  Axial  density  profile  down  tube  centerline  for  14.24  cm  antenna. 
Dashed  lines  indicate  end  rings  of  antenna. 
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Figure  5.26.  Axial  density  profile  down  tube  centerline  for  11.19  cm  antenna. 
Dashed  lines  indicate  end  rings  of  antenna. 


The  discussion  on  the  structure  of  each  profile  begins  with  returning  to  the  uniform 
density  dispersion  relation  as  given  in  Chapter  2  as 


coo/ 

_ pe 

COS  6 


(5.8) 


for  cylindrical  coordinates,  and  m  =  +1  antenna  with  »  k^^ ,  this  becomes 


n  3.83  k. 


a 


(5.9) 


Therefore,  for  a  fixed  tube  radius  ‘a’,  frequency  ‘  nr  =  2;r/  ’,  the  relation  n/ Bq  should  be 
a  constant,  or  more  specifically,  proportional  to  a  constant  kz,  if  kz  is  a  fixed  value. 
Previously,  we  had  discovered  the  radial  density  profiles  (modeled  as  Gaussians)  were 
not  uniform  as  the  plasma  length  was  traversed  away  from  the  antenna.  This  indicated 
the  axial  wavenumber  was  not  constant  and  as  a  result,  n/ Bo  will  not  be  either. 

However,  for  each  antenna,  n  /Bo  is  constant  in  the  region  where  the  magnetic  field  is 
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decaying;  agreeing  with  the  description  provided  for  by  a  uniform  radial  density  profile. 
This  result  is  in  agreement  with  the  radial  profiles  taken  at  distances  away  from  the 
antenna  as  previously  shown  in  Figure  5.4. 

5.3.1  Downstream  Density  Peak 

The  structure  of  the  axial  density  profiles  not  yet  explained  is  the  region  where  the 
density  rises  downstream  of  the  source  antenna;  as  observed  in  Figure  5.25  and  Figure 
5.26  for  the  14.2  cm  and  1 1.2  cm  antennae,  respectively.  Downstream  density  peaks 
have  been  previously  identified  by  two  different  research  groups  [14,  15].  Additionally, 
both  groups  have  indentified  the  majority  of  RF  energy  absorbed  by  the  plasma  is  under 
the  antenna.  In  fact,  this  has  been  estimated  to  be  >  70-80%  of  the  total  energy  absorbed 
within  one  wavelength  of  the  antenna  and  heating  in  the  downstream  region  is  nearly 
negligible  [15]. 

In  the  next  section,  we  provide  an  explanation  for  the  rise  in  the  downstream  density  in 
terms  of  a  simple  energy  balance  within  a  plasma  volume  where  a  temperature  decay 
Te(z)  occurs  in  the  +z  direction.  Initially,  no  assumption  of  helicon  plasma  or  radial 
magnetic  confinement  will  be  considered.  After  performing  an  initial  analysis,  the  model 
will  be  applied  to  the  experimental  electron  temperature  profiles  for  all  three  antennae 
with  specific  emphasis  on  those  which  exhibit  a  downstream  density  peak. 

5. 3. 1.1  Global  Energy  Balance  Model 

The  following  analysis  refers  to  the  plasma  volume  given  in  Figure  5.27. 
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Figure  5.27.  Representation  of  plasma  volume  for  energy  balance. 


Performing  a  simple  energy  balance  for  determination  of  the  electron  density  [16] 


Pahs  enU^AajjrSj. 


(5.10) 


Physically,  this  expression  equates  the  total  power  absorbed  by  the  plasma  Pats  to  the 
total  power  lost.  The  power  lost  is  equivalent  to  an  electron-ion  pair  escaping  with 
velocity  ub  through  an  effective  area  Agff  for  an  energy  loss  Stot-  In  this  case,  the  energy 
loss  can  either  be  due  to  a  collisional  process  or  a  sheath  loss.  The  total  energy  lost  is 
then  defined  as 


^T=^c+^e+^i 


(5.11) 


where  Stot  is  the  total  energy  lost  per  electron-ion  pair.  This  is  the  result  of  energy  lost 
through  Sc  (collisional  processes;  elastic  scattering,  excitation,  ionization),  or  through 
electron  and  ion  kinetic  energy  lost  across  sheaths  near  a  surface  or  boundary,  St  and  Se 
respectively,  which  are  defined  by  [16] 

^.=27;  (5.12) 


(5.13) 
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where  Fs  is  the  sheath  voltage  for  an  insulating  wall  defined  as 


K  =  T„  In 


(5.14) 


substituting  (5.14)  into  (5.13)  beeomes 


{  M  Y  1 
- -  +-Y 

2;rm„  2 


(5.15) 


As  the  gas  in  this  work  is  Argon,  (5.15)  beeomes 


=5.187: 


(5.16) 


The  remaining  energy  loss  quantity  8c  is  the  collisional  energy  loss  per  eleetron-ion  pair 
created.  This  considers  electrons  that  can  elastically  collide,  excite,  or  ionize  the 
background  gas.  Therefore,  the  collisional  energy  for  an  electron-ion  pair  is  defined  as 


(5.17) 


For  Argon,  8ex  =  12.14  Volts  and  8iz=  15.76  Volts.  The  rate  constants  K^x  are  given  in 
Lieberman  [16]  and  approximated  by 


V,  =2.336  xlO^'^Tl'-^^e 


14  1 .609  0.06 1 8(ln  T;  )  -0. 1 1 7 l(ln rj 


(5.18) 


K.  =2.34xl0^'^r°-^®e^‘’-^^'^' 


(5.19) 


K  =2.48xl0^‘^r“^e^'"-^*'^' 


(5.20) 


where  the  rate  constants  Kxx  (m 7s)  are  valid  for  temperatures  in  the  range  -'1-7  eV  [16]. 
Additionally,  the  Bohm  velocity  ub  is  defined  as 


(5.21) 
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Therefore,  the  final  undefined  quantity  in  expression  (5.10)  is  the  effeetive  area  for 
partieles  loss  Aeff.  This  ean  be  done  for  a  cylinder  with  the  radius  ‘a  ’  and  length  ‘T  ’  as 

A^^=2S.A.,,,h,+S.A.^X  (5.22) 

Here,  S.A.  are  the  areas  for  the  respective  cylindrical  surfaces.  The  “2”  in  (5.22) 
references  both  ends  of  the  cylinder.  Therefore,  (5.22)  may  also  be  written  as 


A^jj-  =  2na  +  2n:aLh^ 


(5.23) 


hi  and  ha  are  the  sheath  density  ratios,  which  for  an  intermediate  pressure  where  the  mean 
free  path  is  less  than  the  device  dimensions  X.  <K  (a,T)  can  be  approximated  by  [16] 


h,  =:^«0.86 
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An  intermediate  pressure  is  determined  based  on  a  mean  free  path  for  Argon  given  by 


x.= 


1 


1 


n  (j.  330j!? 


cm  {p  in  Torr) 


(5.26) 


for  10  mTorr  (5.26)  becomes 


X  «  0.3  cm 


(5.27) 


and  knowing  ‘a’  and  ‘T’  are  much  larger  than  ~  0.3  cm,  the  intermediate  pressure  sheath 
density  ratios  will  be  valid. 


Returning  to  the  energy  balance  (5.10)  and  rearranging  for  the  density. 


n  = 


abs 


euB{T^)A,ff{a,L)s^{T^) 


=  n{T^,a,L) 


(5.28) 
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and  for  a  fixed  ‘a’  and  ‘L\  the  density  beeomes  solely  a  function  of  the  electron 
temperature  profile;  albeit  a  very  complex  function  of  electron  temperature  as  illustrated 
by  substitution  for  quantities  with  a  temperature  dependence 


n  =  - 


^  abs 


feT^ 
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5. 3. 1.2  Sample  Temperature  Profiles 

At  this  point,  it  is  instructive  to  assume  a  temperature  profile  in  order  to  observe  the 
behavior  of  the  density  in  the  axial  direction.  To  conform  to  the  experiment  the 
following  quantities  are  specified:  cylindrical  tube  radius  a  =  2.5  cm  ,  length  7.  =  0.5  m, 
neutral  Argon  gas  pressure  =10  mTorr,  with  500  Watts  absorbed  power.  Additionally, 
an  exponentially  decaying  temperature  profile  with  initial  temperature  =6  eVis 
defined  as 

7;(z)  =  7;oA“(^^^“)  (5.30) 

We  further  let  zq  =  0  and  the  decay  constant  a  =  3  so  that  the  axial  temperature  profile  is 
illustrated  in  Figure  5.28. 
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Figure  5.28.  Exponentially  decaying  axial  temperature  profile. 


Using  the  temperature  profile  in  (5.30)  the  axial  density  variation  can  be  solved  at  each 
“z”  location  and  is  shown  in  Figure  5.29. 
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Figure  5.29.  Axial  density  variation  based  on  an  exponentially  decaying 
temperature  profile  illustrating  a  downstream  peak  in  the  plasma  density. 
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The  downstream  density  peak  is  slightly  unexpected  based  on  the  fact  that  no  assumption 
was  made  regarding  an  increased  ionization  rate  or  other  effect  that  takes  place  outside 
the  source  region;  z  =  0.  To  better  understand  why  this  occurs,  the  quantities  that  vary 
with  Te  or  ‘z’  in  (5.28)  are  more  closely  examined. 
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Figure  5.30.  Rate  constants  Kxji  variation  due  to  temperature  profile. 


With  regards  to  Figure  5.30  and  Figure  5.31,  as  the  temperature  deereases,  the  collisional 

rates,  Kei,  Kex,  and  Kiz  all  decrease  thereby  increasing  the  total  energy  loss  downstream. 

The  decaying  or  cooling  electron  temperature  has  the  effect  of  decreasing  the  ion  velocity 
which  causes  the  density  to  rise.  This  effect  is  countered  by  an  increase  in  energy  loss 
from  collisional  and  sheath  processes  forcing  the  density  to  decrease.  However,  the  rate 
of  decay  of  the  electron  temperature  will  initially  cause  the  downstream  density  to  rise 
prior  to  decaying  away  from  the  source. 

Therefore,  dependent  upon  how  the  electron  temperature  decay  is  modeled  as  well  as 
some  specified  initial  conditions,  one  can  obtain  a  reasonable  idea  of  how  the  plasma 
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density  is  going  to  behave.  Furthermore,  the  analysis  to  this  point  has  considered  a 
simple  energy  balance  and  has  not  included  any  effects  such  as  an  external  magnetic  field 
or  wave  heating  to  explain  the  downstream  phenomenon. 


5.3. 1.3  Application  to  the  Present  Data 

Returning  to  the  11.2  cm  antenna,  the  measured  density  and  temperature  profiles  are 

shown  in  Figure  5.32. 
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Figure  5.32.  Measured  density  and  electron  temperature  profiles  for  11.2  cm  antenna. 

The  rise  in  temperature  beyond  about  30  -  35  cm  downstream  is  presumed  to  be  due  to 
insufficient  probe  compensation  at  the  lower  collected  electron  currents;  a  smaller 
collected  electron  current  will  cause  a  rise  in  the  sheath  impedance  according  to  equation 
(4.12)  where  the  probe  impedance  remains  fixed.  At  these  ‘far’  downstream  conditions, 
a  larger  probe  tip  would  have  had  the  effect  of  collecting  a  larger  electron  current 
reducing  the  error  due  to  impedance  ratios.  A  larger  compensation  probe  at  these 
conditions  would  also  have  ‘tracked’  the  plasma  potential  to  a  higher  degree  yielding 
more  precise  electron  temperature  measurements.  The  electron  temperature  profile  is 
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approximated  as  shown  in  Figure  5.33  with  corresponding  error  bars  on  every  other 
measured  data  point.  As  previously  discussed,  the  larger  error  in  the  electron  temperature 
for  the  downstream  locations  is  due  to  a  decrease  in  the  collected  saturation  current  which 
will  force  the  sheath  impedance  to  rise  with  respect  to  the  fixed  probe  impedance. 


Figure  5.33.  Modeled  and  measured  electron  temperature  profiles. 

Using  the  modeled  electron  temperature  and  following  the  same  procedure  for  calculating 
the  density  as  outlined  above,  the  measured  and  expected  density  can  be  compared  for  a 
plasma  length  of  40  cm.  This  is  shown  in  Figure  5.34. 
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Figure  5.34.  Modeled  and  measured  density  profiles  for  11.2  cm  antenna. 
The  top  graph  is  shown  on  the  same  scale  while  the  bottom  scales  are  separated. 


A  discrepancy  between  the  order  of  magnitudes  of  the  calculated  and  measured  plasma 
densities  is  immediately  apparent  from  Figure  5.34.  However,  also  obvious  is  the 
remarkable  agreement  between  the  structures  of  the  two  profiles.  Therefore,  what 
remains  is  to  account  for  the  order  of  magnitude  difference  between  the  two  profiles.  To 
do  this  we  will  have  to  make  some  modifications  to  the  plasma  volume  described  in 
Figure  5.27.  First,  consider  particle  loss  only  across  surfaces  1,  2,  and  3  in  the 
downstream  region.  This  has  the  effect  of  altering  equation  (5.22)  and  (5.23)  so  the  first 
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term  accounts  for  only  the  one  S.A.disk  cross-section;  only  surface  2  should  be  included 
for  this  area  since  small  losses  would  be  expected  in  the  direction  of  surface  4  from  the 
magnetic  field  pointing  in  the  +z  direction.  Equation  (5.22)  becomes 

A,,,=S.A.,,,h,+S.A.^A  (5.31) 

or 

A^^j  =  Tta^h^  +27taLh^  (5.32) 

This  alteration  has  a  negligible  effect  on  the  total  Aeff  because  L  »  a  in  (5.32).  However, 
returning  to  (5.28)  for  the  density,  the  effective  area  for  particle  loss  is  going  to  remain 
the  term  which  has  the  greatest  impact  on  the  order  of  magnitude  for  the  resultant  density. 
Therefore,  we  take  into  consideration  the  axially  applied  Bq  field  which  radially  confines 
the  plasma.  We  can  introduce  a  “loss  factor”  into  (5.31)  such  that 

^eff  ~  XdISK  (5'.^  'disk  )  ZcYL  {S.A.^,1.,)  (5.33) 

where  Xdisk  and  Xcyl  are  identified  as  percent  confinement  quantities,  i.e.,  0<  ;^<1.  For 
the  current  data  we  assume  no  change  to  the  axial  losses  so  that  Zdisk  =  1  -00  ,  but  for  the 
relatively  strong  magnetic  field.  Bo  =  900  Gauss,  we  set  Zcyl  =0.10  stating  that  90%  of 
the  original  plasma  is  now  confined  by  the  magnetic  field.  This  is  also  reiterated  as: 

The  magnetic  field  alters  the  effective  area  for  particles  loss  in  the  radial  direction 
(across  the  surface  area  of  a  cylindrical  boundary)  so  the  effective  area  is  reduced  to 
1/1  (f^  of  the  original  magnitude. 

Utilizing  these  definitions  alters  the  predicted  density  profiles  such  that  close  agreement 
with  the  measured  density  is  achieved.  This  is  illustrated  in  Figure  5.35. 
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A  similar  analysis  is  performed  for  the  14.2  em  antenna  where  a  downstream  density 
peak  was  also  observed.  The  measured  density  and  temperature  profiles  for  this 
condition  are  shown  in  Figure  5.36  where  the  rise  in  temperature  downstream  is  again 
presumed  to  be  the  result  of  insufficient  probe  compensation  (due  to  the  saturation 
current  as  previously  discussed)  for  the  densities  and  currents  measured. 
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Figure  5.36.  Measured  deusity  aud  temperature  profiles  for  14.2  cm  auteuna. 
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The  modeled  temperature  profile  is  again  assumed  to  be  an  exponential  decay.  The 
initial  temperature  and  decay  constants  were  selected  to  provide  the  best  fit  to  the 
measured  density  and  temperature  profiles.  The  plasma  length  was  40  cm  and  the 
cylindrical  plasma  radius  was  2.5  cm.  These  results  are  shown  in  Figure  5.37  (with  error 
bars  for  selected  temperature  data  points)  and  Figure  5.38.  In  this  case,  the  axial  loss 
factor,  Xdisk  -  ^  -00  was  not  adjusted,  however,  the  radial  loss  factor  was  adjusted  to 
XcYL  ~  0-06  in  order  for  the  density  magnitudes  to  reach  agreement. 
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The  final  axial  density  profile  analysis  was  performed  for  the  22.6  cm  antenna.  Although 
there  was  not  a  downstream  density  peak  observed,  the  model  still  accurately  predicted 


the  density  based  off  an  exponential  decaying  temperature  profile  and  a  radial 
confinement  factor  =  0.10  .  The  results  are  shown  in  Figure  5.39  and  Figure  5.40. 
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Figure  5.39.  Measured  aud  calculated  axial  deusity  profiles  for  the  22.6  cm  auteoua. 
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Figure  5.40.  Modeled  axial  electron  temperature;  axial  decay. 

The  power  of  this  simple  model  lies  in  its  global  nature.  It  is  able  to  provide  an 
explanation  for  a  downstream  density  peak  without  assuming  any  applied  magnetic  fields 
or  wave  particle  interactions.  The  model  has  been  adapted  to  account  for  radial 
confinement  which  exists  when  a  magnetic  field  is  applied  by  supplying  a  few  correction 
factors.  However,  it  was  necessary  to  change  the  correction  factor  for  confinement 
between  the  different  antennae  even  though  the  applied  magnetic  field  did  not  change. 
This  suggests  that  while  the  model  is  valid  for  providing  an  axial  density  structure  based 
on  a  cooling  electron  temperature,  the  correction  factors  may  have  to  be  further  refined, 
particularly  the  radial  loss  factor  .  Since  this  factor  represents  the  fraction  of  the 
surface  area  of  the  cylindrical  tube  for  radial  losses,  the  variation  of  this  term  suggests  a 
dependence  on  a  cross  field  diffusion  coefficient.  Should  a  more  detailed  model  of  an 
axial  density  profile  be  required,  this  component  of  the  model  may  require  further 
investigation.  The  objective  of  this  model  was  to  explain  the  formation  of  a  downstream 
density  peak  and  consequently  further  analysis  of  this  point  is  beyond  the  scope  of  the 
current  work. 
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5.4  Additional  Cases  of  Interest:  Conical  Antennae 

Two  additional  cases  of  interest  which  were  tested  in  this  work  were  for  conical  tube 
geometry.  Three  dimensional  antenna  lengths  of  22.6  cm  were  used  for  both  cases  and 
their  results  will  briefly  be  presented. 

5.4.1  Converging  Conical  Antenna 

A  converging  conical  tube  operating  at  500  Watts  in  a  peak  900  Gauss  statie  magnetic 
field  is  studied.  The  result  of  three  dimensional  imaging  of  the  propagating  wave  is 
shown  in  Figure  5.41.  The  edge  of  the  antenna  is  located  at  z  =  0  and  the  magnetic  field 
is  in  the  +z  direction. 


Y 


Figure  5.41.  Converging  conical  axial  bz  scans. 


The  radial  scans  for  each  cross  section  span  ±1.0  cm.  This  differs  from  the  eylindrieal 
tubes  where  seans  were  complete  from  ±  2.4  em.  The  radial  scans  were  limited  due  to 
probe  access  at  the  narrow  end  of  the  conical  tube;  however,  a  detailed  representation  of 
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the  wave  fields  is  given  by  Figure  5.41.  The  bz  seans  indicate  the  wave  terminates 
approximately  20  cm  downstream  of  the  antenna.  This  is  believed  to  be  a  consequence  of 
the  plasma  interaction  with  the  wall  as  propagation  is  toward  more  narrow  radii  whereby 
the  properties  of  the  plasma  dielectric  no  longer  meet  conditions  necessary  for  helicon 
wave  propagation.  This  explanation  is  complemented  by  the  axial  density  profile 
illustrated  in  Figure  5.42.  The  peak  density  values  obtained  are  less  than  any  of  the 


cylindrical  antennae  tested  and  the  axial  density  drop  occurs  in  a  region  closer  to  the 
antenna  as  a  result  of  the  narrowing  tube  radius  moving  downstream.  The  scale  of  the 
probe  size  to  tube  diameter  regarding  interference  with  the  wave  may  also  play  a  role. 
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Figure  5.42.  Converging  conical  tube  axial  density  profile. 
Dashed  lines  indicate  the  end  rinos  of  the  antenna. 


5.4.2  Diverging  Conical  Antenna 

The  second  case  of  interest  was  a  diverging  conical  tube.  The  test  conditions  were  the 
same  at  500  Watts  900  Gauss  and  the  three  dimensional  antenna  length  was  22.6  cm.  It 
is  similar  to  the  converging  tube  and  since  both  ends  of  the  tube  were  attached  to  conflat 
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flanges,  testing  this  condition  simply  required  “flipping”  the  quartz  tube  around.  Three 
dimensional  imaging  of  the  wave  via  bz  scans  are  shown  in  Figure  5.43. 


The  results  of  the  bz  imaging  show  results  familiar  to  the  cylindrical  tube  of  the  same 
antenna  length,  22.6  cm  (Figure  5.17).  The  measured  wavelength  for  the  diverging  tube 
was  25.8  ±  0.5  cm;  slightly  longer  than  the  25.1  cm  cylindrical  counterpart.  The 
difference  lies  in  the  radial  width  component  as  describe  previously. 

Axial  density  measurements  were  taken  and  compared  against  those  of  the  cylindrical 
tube  of  same  antenna  length;  Figure  5.44.  While  the  density  measurements  for  both 
versions  of  the  conical  tube  exhibit  similar  peak  densities,  both  are  still  much  less  than 
the  density  values  obtained  with  a  traditional  cylindrical  quartz  tube.  It  is  for  this  reason 
that  further  investigation  of  the  conical  geometry  was  abandoned.  For  the  same  input 
power,  the  cylindrical  tube  exhibited  better  performance  than  its  conical  counterparts. 
Performance  in  terms  of  the  peak  density  suggests  the  cylindrical  antenna  has  a  higher 
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coupling  efficiency  than  the  conical  geometry;  though  this  measurement  was  not 
quantified. 
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Figure  5.44.  Diverging  conical  tube  axial  density  profile. 
Dashed  lines  indicate  the  end  rings  of  the  antenna. 


5.5  Ionization  Cost  Analysis 

A  similar  analysis  with  regard  to  the  power  balance  performed  in  the  previous  section  can 
be  performed  to  examine  the  ionization  “cosf  ’  in  helicon  plasma.  Ionization  cost  refers 
to  the  total  energy  (or  power)  required  to  produce  an  electron-ion  pair.  We  can  define  the 
ionization  cost  in  terms  of  equating  the  power  absorbed  by  the  plasma  to  the  power  lost 
similar  to  (5.10)  as 

p  (5  34) 

at 

where  the  Pats  =  500  Watts  as  before  and  Stot  is  the  same  total  energy  required  to 
maintain  the  plasma  discharge  due  to  collisional  and  sheath  losses.  dN/dt  is  the  particle 
loss  rate  defined  as 
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(5.35) 


dN 

dt 


=  Affnc, 


Substituting  (5.35)  into  (5.34)  and  solving  for  Stot,  now  referred  to  as  the  ionization  cost, 
yields  the  following  expression 


^TOT 


^  abs 


Aff^c, 


(5.36) 


This  quantity  has  already  been  calculated,  from  equation  (5.11),  for  each  antenna  when 
analyzing  the  axial  density  profiles  and  the  result  is  illustrated  in  Figure  5.45. 


As  before,  the  assumption  has  been  made  that  the  plasma  absorbs  all  500  Watts  of  input 
power  and  there  are  no  other  loss  mechanisms  such  as  transmission  line  losses,  matching 
network  impedance  losses,  or  antenna  radiation  losses.  The  effective  area  Agff  again  has 
the  built-in  confinement  parameter  x  to  account  for  the  radial  particle  confinement  due  to 
the  static  magnetic  field. 


In  each  case,  the  antenna  edge  is  located  at  z  =  0,  with  +z  downstream  of  the  antenna. 
The  useful  value  of  the  ionization  cost  in  each  case  should  be  taken  at  z  =  0  since  this  is 
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the  antenna  edge.  Beeause  almost  all  the  RF  power  is  absorbed  in  the  near  field  region  of 
the  antenna  (less  than  one  helieon  wavelength),  assuming  500  Watts  of  power  is  being 
absorbed  far  downstream  is  unreasonable.  Therefore,  the  values  for  the  ionization  cost  at 
z  =  0  range  from  approximately  75  -  85  eV  for  the  helicon  discharges  investigated  in  this 
work. 

5.6  Chapter  Summary 

This  chapter  has  shown  the  close  agreement  between  measured  bz  profiles  with  magnetic 
induction  probes  and  those  numerically  solved  bz  profiles  when  radial  density  gradients 
are  considered.  The  onset  of  helicon  waves  has  been  determined  through  the  radial  bz 
measurements  over  a  wide  range  of  pressures,  power,  and  magnetic  field  strengths. 
Utilizing  the  measured  radial  density  profiles,  we  were  able  to  solve  for  the  bz  fields  as 
well  as  the  axial  wavelength  along  the  magnetic  field  line.  When  the  single  axis  z- 
direction  measurement  of  the  axial  wavelength  was  performed,  large  disagreement  with 
the  calculate  value  was  observed.  Therefore,  this  chapter  redefined  the  helicon 
wavelength  in  terms  of  the  full  wave  helix.  Three  dimensional  contour  maps  were 
constructed  for  the  first  time  to  visually  illustrate  a  helicon  wave.  When  the  3-d  helix  is 
considered,  agreement  with  predicted  axial  wavelength  values  is  achieved  to  within  35  %, 
contrasted  with  the  >  100%  differences  with  2-d  measured  wavelengths. 

The  downstream  density  peak  observed  has  been  modeled  in  terms  of  a  decaying  axial 
temperature  profile  away  from  the  antenna  /  source.  The  qualitative  structure  of  the 
model  agreed  with  axial  density  measurements  although  it  was  not  until  we  considered  a 
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radial  confinement  factor  that  close  quantitative  agreement  between  measured  and 
ealculated  densities  was  aehieved.  The  radial  confinement  factor  limited  radial  losses  to 
6-10  %  and  were  due  to  the  applied  magnetie  field. 
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CHAPTER  6:  CONCLUSIONS  AND  FUTURE  WORK 


This  work  has  consistently  completed  the  cycle:  theory  -  experiment  -  theory  -  . . . .  This 
chapter  will  follow  the  sequential  format  of  the  dissertation  and  highlight  the  primary 
results  of  three  dimensional  plasma  waves  and  density  peaks  away  from  an  ionizing 
source;  two  subjects  within  the  breadth  of  helicon  work  that  had  yet  to  be  fully  explored 
prior  to  this  research.  In  addition  to  the  innovative  diagnostic  achievements  (b-dot 
frequency  characterization  and  data  analysis),  this  work  has  methodically  compared  the 
non-uniform  radial  density  treatment  of  helicon  waves  with  experimentally  measured  bz 
results.  The  visual  representations  of  the  axial  helicon  wavelength  through  three 
dimensional  contour  mappings  illustrate  new  insights  into  helicon  physics  i.e.,  the  helical 
nature  of  the  propagating  wave  and  the  necessary  azimuthal  component  of  the 
wavelength.  Finally,  discussion  regarding  additional  areas  where  further  research  is 
warranted  will  be  discussed. 

6.1  On  Theory 

Given  that  helicon  waves  are  bounded  whistler  waves,  the  theoretical  treatment  of 
helicons  in  this  research  began  with  the  most  general  assumptions:  infinite  uniform  static 
medium.  The  next  step  was  to  include  a  boundary  while  assuming  wave  propagation  in 
radially  uniform  plasma  immersed  in  an  infinitely  long  magnetic  field.  The  resulting 
dispersion  relation  essentially  coincided  with  the  description  of  whistler  waves  in  free 
space.  The  radially  uniform  density  theory  also  provides  measurable  quantities  in  the 
laboratory;  the  br,  be,  and  bz  plasma  wave  fields.  The  measured  results  agreed  at  certain 
operating  conditions  with  the  uniform  theory  yet  still  failed  to  accurately  represent  all  test 
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conditions.  Therefore,  non-uniform  radial  density  distributions  represented  by  Gaussian 
profdes  were  introduced.  The  wave  fields  br,  be,  and  bz  were  numerically  solved  for 
subject  to  a  vanishing  br,  radial  magnetic  field  component,  on  the  boundary.  However, 
the  Gaussian  density  representations  did  not  precisely  eorrespond  with  the  experimentally 
determined  density  profiles.  While  the  agreement  between  the  density  profiles  and  wave 
fields  is  strong,  the  differences  are  not  necessarily  negligible  and  eannot  simply  be 
attributed  to  measurement  error.  Therefore,  a  suggestion  for  further  analysis  would  be  to 
inelude  a  vacuum  ‘gap’  layer  between  the  insulating  boundary  and  bulk  plasma.  The 
present  work  considered  the  plasma  density  profile  to  extend  to  the  tube  wall  while  in  a 
laboratory  discharge  a  vaeuum  gap  layer  or  sheath  will  exist.  Taking  this  into  account  in 
a  future  numerieal  solution  may  provide  even  further  agreement  with  laboratory 
quantities. 

Theoretieal  consideration  also  showed  that  inclusion  of  finite  eleetron  mass  and 
eollisionality  within  the  plasma  volume  generated  a  second  branch  to  the  dispersion 
relation  for  helieon  waves.  This  seeond  braneh  is  ealled  the  Trivelpiece-Gould  (TG-) 
braneh;  a  radially  direeted  eleetron  eyelotron  wave.  The  radially  direeted  nature  of  this 
wave  was  shown  to  stem  from  the  waves  perpendieular  wavenumber  being  much  larger 
than  the  parallel  wavenumber;  »  ky .  It  is  this  radial  TG-  wave  that  is  believed  to  be 

responsible  for  the  high  degree  of  ionization  within  helieon  plasma  from  whieh  eleetrons 
gain  energy.  The  existenee  of  an  eleetron  wave  in  the  near  field  region  of  the  souree 
antenna  is  consistent  with  reports  that  the  majority  of  RF  energy  is  absorbed  near  the 
antenna  and  not  transported  axially  downstream  [1].  It  is  also  eonsistent  with  the 
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relatively  high  (-'6-7  eV)  eleetron  temperatures  observed  near  the  antenna  region.  This 
braneh  was  additionally  shown  to  be  obscured  at  higher  magnetic  fields  (greater  than  a 
few  hundred  Gauss).  This  was  a  result  of  the  short  perpendicular  wavelength  at  higher 
magnetic  fields  (typically  less  than  1.0  cm).  Therefore,  this  is  an  extremely  difficult 
wave  to  measure  because  of  its  short  radial  nature.  To  date,  there  has  only  been  one 
experimental  measurement  with  a  J-dot  probe  which  measured  TG-  waves  in  a  helicon 
discharge  [2].  Additionally,  this  wave  is  more  readily  observed  at  lower  magnetic  fields 
(<  100  Gauss),  which  were  shown  to  increase  the  radial  wavelength  and  may  possibly 
ease  the  measurement  difficulty.  However,  at  lower  magnetic  fields,  the  peak  densities 
(-10  m'  )  are  significantly  lower  and  the  benefits  of  this  operational  regime  are 
indeterminate.  Though,  if  a  useful  application  is  found,  further  measurement  of  TG- 
waves  may  be  warranted. 

6,2  On  Diagnostics 

Identification  of  helicon  waves  was  demonstrated  through  a  few  different  methods. 
Utilization  of  the  matching  network  power  spectrum  demonstrated  wave  propagation 
when  a  sharp  drop  in  the  reflected  power  was  observed.  However,  this  effect  was  also 
explained  by  the  matching  networks  dynamic  capacitive  limitation,  i.e.,  the  variable 
range  of  the  vacuum  capacitors  within  the  matching  network  were  not  large  enough  for 
the  plasma  impedance  range  to  equate  a  50  Q  resistive  load.  A  future  project  may  desire 
to  report  the  capacitance  values  of  the  matching  network  in  a  very  detailed  fashion  over  a 
wide  power  and  magnetic  field  operating  spectrum.  This  should  readily  demonstrate  the 
validity  of  the  matching  network’s  limitations. 
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Measurements  with  magnetic  field  probes  in  a  high  frequency  plasma  environment  have 
been  shown  to  have  some  difficulty  in  implementation.  Internal  measurements  where 
spatial  characterization  of  the  plasma  parameters  was  desired  faced  multiple  challenges. 
First,  for  both  b-dot  probes  and  RF  compensated  Langmuir  probes,  the  size  of  the  probe 
with  respect  to  the  plasma  volume  demonstrated  probe  interference  near  borderline 
helicon  wave  conditions.  In  the  case  of  b-dot  probes,  the  nominal  probe  size  to  non- 
intrusively  measure  the  plasma  wave  fields  had  to  be  balanced  against  the  desire  for  high 
probe  sensitivity  which  increased  with  the  size  of  the  probe.  This  work  demonstrated  that 
the  most  time  efficient  and  accurate  method  to  characterize  the  entire  probe  system 
(including  line-lengths,  cable  types,  and  coupling  transformers)  was  to  utilize  an 
impedance  analyzer  where  the  entire  frequency  range  of  interest  was  assessed;  40  Hz  - 
110  MHz  in  this  report.  Post  processing  of  data  required  performing  a  Fourier  frequency 
analysis  to  compensate  for  spurious  signals  to  remove  unwanted  harmonics  of  the  driving 
frequency.  If  this  had  not  been  given  thorough  attention,  large  errors  in  quantifying  the 
wave  amplitudes  would  have  been  present. 

The  impedance  analyzer  used  to  characterize  the  b-dot  probes  also  served  to  rapidly 
identify  hundreds  of  impendence  curves  for  the  blocking  filters  used  in  Langmuir  probe 
compensation.  Because  the  plasma  sheath  impedance  can  be  on  the  order  of  a  few  kQ’s, 
the  blocking  filters  required  impedances  much  larger.  This  requirement  is  due  to  the  AC 
signal  being  much  larger  than  the  DC  component  of  interest.  Through  utilizing  large  AC 
impedance  at  the  fundamental  and  second  harmonic  of  the  driving  frequency,  the 
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blocking  filters  acted  as  voltage  dividers  so  that  the  AC  current  being  eolleeted  was 
minimized.  However,  the  capaeitance  from  the  eable  lengths  from  the  probe  tip  to  the 
data  aequisition  (oscilloscope  in  this  work)  caused  the  large  filter  impedance  to  drop  by 
nearly  an  order  of  magnitude.  While  still  mueh  greater  than  the  plasma  sheath 
impedanee,  shorter  line  lengths  (or  possibly  additional  filters)  in  future  work  would 
further  minimize  the  AC  eontribution.  This  effeet  was  most  evident  when  ealculating  the 
eleetron  temperature  at  lower  densities;  or  lower  DC  eurrent  colleetion. 

As  a  final  note  on  the  Langmuir  probe  analysis,  the  voltage  sweeps  performed  were  for  ± 
100  Volts  and  it  may  prove  worthy  to  sweep  to  a  slightly  higher  voltage  in  the  eleetron 
saturation  region.  This  could  not  be  reliably  performed  in  this  work  sinee  the  voltage 
supply  was  ±  1000  Volts  but  saturated  at  40  mA.  This  is  still  quite  a  large  current  in  the 
electron  saturation  region,  however,  either  due  to  inadequate  compensation  at  these  large 
voltages  or  probe  size  being  too  large  (eolleeting  more  current)  the  voltage  supply 
beeame  eurrent  limited  premature  of  +  100  Volts. 

6,3  On  Three  Dimensional  Results 

This  work  provided  the  first  report  on  three  dimensional  imaging  of  a  helieon  wave.  The 
contour  plots  over  a  full  wavelength  for  three  separate  antenna  lengths  provided  an 
insightful  visualization  to  helicon  waves.  They  introdueed  a  radial  eomponent  (that 
deereased  axially)  whieh  had  not  been  eonsidered  before.  This  aided  in  explaining  why 
previous  reports  on  helicon  wavelength  may  have  been  underestimated.  It  also  explains 
why  two  dimensional  measurements  may  disagree  with  theoretical  predictions  more  than 
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anticipated.  The  contour  plots  also  showed  agreement  with  theoretieal  profiles  when  the 


hr  and  be  profdes  were  measured.  Seanning  of  this  measurement  took  place  in  the  y  =  0 
plane  (or  vertical  centerline  of  the  tube).  To  fully  image  all  components  of  the  wave,  hr, 
be,  and  bz,  three  dimensional  scans  of  the  ‘r’  and  ‘6’  components  should  be  taken; 
although  little  additional  information  may  be  revealed  (unless  one  is  investigating  the 
radial  TG-  wave  near  the  antenna).  Consequently,  this  was  one  reason  this  was  not 
performed  in  this  work,  even  though  the  scanning  hardware  allowed  for  -x,  -y,  -z  scans. 
Radial  and  poloidal  scans  would  be  more  difficult  unless  modifieations  were  made. 
However,  full  -x  and  -y  scans  could  be  complete  and  the  converted  to  the  cylindrical 
counterparts,  -r  and  -0,  through  a  eoordinate  transformation. 

It  is  also  an  unfortunate  outcome  that  performing  the  three  dimensional  sweeps  was  so 
time  intensive.  The  only  way  this  proeess  might  be  expedited  is  to  include  additional 
probes  that  ean  acquire  data  simultaneously.  However,  this  may  introduee  further 
intrusion  to  the  magnetie  environment  and  obseure  the  results. 


6.4  On  Downstream  Density  Peak  and  Model 

Axial  measurements  of  the  plasma  density  revealed  an  unexpected  feature;  a  peak 
downstream  from  the  ionizing  souree  (antenna).  The  peak  was  first  observed  when  the 
14.2  cm  antenna  length  was  tested.  To  investigate  whether  the  peak  was  due  to  the 
external  magnetie  field,  it  was  believed  shortening  and  moving  the  antenna  (1 1.2  cm 
length)  would  cause  the  peak  to  shift  axial  loeations.  If  the  peak  had  not  shifted  axial 
locations,  then  the  eause  of  the  downstream  peak  would  have  to  be  related  to  the 
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externally  applied  magnetie  field.  Fortunately,  the  peak  shifted  axial  locations  forcing 
the  possibility  that  the  rise  in  density  was  a  consequence  of  increased  wave  absorption  by 
the  plasma.  However,  since  an  increase  in  electron  temperature  downstream  was  not 
observed  and  due  to  the  multiple  reports  estimating  wave  absorption  to  be  >  70  %  in  the 
near  field  region  of  the  antenna,  an  alternative  explanation  was  sought.  This  work 
demonstrated  a  simple  explanation  for  the  rise  in  downstream  density  that  is  caused  by  a 
decaying  electron  temperature  as  a  result  of  collisional  processes.  In  both  cases  however, 
performing  a  simple  energy  balance  where  particles  escape  a  cylindrical  volume  through 
radial  and  axial  end  losses  underestimated  the  expected  density  by  an  order  of  magnitude. 
To  this  point,  the  density  calculation  had  not  considered  any  applied  magnetic  field.  This 
prompted  the  introduction  of  particle  confinement  factors  to  account  for  the  magnetic 
field.  Because  the  magnetic  field  was  uniform  over  the  regions  where  the  density  peaks 
were  observed,  a  lumped  correction  factor  could  be  applied  uniformly  over  the  axial 
region  where  the  area  for  particle  loss  was  perpendicular  to  the  applied  field. 

The  objective  of  performing  the  energy  balance  was  to  formulate  an  explanation  for  the 
downstream  density  peak.  While  the  model  provided  this  result,  the  basis  for  the  particle 
confinement  factors  has  yet  to  be  defined.  Future  work  should  refine  this  model  to 
explain  the  confinement  factors,  which  are  needed  because  of  the  applied  magnetic  field. 
These  factors  are  most  likely  dependent  on  a  radially  diffusion  coefficient;  such  as  Bohm 
diffusion  in  a  magnetic  field. 
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Additionally,  due  to  the  diffieulty  in  quantifying  the  eleetron  temperature  in  the  axial 
direetion  (whieh  is  a  basis  for  the  density  model),  measurement  of  the  temperature  should 
be  verified  by  a  separate  diagnostie.  The  Langmuir  probe’s  eompensation  at  lower 
densities  should  also  be  enhaneed  to  provide  more  aeeurate  values  at  these  operating 
eonditions. 

6.5  On  Possible  Use  as  an  Electric  Propulsion  Device 

The  ionization  eost  analysis  was  determined  in  part  to  assess  the  viability  of  using  a 
helieon  source  in  an  electric  propulsion  capacity.  The  ~  75  eV  -  85  eV  required  to 
produce  each  electron  -  ion  pair  is  not  as  vast  an  improvement  as  would  have  been 
desired  when  compared  to  state  of  the  art  electric  propulsion  Hall  thrusters  where  the 
ionization  cost  ranges  ~  80  eV  -  90  eV.  However,  this  particular  device  was  not 
specifically  designed  to  optimize  the  ionization  cost. 

With  regards  to  utilizing  helicons  as  a  stand-alone  propulsion  device,  some  groups  have 
reported  claims  of  a  double  layer  formation  in  which  ions  are  axially  accelerated  out  of 
the  device.  However,  the  claim  has  not  been  independently  confirmed  or  shown  to  be 
inherent  to  helicon  generated  plasma.  This  could  in  part  be  due  to  double  layer  formation 
resulting  from  discontinuous  changes  in  geometry.  Since  nearly  all  helicon  work  to  date 
has  consisted  of  a  discharge  chamber  connected  to  a  larger  diameter  diffusion  chamber, 
the  formation  of  a  double  layer  in  these  setups  could  be  the  result  of  sudden  changes  in 
geometry,  which  would  cause  a  discontinuous  change  in  potential  structure. 

Additionally,  the  formation  of  a  helicon  double  layer  had  been  observed  at  relatively  low 
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pressures  (<  1  mTorr)  where  the  mean  free  path  for  eollisions  approaches  device  lengths. 
Results  presented  here  showed  no  evidence  of  quantifiable  helicon  wave  formation  at 
pressures  below  2  mTorr;  the  lower  pressures  consider  at  500  Watts  could  not  produce  a 
high  enough  density  to  allow  wave  propagation. 

Utilizing  helicon  waves  as  a  source  for  a  secondary  acceleration  stage  may  prove  more 
attractive  in  terms  of  overall  efficiency  if  the  ionization  cost  can  be  further  lowered.  This 
could  be  achieved  by  further  increasing  the  density  at  a  given  power  by  utilizing  a 
stronger  static  magnetic  field.  Also,  using  helicon  waves  solely  as  an  ionization  source 
would  remove  the  complexity  of  plasma  detachment  from  a  magnetic  field.  Regardless, 
helicons  do  offer  attractive  features  for  propulsive  devices  such  as:  high  density 
formation,  electrode-less  operation,  and  the  ability  to  operate  on  nearly  any  gaseous 
propellant. 

6, 6  Final  Summary 

This  work  has  provided  some  valuable  insight  into  the  significance  of  the  helicon 
wavelength.  Three  dimensional  imaging  of  the  waves  illustrated  the  necessary 
components  toward  reevaluating  the  full  helix  definition  of  wavelength;  as  opposed  to 
two  dimensional  measurements.  Propagating  waves  into  a  dense  (~  10  m'  )  plasma 
medium  demonstrated  the  wavelength  decreases  by  nearly  two  orders  of  magnitude  with 
respect  to  the  free  space  wavelength;  f =13.56  MHz  ^  22  meters.  The  antenna  length 

did  not  appear  to  play  a  significant  role  on  the  length  of  helicon  waves.  The  wavelength 
of  helicons  was  shown  (both  theoretically  and  experimentally)  to  be  determined  by  the 
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radial  density  profile  of  the  medium  as  well  as  by  the  eonstaney  of  this  radial  profile  in 
the  axial  direetion.  As  a  eonsideration  for  future  work,  determining  the  optimal  antenna 
length  would  be  valuable.  This  aspeet  was  demonstrated  by  neither  the  longest  (24.6  em) 
nor  the  shortest  (1 1.2  em)  length  antennae  generating  the  greatest  densities.  The  14.2  em 
antenna  ereated  the  most  dense  medium  for  helieon  waves  and  was  also  shown  to 
propagate  the  shortest  axial  helieon  wave.  This  “optimization”  effeet  is  most  likely  due 
to  a  eoupling  effieieney  between  antenna  and  plasma  and  eould  also  be  related  to  the 
geometrieal  radius  of  the  cylindrieal  diseharge  tube. 
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APPENDIX  A:  PLASMA  WA  VES  MODEL 


This  Appendix  will  follow  the  development  of  the  eold  plasma  dispersion  relation 
(CPDR)  to  obtain  the  plasma  dieleetrie  tensor.  Due  to  the  eomplex  mathematieal  nature 
of  plasma  waves,  the  best  description  of  the  results  is  given  schematically  by  a 
Clemmow-Mullaly-Allis  (CMA)  diagram.  The  derivation  here  follows  that  found  in 
Swanson  [1]  and  Stix  [2]  and  the  important  features  will  be  reproduced. 


A.1  Derivation  of  the  CPDR 

The  derivation  of  the  CPDR  begins  with  the  equation  of  motion  for  a  single  species  ‘y’  in 
an  electromagnetic  field: 

dv, 

m^^  =  q^-{E  +  v^xB)  (A.l) 


with  Maxwell’s  equations 


VxE  =  - 


dt 


VxB  =  jUq 


dE 


here,  the  expression  for  the  total  plasma  current  is  given  by 


(A.2) 

(A.3) 


(A.4) 


The  plasma  is  presumed  to  be  uniform  and  homogeneous  in  space  and  time  and  so 
solutions  will  take  the  form. 


E  =  EN 


i{k-r-a)t) 


.6  =  e‘ 


i{k-r-cot) 


(A.5) 
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Here,  Bo  is  the  applied  static  field  in  the  z-direction  and  is  assumed  .  With  the 


Fourier  transforms  according  to 


V  ik 


- >  -ICO 

dt 


the  equation  of  motion  now  becomes, 

-iconijV^j  =  qj(E^  +  x  B^) 


(A.6) 


Solving  for  the  velocity 


iq 


j 


m^{co^-col.) 

m  ( —  ro^. 

CJJ 


(coE^+iSjCo^jE^,) 


m.(co  -coj 
^  E 


{-isw.E^  +  coE^ 


(A.7) 


V' 

iq. 


rrijO) 


Here,  Sj  =  qj/  \qj\  is  the  sign  of  the  charge  for  the  species  and  cOcj  is  the  cyclotron 
frequency  given  by. 


%=■ 


(Ij  I  B, 

m. 


(A.8) 


If  rotating  coordinates  are  now  defined  as. 


E±=E,±iEy 


(A.9) 


along  with  the  definition  of  the  plasma  frequency. 


2 


(A.10) 


Then  the  current  density  can  be  written. 
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(A.  11) 


7±  ='^oS- 


CO 


PJ 


r-" 

7z  ='Ao^- 


a/. 

jle 


j 


m 


Now  define  the  eleetrical  eurrent  density  by, 

j  =  <yE 


(A.  12) 


and  the  eleetric  displacement  vector  by, 


D  =  s^Y.E  (A.  13) 

where  the  dielectric  tensor  K  is  given  by, 

K  =  1 — ^  (A.14) 

ICOSq 

Finally,  to  obtain  the  equivalent  dielectric  tensor  K,  the  plasma  current  (A.  15)  and  the 
displacement  current  (A.  16)  must  be  combined  such  that 

j  -  icos^E  =  -IcoSq  K  £■  (A.  1 7) 


where  K  is  now  given  by. 


-iD 

0^ 

0  ^ 

K  = 

iD 

1 

S 

0 

= 

A, 

0 

10 

0 

Pj 

1  0 

0 

A3. 

and  each  element  is  defined  as, 

1  ^  a/. 

2  j  Q)  -Q}^ 


(A.  18) 


(A.19) 


1  ^  S  CO  CO 

{R  L)  =  y 

2  jCo{co-colj) 

(A.20) 

K,^P  =  1-Y^ 

j  CO 

(A.21) 
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(A.22) 


ZO) 

- PI - 

i  CO{C0  +  SjC0^j) 

^  q/. 

K,-iK,=L  =  S-D  =  \-y - -  (A.23) 

y  o}(o}-Sj(o^j) 

Here,  R,  L,  S,  D,  and  P  refer  to  the  notation  that  deseribes  right,  left,  sum,  differenee,  and 
plasma  wave  terms  in  referenee  to  the  dieleetrie  tensor  elements. 

A.2  The  CMA  Diagram 

At  this  point,  it  is  of  interest  to  analyze  the  different  types  of  waves  that  ean  propagate 
within  plasma  dieleetrie.  To  do  this,  boundaries  ean  be  eonstrueted  that  represent  the 
eutoffs  and  resonanees  of  eaeh  wave  by  setting  R,  L,  and  S  equal  to  zero  and  infinity, 
resulting  in  i.e., 

R^cc  {electron  cyclotron  resonance) 
f  ^  00  {ion  cyclotron  resonance) 

5^0  {upper  &  lower  hybrid  frequencies) 

2 

CO  CO 

and  then  graphing  the  ehange  in  — ^  with  respeet  to  — ^  where  cOce  is  the  eleetron 

CO  CO 

eyelotron  frequeney,  cOp  is  the  eleetron  plasma  frequeney,  and  co  is  the  driving  frequeney. 

When  eonsidering  the  eutoffs  and  resonanees  for  a  two  speeies  plasma,  namely  ions  and 
eleetrons,  we  must  eonstruet  the  diagram  using  equations  A.17-A.21  with  the  substitution 
for  the  ion  mass  as  mi=1836  This  redefines  the  ion  plasma  and  eyelotron  frequeneies 
sueh  that 
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(A.24) 


^2  _ 

Sf^m.  \S36Sf^m^ 


and 


_  g-gp  _  _  (^ce 

m-  1836^^  i^i/ 


(A.25) 


Solving  A.17-A.21  for  the  cutoffs  and  resonances,  we  obtain  the  CMA  diagram.  This  is 


fTt  / 

shown  in  Figure  A.  1  where  the  mass  ratio  is  taken  to  be  '/  =5  to  illustrate  the  cutoffs 


and  resonances. 


Figure  A.l.  CMA  diagram  with  mi  /  me  =  5. 


At  this  point,  to  better  understand  the  influence  the  mass  ration  has,  we  can  consider  a 
few  different  gases  such  as  Fielium,  Argon,  and  Xenon  driven  at /=  13.56  MHz.  The 
CMA  diagram  can  then  be  transformed  into  a  n  V5.  Bo  graph  as  shown  in  Figure  A.2 
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where  the  location  of  the  various  ion  and  electron  resonances  are  given  in  terms  of 
magnetic  field.  It  is  no  surprise  that  the  ion  cyclotron  resonant  B-field  increases  with 
increased  mass,  however,  the  point  of  interest  for  helicon  waves  (usually  defined  between 
the  lower  hybrid  and  electron  cyclotron  frequencies)  is  that  the  'Bo-space  ’  decreases  for 
lighter  gases,  i.e.,  helicon  waves  for  Helium  should  range  ~  5  Gauss  <  ~  295  Gauss . 
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Figure  A.2.  Density  vs.  Magnetic  field  graphs  for  Helium,  Argon,  and  Xenon  driven  at  13.56  MHz. 
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APPENDIX  B:  CODE  FOR  NON-UNIFORM  DENSITY  WA  VE  FIELDS 


This  Appendix  provides  the  eode  implemented  in  MATHEMATICA  for  solving  the  non- 
uniform  radial  density  magnetie  field  profiles  deseribed  in  Chapter  2.  The  primary 
outputs  are  graphs  and  data  tables  for  the  bz  and  bz  ’  profiles  as  well  as  the  wavenumbers 
and  wavelengths  that  satisfy  the  boundary  eonditions  for  eaeh  partieular  profile. 


B. I  MATHEMATICA  Code 


Clear[c,e,freq,w,uAlA2A3A4A5A6,kl,k2,k3,k4,k5,k6,k0,a,B0,n0,e,m,yl,Y2,Y3,y4,y5,Y6,wl,w2,w3,w 

4,w5,w6,r,al,a2,a3,a4,a5,a6,/31,/32,/33,/34,/35,/36,sl,s2,s3,s4,s5,s6,fl,f2,G,f4,f5,f6,gl,g2,g3,g4,g5,g6,nl,n 

2,n3,n4,n5,n6,sll,s22,s33,s44,s55,s66] 


(*  Clears  the  variables  used  in  the  program  *) 


c=299792458; 

6=8.85418782*10^-12; 
freq=13.56*10''6; 
w=2  7T  freq; 

U=1.25663706*10^-6; 

(*  Defines  constants  *) 

kl=23.710; 

k2=23.3013155; 

k3=22.38015; 

k4=2 1.09 145; 

k5=19.60375; 

k6=16.9723; 

(*  Values  for  wave  numbers  that  satisfy  the  boundary  condition  *) 

Al=2  7T/kl 
A2=2  7T  /  k2 
A3=2  7T/k3 
A4=2  7T  /  k4 
A5=2  7T/k5 
A6=2  7T  /  k6 

(*  Calculates  the  wavelength  for  each  profile  *) 

k0=w/c; 

a=0.030; 

B0=0.09; 

e=1.60217646*10''-19; 

m=l; 


166 


(*  Some  more  constants  specific  to  the  experiment  *) 

yl=l-(kO/kl)^; 

T2=l-(k0/k2)^; 

y3=l-(k0/k3)^; 

y4=l-(k0/k4)^; 

y5=l-(k0/k5)^; 

y6=l-(k0/k6)^; 

(*  Defines  y;  term  that  accounts  for  displacement  current  *) 
n0=1.0*10''19; 

(*  Initial  peak  centerline  density  value  *) 

wl=0.005; 

w2=0.010; 

w3=0.015; 

w4=0.021; 

w5=0.030; 

w6=l; 

(*  Exponents  for  Gaussian  profiles  *) 

nl[r_]~nO  e(-(r/wl)'^2) 
n2[r_]~n0  e(-(r/w2)'^2) 
n3[r_]~n0  e(-(r/w3)'^2) 
n4[r_]~n0  e(-(r/w4)'^2) 
n5[r_]~n0  e(-(r/w5)'^2) 
n6[r_]~n0  e(-(r/w6)'^2) 

(*  Radial  density  profiles  *) 

al[r_]:=w/kl  (e  iJ  nl[r]  )/B0 
a2[r_]:=w/k2  (e  lu  n2[r]  )/B0 
a3[rj  :=w/k3  (e  ^  n3[r]  )/B0 
a4[r_]:=w/k4  (e  lu  n4[r]  )/B0 
a5[r_]:=w/k5  (e  iJ  n5[r]  )/B0 
a6[r_]:=w/k6  (e  iJ  n6[r]  )/B0 

/31[rJ:=al[r]^-klVl^ 

I32[r_]:=a2[rf-k2^  y2^ 

/33[rJ:=a3[r]^-k3V3^ 

I34[r_]:=a4[rf-k4^  y4^ 

/35[rJ:=a5[r]^-k5^  y5^ 

/36[rJ:=a6[r]^-k6^  y6^ 

(*  Defining  wavenumbers  *) 

fl[rJ-l/r-(2  al[r]  al'[r])//31[r] 
f2[rJ-l/r-(2  a2[r]  a2'[r])//32[r] 

G[rJ-l/r-(2  a3[r]  a3'[r])//33[r] 
f4[rJ-l/r-(2  a4[r]  a4'[r])//34[r] 
f5[rJ-l/r-(2  a5[r]  a5'[r])//35[r] 
f6[rJ-l/r-(2  a6[r]  a6'[r])//36[r] 
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gl[rJ:=/31[r]/Yl-mV-m/kl  al'[r]/(  yl  r)  (l+(2  kl^  yl^)//31[r]) 
g2[rJ:=/32[r]/y2-mV-m/k2  a2'[r]/(  y2  r)  (l+(2  k2^  y2^)//32[r]) 
g3[rJ:=/33[r]/y3-mV-m/k3  a3'[r]/(  y3  r)  (l+(2  k3^  y3^)//33[r]) 
g4[rj-/34[r]/y4-m^/r^-m/k4  a4'[r]/(  y4  r)  (l+(2  k4^  y4^)//34[r]) 
g5[rj :=/35[r]/y5-mV-m/k5  a5'[r]/(  y5  r)  (l+(2  k5^  y5^)//35[r]) 
g6[r_]:=/36[r]/y6-m^/r^-m/k6  a6'[r]/(  y6  r)  (l+(2  k6^  y6^)//36[r]) 

(*  Defining  f(r)  and  g(r)  from  defined  in  Chapter  2  *) 

Solve[al[rsl]==yl  kl,rsl]; 

Solve[a2[rs2]==y2  k2,rs2]; 

Solve[a3[rs3]==y3  k3,rs3]; 

Solve[a4[rs4]==y4  k4,rs4]; 

Solve[a5[rs5]==y5  k5,rs5]; 

Solve[a6[rs6]==y6  k6,rs6]; 

(*  Solves  for  location  of  singularity  *) 

sl=NDSolve[{Bzl"[r]+fl[r]Bzl'[r]+gl[r]Bzl[r]==0,mal[a]  Bzl[a]+akl  yl 
Bzl'[a]==0,Bzl'[.0005]==l},Bzl[r],{r, 0.0001, a}] 

sll=NDSolve[{Bzl"[r]+fl[r]Bzl'[r]+gl[r]Bzl[r]==0,mal[a]  Bzl[a]+akl  yl 
Bzl'[a]==0,Bzl'[.0005]==l},Bzl'[r],{r, 0.0001, a}] 

s2=NDSolve[{Bz2"[r]+f2[r]Bz2'[r]+  g2[r]Bz2[r]=:=0,m  a2[a]  Bz2[a]+a  k2  y2 
Bz2'[a]==0,Bz2'[.0001]==l},Bz2[r],{r,0.0001,a}] 

s22=NDSolve[{Bz2"[r]+f2[r]Bz2'[r]+  g2[r]Bz2[r]==0,m  a2[a]  Bz2[a]+a  k2  y2 
Bz2'[a]  ==0,Bz2'[.0001]  ==  1 }  ,Bz2'[r],  {r,0.000 1  ,a}] 

s3=NDSolve[{Bz3"[r]+G[r]Bz3'[r]+  g3[r]Bz3[r]==0,m  a3[a]  Bz3[a]+a  k3  y3 
Bz3'[a]==0,Bz3'[.0001]==l},Bz3[r],{r,0.0001,a}] 

s33=NDSolve[{Bz3"[r]+G[r]Bz3'[r]+  g3[r]Bz3[r]==0,m  a3[a]  Bz3[a]+a  k3  y3 
Bz3'[a]==0,Bz3'[.0001]==l},Bz3'[r],{r,0.0001,a}] 

s4=NDSolve[{Bz4"[r]+f4[r]Bz4'[r]+  g4[r]Bz4[r]==0,m  a4[a]  Bz4[a]+a  k4  y4 
Bz4'[a]  ==0,Bz4'[.0001]  ==  1 }  ,Bz4[r],  {r,0.000 1  ,a}] 

s44=NDSolve[{Bz4"[r]+f4[r]Bz4'[r]+  g4[r]Bz4[r]==0,m  a4[a]  Bz4[a]+a  k4  y4 
Bz4'[a]  ==0,Bz4'[.0001]  ==  1  },Bz4'[r],  {r,0.0001,a}] 

s5=NDSolve[{Bz5"[r]+f5[r]Bz5'[r]+  g5[r]Bz5[r]==0,m  a5[a]  Bz5[a]+a  k5  y5 
Bz5'[a]==0,Bz5'[.0001]==l},Bz5[r],{r,0.0001,a}] 

s55=NDSolve[{Bz5"[r]+f5[r]Bz5'[r]+  g5[r]Bz5[r]==0,m  a5[a]  Bz5[a]+a  k5  y5 
Bz5'[a]  ==0,Bz5'[.0001]  ==  1  },Bz5'[r],  {r,0.0001,a}] 

s6=NDSolve[{Bz6"[r]+f6[r]Bz6'[r]+  g6[r]Bz6[r]==0,m  a6[a]  Bz6[a]+a  k6  y6 
Bz6'[a]  ==0,Bz6'[.000 1]  ==  1 }  ,Bz6[r],  {r,0.000 1  ,a}] 

s66=NDSolve[{Bz6"[r]+f6[r]Bz6'[r]+  g6[r]Bz6[r]==0,m  a6[a]  Bz6[a]+a  k6  y6 
Bz6'[a]  ==0,Bz6'[.0001]  ==  1  },Bz6'[r],  {r,0.0001,a}] 

(*  Solves  differential  equations  for  each  profile  and  stores  the  output  array  in  “sx”  and  “sxx” 

Plot[{Bzl[r]/.sl,Bz2[r]/.s2,Bz3[r]/.s3,Bz4[r]/.s4,Bz5[r]/.s5,Bz6[r]/.s6},{r,0.0005,a},PlotRange^{{0,a},{- 
.0002,. 0015}},  AxesLabeW  {Radius,Bz}] 
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Plot[{Bzl'[r]/.sll,Bz2'[r]/.s22,Bz3'[r]/.s33,Bz4'[r]/.s44,Bz5'[r]/.s55,Bz6'[r]/.s66},{r, 0.0005, a}, PlotRange^ 
{{0, a}, {-.0002, .0015}},  AxesLabeW  {Radius,Bz'}] 


Plot[{nl[r],n2[r],n3[r],n4[r],n5[r],n6[r]},{r,0,a},AxesLabeW{Radius,Density},PlotRange^{{0,a},{0,n0}}] 
Plot} {n'[r]},{r,0, a}, AxesLabel^ (Radius,  Density  Prime}]; 

Plot}  {a}r] } ,  {r,0,a}  ,AxesLabel^  (Radius,  Alpha}  ] ; 

Plot}  {a'}r] } ,  {r,0,a} , AxesLabeW  (Radius, Alpha  Prime }  ] ; 

Plot}  {/3}r] } ,  {r,0,a}  ,AxesLabel^  {Radius,Beta}  ] ; 

Plot}  (f}r] } ,  {r,0,a}  ,AxesLabeW  {Radius,!}  ] ; 

Plot}  {g}r] } ,  {r,0,a}  ,AxesLabel^  {Radius,g}  ] ; 


(*  Plots  of  radially  varying  quantities  *) 

Export}"Bzl.esv",Table}Bzl}r]/.sl,{r,0.0005,a,.00005}]]; 

Export}"Bz2.esv",Table}Bz2}r]/.s2,{r,0.0005,a,.00005}]]; 

Export}"Bz3.esv",Table}Bz3}r]/.s3,{r,0.0005,a,.00005}]]; 

Export}"Bz4.esv",Table}Bz4}r]/.s4,{r,0.0005,a,.00005}]]; 

Export}"Bz5.esv",Table}Bz5}r]/.s5,{r,0.0005,a,.00005}]]; 

Export}"Bz6.esv",Table}Bz6}r]/.s6,{r,0.0005,a,.00005}]]; 

Export}"Bzll.esv",Table}Bzl'}r]/.sll,{r,0.0005,a,.00005}]]; 

Export}"Bz22.esv",Table}Bz2'}r]/.s22,{r,0.0005,a,.00005}]]; 

Export}"Bz33.esv",Table}Bz3'}r]/.s33,{r,0.0005,a,.00005}]]; 

Export}"Bz44.esv",Table}Bz4'}r]/.s44,{r,0.0005,a,.00005}]]; 

Export}"Bz55.csv",Table}Bz5'}r]/.s55,{r,0.0005,a,.00005}]]; 

Export}"Bz66.csv",Table}Bz6'}r]/.s66,{r,0.0005,a,.00005}]]; 

(*  Exports  bz  and  bz’  to  table  arrays  of  numerie  data  *) 
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APPENDIX  C:  IGOR  CODE  FOR  EFT  OF  B-DOT  SIGNALS 


The  following  code  was  implemented  in  IGOR  for  rapid  analysis  of  b-dot  probe  signals. 
The  outputs  are  5  graphs:  2  raw  data  signals,  2  FFT  signals  corresponding  to  the  raw  data, 
and  1  trace  corresponding  to  the  phase  delay  between  both  signals.  The  outputs  are 
shown  in  Figure  C.  1  at  the  start  of  this  Appendix  with  the  code  to  follow.  Additionally, 
the  wave  amplitudes  for  each  probe  obtained  from  the  FFT  analysis  and  the  numerical 
value  of  the  phase  between  the  two  signals  is  computed.  The  naming  convention  of  the 
two  different  probes  is:  “Probe”  and  “Reference”.  For  this  example  trace,  the  output 
values  are 


Frequency  of  Maximum:  13.33  MHz 
Probe  FFT  Value  at  the  Max:  46606.2 
Probe  FFT  Value  at  13.56  MHz:  46072.3 

Frequency  of  Maximum:  13.48  MHz 
Reference  FFT  Value  at  the  Max:  4980.33 
Reference  FFT  Value  at  13.56  MHz:  4972.04 

Time  Shift  between  Signals:  -2.9408e-08 

Probe  Signal:  3.686  V 
Reference  Signal:  0.398  V 
Phase  in  Degrees:  -143.558  degrees 
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Probe :  Raw 


Probe :  FFT 


Reference :  Raw  Reference :  FFT 


Phase  Delay 


Figure  C.l.  Output  from  IGOR  code  of  aualysis  of  b-dot  probe  siguals. 
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C.l  Code  Written  and  Implemented  in  IGOR 

#pragma  rtGlobals=l 

//  Use  modern  global  aeeess  method. 

Function  all() 

variable  /G  psig 
variable  /G  Phase 
variable  /G  rsig 

SetScale/P  x  0,1.6e-ll,"s",  X_0_ 

SetScale  d  0,0, "V",  X_0_ 

//Set  Scale  for  Probe 

SetScale/P  x  0,1.6e-ll,"s",  X_l_ 

SetScale  d  0,0,"V",  X_l_ 

//Set  Scale  For  Reference 

FFT/OUT=3/DEST=Probe_FFT  X_0_ 
FFT/OUT=3/DEST=Reference_FFTX_l_ 

//EET  of  Probe  and  Reference 

Display  /W=(l 6.5,39.5,41 1,248)  X_0_ 

Display  /W=(15.75,268.25,410.25,476.75)  X_l_ 
Display /W=(416.25,41,810.75,249.5)  Probe_PPT 
SetAxis  bottom  0,le+08 
ModifyGraph  mode=2,lsize=5 
Showinfo 

Display /W=(415. 5,268.25,810,476.75)  Reference_PPT 
SetAxis  bottom  0,  le+08 
ModifyGraph  mode=2,lsize=5 
Showinfo 

//Graph  Size  and  Eocations  and  Eormat 

//Activate  Graph  and  Place  Cursors  on  the  Graph 

Do  Window  /E  Graph2 
Showinfo  /W=Graph2 
Cursor  /P  A,Probe_PPT,  3 
Cursor  /P  B,Probe_PPT,  7 

CurvePit  /q  gauss  Probe_EFT[pcsr(A),pcsr(B)]  /D 
//Eit  Reference  EET  with  Gaussian 
wave  fit  Probe  FEE 


172 


//Statistics  for  Fit 

wavestats  /q  fit_Probe_FFT 

printf  "Frequency  of  Maximum:  %.2WlPHz\r",  V_maxloc 

//Output  Max  Value  and  Frequency 

printf  "Probe  FFT  Value  at  the  Max:  %g\r",  V_max 
eursor  /P  A  ,fit_Probe_FFT,  (  x2pnt(fit_Probe_FFT,  13.56e6)) 

//Output  Value  and  13.56  MHz 

printf  "Probe  FFT  Value  at  13.56  MHz:  %g\r",  vcsr(A) 

psig=vesr(A)*2/numpnts(X_0_) 

Do  Window  /F  Graph3 

//Activate  Graph  and  Place  Cursors  on  the  Graph 

Showinfo  /W=Graph3 
Cursor  /P  A,Reference_FFT,  3 
Cursor  /P  B,Reference_FFT,  7 

CurveFit  /q  gauss  Referenoe_FFT[pcsr(A),posr(B)]  /D 
//Fit  Reference  FFT  with  Gaussian 
wave  lit_Referenee_FFT 
//Statistics  for  Fit 

wavestats  /q  fit_Referenee_FFT 

printf  "Frequeney  of  Maximum:  %.2WlPHz\r",  V_maxloc 

//Output  Max  Value  and  Frequency 

printf  "Reference  FFT  Value  at  the  Max:  %g\r",  V  max 

eursor /P  A  ,lit_Reference_FFT,  (  x2pnt(fit_Referenoe_FFT,  13.56e6)) 

printf  "Referenee  FFT  Value  at  13.56  MHz:  %g\r",  vcsr(A  ) 

//Output  Value  and  13.56  MHz 

rsig=vcsr(  A)  *  2/ numpnts(X_  1  _) 

//Converting  FFT  Output  to  Volts 


173 


Duplicate/0  X_l_,W_Correlation;DelayUpdate 

//Using  correlation  to  find  phase  between  probe  and  referenee  signal 

Correlate  X_0_,  W_Correlation;DelayUpdate 
Display  /W=(16. 5, 497, 41 1,705.5)  W_Correlation 
wavestats  /q  W_Correlation 
print  "Time  Shift  between  Signals:",  V  maxloc 
Phase=(3  60  *  Vmaxloe/ (1/13. 5  6e6)) 

printf'V" 

printf  "Probe  Signal:  %8.3f',  psig 
//Display  Outputs  to  use  in  Exeel 
printf"  VV 

printf  "Referenee  Signal:  %8.3f',  rsig 
printf"  VV 

printf  "Phase  in  Degrees:  %8.3f'.  Phase 

printf  "  degrees  \r" 

printf'V 

Do  Window  /H 

KillWindow  Graph4 
KillWindow  Graph3 
KillWindow  Graph2 
KillWindow  Graph  1 
KillWindow  GraphO 

//Closing  Windows 

KillWaves  /A/Z;  KillVariables  /A/Z;  KillStrings  /A/Z 
//Clearing  Data 
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APPENDIX  D:  SAMPLE  IV  CHARA  CTERISTIC  AND  ANAL  YSIS 


This  Appendix  will  analyze  a  representative  IV  eharaeteristie  obtained  in  RF  helieon 
plasma  at  10  mTorr,  500  Watts,  900  Gauss,  and  13.56  MHz.  The  analysis  will  follow  the 
method  by  Ruzic  in  “Eleetrie  Probes  for  Low  Temperature  Plasmas”  [1]. 

D.I  Constructing  an  IV  Characteristic 

Langmuir  probe  analysis  of  an  IV  eharaeteristie  begins  by  measuring  the  eurrent  to  a 
probe  tip  while  sweeping  the  voltage  between  the  ion  and  eleetron  saturation  regions. 

The  ion  saturation  region  of  an  IV  eharaeteristie  oeeurs  when  the  probe  tip  is  biased  at 
sueh  a  high  negative  voltage  with  respeet  to  the  floating  potential  that  eleetrons  are 
repelled  from  the  probe  tip  and  only  ions  are  eolleeted.  The  eorollary  is  the  eleetron 
saturation  region  where  the  probe  tip  is  biased  at  sueh  a  high  positive  voltage  with  respeet 
to  floating  potential  that  only  eleetrons  are  eolleeted  and  nearly  all  ions  are  repelled.  The 
floating  potential  referred  to  here  is  the  voltage  at  whieh  an  equal  eurrent  due  to  ions  and 
eleetrons  are  eolleeted  at  the  probe  tip  and  therefore  the  total  eurrent  to  the  probe  is  zero. 
These  regions  are  shown  below  in  a  representative  IV  eharaeteristie. 
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ion  saturation  electron  electron  saturation 

retardation 


Figure  D.l.  Representative  IV  Characteristic  with  clearly  defined  regions  and  transitions. 

Typical  construction  of  an  IV  characteristic  is  accomplished  through  independent 
measurement  of  1)  the  swept  voltage  to  the  probe  tip  and  2)  the  current  drawn  by  the 
probe.  These  two  traees  are  then  recorded  on  an  oseilloscope:  voltage  vs.  time;  current 
vs.  time.  The  voltage  and  current  vs.  time  traces  are  shown  in  Figure  D.2. 
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The  raw  voltage  vs.  time  and  eurrent  vs.  time  traees  can  then  be  used  to  construct  the  IV 
characteristic.  However,  because  of  residual  RF  noise  present  in  the  traces  as  seen  as 
small  oscillations  in  Figure  D.2  and  to  make  the  25,000  data  point  sample  set  more 
manageable  for  processing,  a  1000  point  cubic  spline  was  fit  to  the  data.  Both  the  raw  IV 
trace  and  the  cubic  spline  fit  are  shown  in  Figure  D.3. 


Figure  D.3.  Raw  IV  Characteristic  with  1000  point  cubic  spline  fit. 
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As  a  matter  of  preference,  the  current  here  at  large  negative  potentials  is  positive  because 
it  is  due  to  ion  currents  while  the  current  is  negative  at  large  positive  potentials  resulting 
from  electrons.  The  differences  between  the  well  defined  characteristic  of  Figure  D.l 
and  the  measured  IV  characteristic  Figure  D.3  primarily  lay  in  the  electron  saturation  and 
electron  retardation  regions  of  the  curve;  sometimes  referred  to  as  the  transition  region. 
The  remainder  of  this  appendix  will  focus  on  obtaining  the  relevant  quantities  from  the 
IV  characteristic  in  order  to  accurately  quantify  the  plasma. 

D.2  Analysis  of  IV  Characteristic 

In  order  to  begin  the  analysis,  we  first  need  to  determine  what  regime  the  plasma  is  in, 
i.e.,  collisionless  vs.  collisional  and  thin  vs.  thick  sheath.  This  determination  is  based  on 
the  radius  of  the  measuring  probe,  the  Debye  length,  and  the  mean  free  path  for  electron- 
neutral  collisions.  For  the  majority  of  the  analysis  performed  in  this  work,  the  neutral 
pressure  was  10  mTorr.  Using  the  definition  of  the  mean  free  path  for  electron-neutral 
collisions 

^=—  (D.l) 

19  2 

where  the  collisional  cross  section  g  is  approximated  as  5  x  10"  m  ,  the  mean  free  path 
reduces  to 

— ^  (D2) 

p(mTorr) 

and  from  p  =  10  mTorr,  Xq  =  0.61  cm. 
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The  cylindrical  probe  radius  (tungsten  wire)  was  rp  =0.01443  cm  and  the  last  quantity  to 
determine  is  the  Debye  length  given  as 


^.7430  IMl 

V  e  V  ] 


(D.3) 


However,  to  get  an  estimate  for  the  Debye  length  and  determine  the  plasma  regime,  we 
require  an  estimate  for  Te  and  ng.  As  a  first  approximation  for  helicon  plasma  where 
densities  are  approximately  10  -  10  m'  and  assuming  a  temperature  of  5  eV,  the 
Debye  length  will  be  approximately,  0.0005  -  0.0016  cm  and  4Xd  0.0021  -  0.0066  cm. 
Therefore,  for  the  plasma  in  consideration  here. 


4/1^  <  <  Aq 

0.0066  cm  <  0.01443  cm  <  0.61  cm 


(D.4) 


the  regime  is  collisionless  thin  sheath  [1].  However,  as  long  as  we  know  that  the  sheath 
is  collisionless,  a  thin  sheath,  thick  sheath,  or  transitional  between  the  two  can  be  solved 
by  the  Lafromboise  method.  This  is  a  solution  to  Poisson’s  equation  based  on  the  ratio  of 
the  probe  radius  to  Debye  length.  In  this  case  the  ion  current  to  the  probe  is  given  by  [1] 


n.[m-’]  =  L6.lO"  /-[‘-'""I  . 

V  neV] 


(D.5) 


where  p  =  40  for  Argon,  Aprobe  =  2.72  x  10'^  m^,  and  J.  is  a  dimensionless  quantity 

which  accounts  for  ion  acceleration  in  the  pre-sheath  and  is  graphically  determined  in 
Figure  D.4. 
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e(V  -  V)  /  kT 

plasma  e 


Figure  D.4.  Lafromboise  plot  for  various  values  of  probe  radius  to  Debye  length. 

Plot  taken  from  Ruzic  (1994). 


However,  since  j*  is  a  function  of  Debye  length,  which  is  a  function  of  density,  the 
equation  for  the  density  becomes  transcendental  and  requires  an  iterative  solution.  In 
Figure  D.4,  the  sheath  is  considered  “thick”  for  rp  /  A-d  0  and  “thin”  for  rp  /  A,d  100. 
Additionally,  an  approximation  to  the  curves  in  terms  of  X*  for  j*  can  be  utilized.  For 
example,  in  terms  of  a  thick  sheath  approximation,  the  curve  fit  is  given  as 


and  since 


7, 


1.127(v;f 


(D.6)  and  (D.7)  are  inserted  into  (D.5)  to  give 


1.42x10'" 


I[amps]^J  ju[amu] 


A  probe  W 


(D.6) 


(D.7) 


(D.8) 
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which  is  independent  of  eleetron  temperature.  Furthermore,  while  the  thick  sheath  ease 
was  selected  as  a  starting  point  where  rp  /  A,d  0,  the  curve  could  have  been 
approximated  for  any  of  the  ratios  in  Figure  D.4  and  the  only  ehange  would  be  the 
numerieal  constant  in  (D.8).  Additionally,  because  an  iterative  method  to  solution  is 
being  utilized,  the  only  effeet  that  the  initial  constant  has  is  on  the  initial  convergence, 
i.e.,  the  initial  ratio  of  rp  /  A,d  will  converge  more  quickly  if  the  initial  estimate  for  the 
density  is  “elose”  to  the  final  solution.  Approximations  to  the  eurves  in  terms  of 
constants  which  can  be  used  in  (D.6)  ean  be  seen  in  Figure  D.5.  This  representation  is 
only  used  to  illustrate  that  an  approximation  to  remove  the  electron  temperature  in  terms 
of  X*  was  originally  given  in  the  “thiek”  sheath  ease,  though  other  approximations  could 
be  used  whieh  would  not  alter  the  final  solution  for  the  density. 


Figure  D.5.  Approximate  constant  values  for  probe  radius  to  Debye  length  ratios. 

Therefore,  aecording  to  (D.5)  or  (D.8),  the  quantities  to  calculate  the  plasma  density  have 
been  identified.  As  is  the  situation  in  most  RF  plasma  IV  oharaeteristies,  the  diffieult 
quantity  to  identify  is  the  plasma  potential.  This  also  makes  accurate  measurement  of  the 
electron  temperature  difficult. 
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D.3  Plasma  Potential,  Electron  Temperature,  and  Density 

Typically,  the  plasma  potential  ean  be  indentified  in  two  ways.  Primarily  it  is  done 
through  identifying  the  “knee”  in  the  IV  trace,  although,  as  seen  from  Figure  D.3,  this 
loeation  can  be  ambiguous.  Therefore,  the  next  step  is  to  ealeulate  the  derivative  of  the 
eurrent  with  respeet  to  voltage  and  find  the  maximum  /  minimum  (depending  on  whether 
the  eurrent  was  taken  to  be  positive  or  negative).  However,  the  derivative  method  toward 
obtaining  the  plasma  potential  is  often  obscured  in  the  transition  region  where  there  is  a 
large  amount  of  noise  as  is  often  the  ease  in  RF  plasma.  Therefore,  the  method  used  in 
this  work  is  to  first  remove  the  ion  eurrent  by  approximating  the  ion  current  in  the  ion 
saturation  region  (through  a  parabola  apex  at  the  plasma  potential  Vpiasma)  and  subtraeting 
it  from  the  total  current.  This  leaves  a  trace  for  the  electron  current  to  the  probe  as  a 
funetion  of  voltage.  Then  the  log  of  the  eleetron  eurrent  vs.  voltage  ean  be  plotted  as 
shown  in  Figure  D.6. 


Figure  D.6.  Plot  of  Ln(Ie)  vs.  Voltage  for  determine  the  plasma  potential. 

Approximating  the  linear  region  before  the  plasma  potential  and  the  linear  region  in  the 
eleetron  saturation  region,  the  plasma  potential  is  given  by  the  intersection  of  the  two 
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lines  as  shown  in  Figure  D. 6  to  be  Vpiasma  =  56.03  Volts.  The  floating  potential  (location 
where  total  current  is  zero;  equal  parts  electron  and  ion  current)  is  found  from  the  IV 
characteristic,  Figure  D.3,  as  Vfioating  =  52.87  Volts.  The  approximation  used  to  calculate 
the  plasma  potential  is  based  on  the  correct  subtraction  of  the  ion  current  from  the  total 
current  dependent  on  identifying  the  correct  plasma  potential.  The  values  used  to 
calculate  the  ion  current  can  be  adjusted,  which  will  change  the  calculated  electron 
current  and  alter  the  Ln  (le)  graph.  Figure  D.6.  Therefore,  altering  the  approximation  to 
the  ion  current  will  affect  the  value  of  the  plasma  potential  and  in  turn  the  density; 
however,  after  a  few  iterations  on  subtracting  the  ion  current,  the  plasma  potential  usually 
converges  to  a  single  value  and  has  little  affect  on  changing  the  value  of  the  density. 
However,  there  will  be  a  much  more  significant  effect  on  using  the  plasma  potential  to 
obtain  the  electron  temperature.  The  electron  temperature  is  found  from  taking  the 
inverse  slope  of  a  line  drawn  through  the  region  between  the  floating  and  plasma 
potentials;  also  shown  in  Figure  D.6.  In  this  case,  the  slope  of  the  line  for  this  plasma 
potential  results  in  an  electron  temperature  of  6.96  eV;  since  this  was  a  single  data  trace, 
the  only  error  considered  was  the  sheath  /  probe  impedance  for  the  collected  saturation 
current.  Utilizing  the  method  described  at  the  end  of  Chapter  4,  the  error  calculated  was 
7.8  %  resulting  in  a  Te  of  6.96  ±  0.54  eV. 

The  next  step  in  calculating  the  density  is  to  insert  the  corresponding  values  with  this 
temperature  into  (D.5).  Using  the  plasma  potential  and  electron  temperature,  along  with 
an  appropriate  value  for  the  ion  saturation  current  at  approximately 

(D.9) 
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where  letting  X*  =  15,  so  that  Vsat  S  -50  Volts  and  rp  /  =  7  so  that  j*=  3,  the  initial 

value  for  the  plasma  density  aeeording  to  (D.5)  is  2.69  x  10  m'  .  The  iterative  method 
then  proeeeds  based  on  the  new  value  for  the  Debye  length  ealeulated  from  this  density 
so  that  the  new  rp  /  =  12  and  7*=  2.8.  After  4  iterations,  7*  ehanges  by  only  a  small 

18  3 

amount  and  the  final  density  is  given  by  n  S  4.04  x  10  m'  . 


An  alternative  method  of  obtaining  a  first  oaleulation  for  the  density  is  to  use  the  slope  in 
the  line  of  the  I  vs.  AV  graph  in  accordance  with  (D.8).  . 


The  slope  can  then  be  used  in  (D.8)  as 


1.42x10 


15 


^ju[amu] 

Aprohei’^^] 


slope 


(D.IO) 


where  msiope  = 


Using  this  method,  msiope  =  5.26  x  10'"^  and  n  =  1.74  x  10^*  m'^. 
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This  method  gives  a  density  (for  this  partieular  ease)  which  is  35%  lower  than  that 
calculated  with  the  iterative  Lafromboise  method.  The  difference  is  explained  by  the  fact 
the  expression  (D.8)  uses  the  limit  for  rp  /  A,d  0  which  is  for  a  very  “thick”  sheath, 
whereas  using  the  Lafromboise  iterative  method,  the  result  for  the  density  is  based  off  a 
converged  value  of  rp  /  A,d  14.43  and  more  accurately  accounts  for  the  effective  ion 

collection  area  as  opposed  to  the  geometric  area  of  the  probe. 

This  method  was  implemented  as  obtained  in  an  EXCEL  workbook  from  Dr.  Martin 
Neumann.  Small  variations  were  made  to  more  efficiently  analyze  IV  characteristics 
such  as  using  the  VLOOKUP  function  in  EXCEL.  However,  all  subjective  user  input 
data  and  graphical  interpretation  lies  with  the  current  author  as  it  was  the  author’s 
responsibility  to  ensure  accurate  data  interpretation  in  the  analysis  program. 
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APPENDIX  E:  ^iWAVE  INTERFEROMETER:  CALIBRATION  -  ANALYSIS 

This  Appendix  will  address  the  analysis  and  derivation  for  obtaining  the  plasma  density 
as  a  function  of  phase  change  measured  by  the  interferometer  signal  which  traverses  the 
plasma  column  of  chord  length  7’.  The  analysis  as  it  is  presented  here  is  only  unique  in 
the  sense  of  how  the  interferometer  is  setup.  It  should  be  widely  applicable  due  to  the 
general  analysis  as  far  as  the  separation  of  both  interferometer  branches  referenced  here 
as:  I-  and  Q-. 


E,1  AFRL ’s  Microwave  Interferometer 

The  general  layout  of  AFRL’s  microwave  interferometer  is  shown  in  Figure  E.l 


Standard  Gain  Conical  Horn 
Antenna 


Q  I 


Figure  E.l,  Schematic  of  the  millimeter  I/Q  Phase  Bridge  Interferometer, 
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The  Gunn  Oscillator  launches  a  90.0  GHz,  17dB  (50mW)  sin  wave  with  a  tuning 
bandwidth  of  8  GHz  and  a  frequency  stability  of  5  MHz.  As  the  wave  passes  through  the 
short  slot  coupler,  it  is  split  into  two  separate  waves  that  are  separated  in  phase  by  90°. 
The  first  wave  or  the  upper  wave  in  Figure  E.l  is  transmitted  at  approximately  13  dB 
through  the  short  slot  coupler  and  phase  shifted  by  90°.  This  is  the  leg  of  the 
interferometer  that  will  pass  through  the  plasma.  After  being  transmitted  across  the 
plasma  length  7’,  the  signal  is  again  shifted  90°  through  a  waveguide  twist  and  separated 
through  another  short  slot  coupler.  This  is  to  ensure  that  the  signal  arriving  at  the  right 
hand  side  of  both  the  “I”  and  “Q”  mixers  in  Figure  E.l  are  at  the  same  phase  angle.  This 
is  due  to  the  90°  shift  from  the  first  short  slot  coupler  and  the  additional  90°  shift  from  the 
waveguide  twist,  i.e.,  ideally  the  phase  change  in  the  wave  on  the  right  hand  side  of  the 
diagram  is  due  to  traversing  the  plasma  and  is  only  out  of  phase  with  respect  to  which  it 
was  originally,  6int. 

The  bottom  leg  of  the  interferometer  is  only  attenuated  in  power  through  the  first  short 
slot  coupler.  No  phase  shift  is  initially  introduced  until  the  signal  is  again  separated 
through  the  second  short  slot  coupler.  The  wave  arriving  at  the  “Q”  mixer  is  shifted  90° 
by  the  second  short  slot  coupler  and  the  wave  arriving  at  the  “I”  mixer  is  again  unaffected 
in  phase  but  attenuated  in  power. 

The  resulting  waves  arriving  at  the  “I”  mixer  both  have  a  zero  degree  phase  shift.  As  a 
result  the  output  voltage  of  this  wave  is  given  by: 

VouT,  =  A  (E.l) 
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where  Voff.i  is  some  offset  voltage  introdueed  by  supplying  power  to  the  mixer  and  ^4/  is 
the  arriving  amplitude  of  the  eombined  waves  due  to  their  power  attenuation  from 
traveling  through  the  short  slot  eouplers  and  aeross  the  plasma. 

Similarly,  for  the  “Q”  mixer,  where  the  arriving  waves  are  90°  out  of  phase,  the  output 
voltage  is  given  by: 

+  +  (E.2) 

or 

VouT^  =  4  +  ^OFF^  (E.3) 

where  Aq  is  again  the  arriving  amplitude  of  the  eombined  waves  and  VoFFgis  the  voltage 
offset  due  to  the  supply  power. 

It  should  be  explained  at  this  point  that  “0”  is  the  phase  ehange  due  to  all  external  faetors 
ineluding  those  ehanges  due  to  traversing  the  plasma.  Other  faetors  or  sourees  that  eould 
introduee  phase  errors  result  from  non-uniformity  in  the  path  lengths  between  the  two 
arms,  vibrations  in  the  room,  temperature  gradients  throughout  the  waveguide,  and  any 
medium  other  than  plasma  between  the  eonieal  gain  horns  sueh  as  a  quartz  tube  or  other 
obstruetion.  The  motivation  is  now  to  isolate  the  phase  ehange  that  is  only  due  to  the 
plasma.  Onee  this  is  obtained,  then  following  a  suitable  analysis,  the  radially  integrated 
eleetron  density  ean  be  obtained.  However,  beeause  is  due  to  all  the  sourees  of 
phase  change, 

^  ^plasma  ^ M^athLength  ^Vibrations  ^etc .  (E.4) 
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we  can  rewrite  this  as, 


O  =  (E.5) 

where  6o  is  the  phase  error  due  to  all  external  sources  which  can  be  obtained  through 
proper  calibration. 


At  this  point  we  can  now  solve  (E.l)  and  (E.3)  for  the  total  phase  change  “0”.  (E.l) 
becomes, 

sin(^)  =  ^OFF, 

4 


and  (E.3)  becomes, 

(E.7) 

Dividing  (E.6)  by  (E.7)  we  obtain, 

tan(0)  =  ^^E17^J^  =  ^  (E.8) 

AVour,-Vo,,, 

We  can  now  rewrite  the  total  phase  shift  in  terms  of  its  component  parts  as  given  by  (E.5) 
and  separate  out  tan(O), 

^^^i0 plasma^  0q) 

Here,  6i  and  6q  are  the  phase  errors  introduced  to  each  mixer  “I”  and  “Q”  respectively, 
not  inclusive  of  the  phase  change  due  to  the  plasma. 


Expanding  the  sin  and  cos  in  terms  of  their  summation  identities. 
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plasma  )  COS(0,  )  +  )  Sin(0,  ) 


(E.IO) 


+^e)  =  cos(0^,^,„Jcos(^g)-sin(0^,^,„Jsin(0g)  (E.ll) 

Inserting  (E.IO)  and  (E.l  1)  into  equation  (E.9)  we  obtain, 

plasma  )  I  )  +  plasma  )  ^^0  1 )  _ 

plasma )  )  sin(0g  )  (E.  1 2) 


Multiplying  the  left  hand  side  of  (E.12)  by 


1/ 

plasma 

1/ 

plasma 


) 

) 


we  obtain, 


tan(^p/..^Jcos(^,)  +  sin(^,) 
oos(6'g  )  -  tan(0^,^,^^ )  sin(0g  ) 


(E.13) 


Rearranging  (E.13)  we  have, 

plasma )  cos(0, )  +  sin(0, )  =  ^  eos(0g  )  -  P  tan(0^,^,„^ )  sin(0g  )  (E.  1 4) 


Einally,  solving  in  terms  of  ian{6piasma),  we  obtain. 


plasma)  = 


y5cos(^g)-sin(^^) 
oos(^/)  +  y0sin(^g) 


and  the  phase  ehange  due  to  the  plasma  is  given  by. 


^plasma  = 


y0eos(^Q)-sin(0^) 

1 - - - r 

eos(0^)  +  y0sin(^(5) 


(E.15) 


(E.16) 


Therefore,  having  an  expression  for  the  phase  ehange  due  to  the  plasma,  the  remainder  is 
left  to  ealibration  and  measurement.  Measurement  will  come  in  the  form  of  a  voltage 
deflection  due  to  the  presence  of  plasma  or  any  other  material  between  the  conical  gain 
horns,  so  at  this  point  it  is  a  matter  of  calibrating  the  phase  errors  and  amplitude 
differences  between  both  I-  and  Q-  branches. 
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E.2  Calibration  and  Analysis 

Calibration  begins  with  voltage  measurement  of  the  I-  and  Q-  signals  when  there  is  no 
plasma  present.  In  the  current  calibration,  the  quartz  tube  is  within  the  conical  gain  horns 
in  order  to  get  a  baseline  calibration  so  that  later  the  only  change  between  the  horns  will 
be  the  presence  of  Argon  plasma.  Calibration  encompasses  inserting  a  Mica  sheet  within 
the  ‘top’  branch  of  the  waveguide.  By  slowly  moving  a  well-known  material 
incrementally  further  into  the  waveguide,  this  causes  the  traversing  wave  to  be  phase 
shifted,  i.e.,  simulating  a  plasma  deflection  although  in  this  case,  we  are  performing  a 
calibration  over  a  full  360  degree  rotation  accomplished  by  moving  the  Mica  sheet.  The 
raw  data  signal  of  performing  this  is  shown  in  Figure  E.2  where  the  vertical  line 
represents  the  location  where  data  will  be  taken,  i.e,  after  the  Mica  sheet  has  induced  a 
full  360  degree  wave  phase  shift,  the  Mica  was  then  reinserted  to  a  known  ‘phase 
location’.  This  location  was  selected  as  a  point  where  both  branches  were  nearly  90 
degrees  out  of  phase  with  respect  to  each  other  so  as  to  minimize  large  phase  errors  that 
could  otherwise  be  introduced. 
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Combined  Phase  Shifter  Position  (|u,m) 
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Figure  E.2.  Raw  signals  from  both  I-  and  Q-  branch  as  the  Mica  sheet  is 
incrementally  inserted  further  into  the  waveguide  causing  a  phase  shift. 


The  next  step  is  to  calculate  the  signal  amplitude.  This  is  most  easily  accomplished  by 
removing  the  average  voltage  offset  and  V^pp^  then  taking  V2  the  difference  between 

the  maximum  and  minimum  values.  In  this  case,  the  wave  amplitudes  and  offsets  for  use 
in  (E.8)  for  calculation  of  (3  are  Ai  =  21.3  mV,  Aq  =  8.7  mV,  V^pp^  =  410.34  mV,  and 


Vqpp^  =  602.50  mV.  The  amplitudes  are  not  equal  because  both  signals  are  attenuated 

differently  as  they  traverse  their  respective  path  lengths  and  short  slot  couplers  while  the 
voltage  offsets  are  also  not  equal  due  to  small  differences  in  the  mixer  electronics,  most 
likely  circuit  impedances. 


The  final  step  is  to  calibrate  the  phase  shifts  in  each  branch  which  are  the  result  of  path 
length  differences,  presence  of  quartz  tube,  thermal  effects,  mechanical  vibrations,  and 
any  other  such  external  influences.  In  the  ideal  case  where  there  would  be  no  Mica  sheet 
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for  calibration,  data  would  be  taken  at  the  zero  location  as  shown  with  the  normalized 


data  in  Figure  E.3. 


Combined  Phase  Shifter  Position  (|um) 


Figure  E.3.  Normalized  and  horizontally  shifted  I-  and  Q-  branches  as  the  Mica  sheet  is 

traversed  across  the  waveguide. 

Clearly,  these  waves  at  the  zero  location  are  not  pure  sin  or  cos  waves  as  would  be 
expected.  In  fact,  at  the  zero  location  they  are  significantly  shifted  from  their  ideal 
behavior.  The  I-branch  {sin)  does  not  start  to  behave  like  a  {sin)  wave  until  the  Mica 
sheet  is  inserted  approximately  680  |um  into  the  waveguide  while  the  Q-branch  {cos)  does 
not  behave  ideally  until  approximately  the  same  location.  Therefore,  it  was  decided  to 
move  the  Mica  sheet  to  a  location  (fairly  arbitrary  yet  easily  quantifiable)  in  the  region 
where  the  wave  behavior  more  nearly  resembles  the  ideal  behavior.  To  do  this,  both 
normalized  waves  were  2ti  scaled  and  the  location  corresponding  to  the  Mica  sheet 
insertion  was  obtained.  The  result  is  shown  in  Figure  E.4. 
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Figure  E.4.  Normalized  and  In  scaled  location  of  the  Mica  sheet  insertion  for  the 
point  where  interferometry  data  was  taken. 

The  resulting  shift  in  the  I-braneh  (sin)  was  0i  =  62.08  degrees  while  the  Q-branch  (cos) 
0Q  =  60.79  degrees. 

With  the  calibration  values  now  obtained  (Aj,  Aq,  0i,  and  0q)  for  use  in  (E.16),  a  plasma 
deflection  will  induce  a  voltage  difference  so  that  Opiasma  can  be  calculated.  As  an 
example,  consider  the  voltage  deflection  due  to  500  Watt  inductive  plasma  as  shown  in 
Figure  E.5.  The  average  voltage  deflections  are  measured  as  =  2.31  mV  and 

^^Qavg  ~  above  calibration  values,  (E.16)  can  then  be  used  to 

obtain  =7.88  degrees  =  0.1375  radians  .  This  can  then  be  used  to  calculate  the 
chord-averaged  plasma  density  in  the  usual  manner  i.e.. 
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-  An:m„SnC^ 


He  = 


plasma 


AJ 


(E.17) 


where  is  the  wavelength  corresponding  the  90  GHz  frequency  ~  3.33  mm  and  7’  is  the 
path  length  of  the  plasma,  in  this  case  ~  5  cm.  Substituting  the  values  yields  a  plasma 
density  of  ~  2.93  x  10^^  m'^. 


Plasma  Voltage 
Deflection 


607.5 

607.4 

607.3 

607.2 

£ 

607.1 

U 

a 

607.0 

CQ 

.fi 

606.9 

O' 

606.8 

606.7 

606.6 

606.5 

-Q  -  branch 
-I-branch 


E 


427.0  g 

CQ 
u 
.fi 


=0.42  mF 

^AVG 


AH  =2.31  mF 

^AVG 


“I - 1 - 1 - 1 - 1 - 1 - 1 - r 


0  10  20  30  40  50  60  70  80  90  100  110  120 

Time  [s] 

Figure  E.5.  Voltage  deflections  in  the  I-  and  Q-  branches  of  the  microwave  interferometer  used 

to  obtain  p. 
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APPENDIX  F:  RADIAL  bz  FIELDS  -  PRESSURE,  POWER,  MAGNETIC  FIELD 


This  Appendix  compares  the  differences  between  the  measured  radial  bz  profiles  and 
those  predicted  by  considering  the  numerical  solutions  to  a  non-uniform  radial  density 
profile.  In  each  case  shown  in  Figure  F.l  -  Figure  F.6,  the  GREEN  border  cases 
demonstrate  strong  agreement  with  the  predicted  fields.  The  “w”  value  for  each 
predicted  profile  refers  to  width  of  the  input  Gaussian  density  profile.  YELLOW  border 
cases  indicate  what  would  be  interpreted  as  a  helicon  wave  profile  had  the  probe  not 
interfered  with  the  wave  as  the  probe  was  radially  scanned.  All  scans  took  place 
underneath  the  end  ring  of  the  “downstream”  antenna.  RED  border  cases  illustrate  a  flat 
profile  where  power  (density)  and  /  or  magnetic  field  conditions  were  insufficient  for  the 
medium  to  support  helicon  wave  propagation.  Additional  conditions  were  considered  for 
the  case  of  0.5  mTorr  primarily  as  an  investigation  into  reported  helicon  double  layers. 
All  conditions  investigated  at  0.5  mTorr  demonstrated  no  wave  propagation  indicative  of 
the  increased  difficulty  in  propagation  of  helicon  waves  at  lower  pressures;  difficulty  in 
terms  of  power  and  /  or  magnetic  field  requirements  as  compared  to  higher  pressure 
cases.  Not  all  conditions  within  each  power  -  magnetic  field  matrix  were  investigated 
after  the  determination  of  “no-wave  propagation”  below  a  minimum  power  or  magnetic 
field  was  made. 
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300  Gauss  600  Gauss  900  Gauss  300  Gauss  600  Gauss  900  Gauss 


500  Watts  400  Watts  300  Watts 


Figure  F.l.  profiles  at  10  mTorr  and  various  input  power  and  magnetic  field  strengths. 


500  Watts  400  Watts  300  Watts 


Figure  F.2.  profiles  at  8  mTorr  and  various  input  power  and  magnetic  field  strengths. 
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300  Gauss  600  Gauss  900  Gauss  300  Gauss  600  Gauss  900  Gauss 


500  Watts  400  Watts  300  Watts 


Figure  F.3.  profiles  at  6  mTorr  and  various  input  power  and  magnetic  field  strengths. 


500  Watts  400  Watts  300  Watts 


Figure  F.4.  profiles  at  4  mTorr  and  various  input  power  and  magnetic  field  strengths. 


198 


0  Gauss  300  Gauss  600  Gauss  900  Gauss  300  GauSS  600  GauSS  900  GauSS 


500  Watts  400  Watts  300  Watts 


Figure  F.5.  profiles  at  2  mTorr  and  various  input  power  and  magnetic  field  strengths. 
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Figure  F.6.  bj  profiles  at  0.5  mTorr  and  various  input  power  and  magnetic  field  strengths. 
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